Thursday, December 12, 2019

RecurrenceTable Examples 3


CRTAP[b_,c1_,c2_,n_,s_,col_]:=Module[{c,init,rt,data},
c=-0.122561^c1+0.744862^c2*I;
init=Table[x+I*y,{x,-b,b,2b/249},{y,-b,b, 2b/249}];
rt=RecurrenceTable[{z[k]==z[k-1]^3+c,z[0]==init},z,{k,n}];
data=Map[If[Max[Abs[#]]<=2,1,0]&,Flatten[rt,{{2},{3},{1}}],{2}];
MatrixPlot[data,ColorFunction->col,ImageSize->s]];
{CRTAP[1.4,3,-1,5,300,"CoffeeTones"],
CRTAP[1.4,1,1,25,300,"BrightBands"]}

No comments:

Post a Comment