M=1.5;
N=1.5;
limit=2;
X=linspace(-M,M,1600);
Y=linspace(-N,N,1600);
Maxcount=64;
count=1;
for ii=X
for jj=Y
c=ii+jj*i;
z=0;
kk=0;
while abs(z)<limit && kk~=Maxcount
z=sin(z^2+conj(2*z))+sin(c);
++kk;
endwhile
F(count++)=1-kk/Maxcount;
endfor
endfor
F=reshape(F,[length(X) length(Y)]);
imwrite(F,'Mandelbrotvariant.tiff','tiff');
కామెంట్లు లేవు:
కామెంట్ను పోస్ట్ చేయండి