「シダの葉グラフィクス」

 2017-10-20 には、CASIO fx-CG50が国内発売される、との事で、カラーグラフィクスの表示できる高機能電卓向けに、「シダの葉グラフィクス」のプログラムを実行してみました。

シダの葉グラフィクスについては、つぎを御参照。

反復関数系 - wikipedia
https://ja.wikipedia.org/wiki/%E5%8F%8D%E5%BE%A9%E9%96%A2%E6%95%B0%E7%B3%BB


当方、未だカラーグラフィクスの高機能電卓を持っていないので、HP Primeのシミュレータで動作するコードを作成してみました。
以下の版は、モノクログラフィクスでの表示ですが、RGB()のパラメタを少し変えるだけで、カラーになります。

プログラム「draw_shida.ppl」

//  title : draw_shida.ppl
//  begin : 2017-10-03 15:18:44 
//  

//  
width  := 320;
height := 240;

//  initial ooord
x := 0;
y := 0;
//  loop counte
cnt :=  30000;
//  working
tmp := 0;

xmin := -3;
xmax :=  5;
ymin :=  0;
ymax :=  8;

//  functions

//  f1()
f1()
BEGIN
  tmp := 0;
  y := 0.16*y;
  x := tmp;
END;

//  f1()
f2()
BEGIN
  tmp := 0.2*x  - 0.26*y;
  y := 0.23*x + 0.22*y + 1.6;
  x := tmp;
END;

//  f1()
f3()
BEGIN
  tmp := -0.15*x + 0.28*y;
  y :=  0.26*x + 0.24*y + 0.44;
  x := tmp;
END;

//  f1()
f4()
BEGIN
  tmp :=  0.85*x + 0.04*y;
  y :=  -0.04*x + 0.85*y + 1.6;
  x := tmp;
END;


//  MAIN 
EXPORT DRAW_SHIDA()
BEGIN
// Clean the screen (G0):
RECT();

  LOCAL  i, r, col;
  LOCAL  xp, yp;

  FOR i FROM 1 TO cnt DO
    r := random();
    IF (r < 0.01)               THEN f1(); col := RGB(0, 0, 0); END;
    IF (r >= 0.01 AND r < 0.07) THEN f2(); col := RGB(0, 0, 0); END;
    IF (r >= 0.07 AND r < 0.15) THEN f3(); col := RGB(0, 0, 0); END;
    IF (r >= 0.15)              THEN f4(); col := RGB(0, 0, 0); END;
    xp := (x-xmin)/(xmax-xmin)*width;
    yp := (y-ymin)/(ymax-ymin)*height;
    PIXON_P(yp, xp, col);
  END;


//  key wait loop
REPEAT UNTIL GETKEY() == -1;
FREEZE;

END; 

小汚いコードですが、 移植は簡単だと思います。