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