Fern

Joe Riel — Fri, 27 Nov 2009 00:53:16 Z

Fern := proc(n::posint
             , x::Array(datatype=float[8])
             , y::Array(datatype=float[8])
             , r::Array(datatype=float[8])
            )
local i, pick;
option autocompile;
    for i to n-1 do
        pick := r[i];
        if pick < 0.01 then
            x[i+1] := 0;
            y[i+1] := 0.16*y[i];
        elif pick < 0.08 then
            x[i+1] := 0.2*x[i] - 0.26*y[i];
            y[i+1] := 0.23*x[i] + 0.22*y[i] + 1.6;
        elif pick < 0.15 then
            x[i+1] := -0.15*x[i] + 0.28*y[i];
            y[i+1] := 0.26*x[i] + 0.24*y[i] + 0.44;
        else
            x[i+1] := 0.85*x[i] + 0.04*y[i];
            y[i+1] := -0.04*x[i] + 0.85*y[i] + 1.6;
        end if;
    end do;
    NULL;
end proc:

n := 10^5:
r := rtable(1..n, frandom(0..1,1), 'datatype'=float[8]):
x := Array(1..n, 'datatype'=float[8]):
y := Array(1..n, 'datatype'=float[8]):

(Fern(n,x,y,r));
plot(x,y, 'style=point', 'color=green');

Thanks very much! Also gives

Deadstar — Fri, 27 Nov 2009 01:15:42 Z

Thanks very much!

Also gives me a basis to create my own ones.

MaplePrimes - answers and comments on Question, Fractal Fern question

Fern

Thanks very much! Also gives