RALP3[a_,b_,c_,d_,n_,s_,col_]:=Module[{rt,a1,b1,c1,d1},
a1=.85+.001a;b1=.9+.01b;c1=.4+.01c;d1=9+.1d;
rt=RecurrenceTable[{x[i+1]==
a1+b1(x[i]Cos[c1-d1/(1+x[i]^2+y[i]^2)]-y[i]Sin[c1-d1/(1+x[i]^2+y[i]^2)]),
y[i+1]==b1(x[i]Sin[c1-d1/(1+x[i]^2+y[i]^2)]+y[i]Cos[c1-d1/(1+x[i]^2+y[i]^2)]),
x[0]==.4,y[0]==.5},{x,y},{i,1,n}];
ListPlot[rt,ColorFunction->Function[{x,y},col[Cos[16x+25y]^2]],
Axes->False,Background->RGBColor["black"],
ImageSize->s,AspectRatio->1]];
{RALP3[0,0,0,0,25000,400,Hue],
RALP3[-72,10,75,-50,15000,400,ColorData["BrightBands"]]}
{RALP3[0,0,0,0,25000,400,Hue],
RALP3[-72,10,75,-50,15000,400,ColorData["BrightBands"]]}
No comments:
Post a Comment