h = figure;
set(h, 'Position', [0 0 1000 1000])

center = [-.73659007459575,-.69270903692418];
%frame = 43
for frame = 1:.1:43
     scale = (1/2^frame)
     [X,Y] = meshgrid(linspace(center(1)  - scale,center(1) + scale,100),linspace(center(2) - scale,center(2) + scale,100));
     Z = ones(size(X));

    for x_ = 1: size(X,1)
        for y_ = 1: size(X,2)
            u = [X(x_,y_),Y(x_,y_)];
            for i = 1:1000
    %             julian set
%                   v = [u(1) * u(1) - u(2) * u(2) +  frame, u(2) * u(1) + u(2) * u(1)+ frame];

    %           mandelbrot set
%                    v = [(u(1) * u(1) - u(2) * u(2)) +  X(x_,y_) , (u(2) * u(1) + u(2) * u(1)) + Y(x_,y_)];
                   v = [(u(1) * u(1) - u(2) * u(2)) + .3, (u(2) * u(1) + u(2) * u(1)) ];

                if((v(1) * v(1) + v(2) * v(2)) > 4)
                    break
                end
                u = v;
                Z(x_,y_) = i;
            end
        end
    end
    fx = X(:);
    fy = Y(:);
    imagesc(fx,fy,Z)
    set(gca,'xtick',[])
    set(gca,'ytick',[])
    drawnow


      result = getframe(h);
      im = frame2im(result);
      [imind,cm] = rgb2ind(im,256);

    if frame == 1
        imwrite(imind,cm,'final.gif','gif', 'Loopcount',inf);
    else
         imwrite(imind,cm,'final.gif','gif','DelayTime',0.1,'WriteMode','append');
    end

end