ex4_21.m

源代码集合

x=linspace(-1.5,1.5,300); y=x; [X,Y]=meshgrid(x,y);
n_iter=30; c=-0.7;
tic, W=julia(X,Y,c,n_iter); toc

pcolor(X,Y,W), shading flat; axis('square'); colormap prism(256);

n_iter=300; tic, W1=julia(X,Y,c,n_iter); toc
pcolor(X,Y,W1), shading flat; axis('square'); colormap prism(256);

max(max(abs(W-W1)))

n_iter=400; tic, W2=julia(X,Y,c,n_iter); toc
pcolor(X,Y,W2), shading flat; axis('square'); colormap prism(256);

max(max(abs(W2-W1)))

colormap flag(256);

n_iter=300; c=-1+0.1j; tic, W1=julia(X,Y,c,n_iter); toc
pcolor(X,Y,W1), shading flat; axis('square'); colormap prism(256);