%script file to make of movie of filled julia sets 
%varying w/ a parameter

%create real and imaginary value arrays
[X,Y]=meshgrid([-2:0.01:2]);
[m,n]=size(X);
maxiters=30;
clear i;
%vary your parameters as you see fit, fix one and vary the other
b=0+i*0;
ar=0:0.2:3;
for p=1:length(ar)
    a=-ar(p)+i*0;
    escrit=[abs(b),sqrt(abs(a)+2)];     %escape criteria must be
    escriter=max(escrit);              %calculated each time through
    T=maxiters*ones(m,n);
    for j=1:m
        fprintf ('on iterate %d\n',j);
        for k=1:n
            z=X(j,k)+i*Y(j,k);
            for t=1:maxiters
                y=z^3+a*z+b;
                z=y;
                if abs(z)>escriter
                T(j,k)=t;
                break;
                end
            end
        end
    end
imagesc(T);
pause(0.1);
%Need to get T as a movie frame
M(:,p)=getframe;
end
movie(M,3);