%[Z,X,Y]=makegrid(-2+2*i,2-2*i,300);
%Note kludge reversal of imag parts, imwrite reverses top/bottom

%imwrite(showcplxm(X,Y,Z),'plane.png','png','Description','z, adj1');
%imwrite(showcplxm(X,Y,Z,0.5),'planep5.png','png','Description','z, adj.5');
%imwrite(showcplxm(X,Y,Z,5),'plane2.png','png','Description','z, adj3');

%Z=Z.*(Z+1).^2./(Z-1);
%imwrite(showcplxm(X,Y,Z,0.65),'ratfn.png','png','Description','z(z+1)^2/(z-1), adj.65');

[Z,X,Y]=makegrid(-2+2*i,2-2*i,300);
W=Z.^3-1;
imwrite(showcplxm(X,Y,W,0.65),'cubic.png','png','Description','z^3-1, adj.65');
%Z=Z-(Z.^3-1)./(3*Z.^2);
%imwrite(showcplxm(X,Y,Z,0.4),'ncubic1.png','png','Description','N1:z^3-1,adj.4');
%Z=Z-(Z.^3-1)./(3*Z.^2);
%imwrite(showcplxm(X,Y,Z,0.4),'ncubic2.png','png','Description','N2:z^3-1,adj.4');
%for k=3:5, Z=Z-(Z.^3-1)./(3*Z.^2); end;
%imwrite(showcplxm(X,Y,Z,0.4),'ncubic5.png','png','Description','N5:z^3-1,adj.4');
%for k=6:20, Z=Z-(Z.^3-1)./(3*Z.^2); end;
%imwrite(showcplxm(X,Y,Z,0.4),'ncubic20.png','png','Description','N20:z^3-1,adj.4');

[Z,X,Y]=makegrid(-2+2*i,2-2*i,300);
Z=Z.^2+i;
imwrite(showcplxm(X,Y,Z,0.1),'juld1.png','png','Description','z^2+i,adj.1');
for k=2:3, Z=Z.^2+i; end;
imwrite(showcplxm(X,Y,Z,0.1),'juld3.png','png','Description','I3:z^2+i,adj.1');
for k=4:6, Z=Z.^2+i; end;
imwrite(showcplxm(X,Y,Z,0.1),'juld6.png','png','Description','I9:z^2+i,adj.1');
for k=7:9, Z=Z.^2+i; end;
imwrite(showcplxm(X,Y,Z,0.1),'juld9.png','png','Description','I9:z^2+i,adj.1');
for k=10:12, Z=Z.^2+i; end;
imwrite(showcplxm(X,Y,Z,0.1),'juld12.png','png','Description','I9:z^2+i,adj.1');

[Z,X,Y]=makegrid(-2+2*i,2-2*i,300);
Z=Z.^2+0.3*i;
imwrite(showcplxm(X,Y,Z,0.1),'julc1.png','png','Description','z^2+.3i,adj.1');
for k=2:3, Z=Z.^2+0.3*i; end;
imwrite(showcplxm(X,Y,Z,0.1),'julc3.png','png','Description','I3:z^2+.3i,adj.1');
for k=4:6, Z=Z.^2+0.3*i; end;
imwrite(showcplxm(X,Y,Z,0.1),'julc6.png','png','Description','I3:z^2+.3i,adj.1');
for k=7:9, Z=Z.^2+0.3*i; end;
imwrite(showcplxm(X,Y,Z,0.1),'julc9.png','png','Description','I9:z^2+.3i,adj.1');