%% Fractals
clc; clear all; close all
Inc=.0025
[xMat,yMat]=meshgrid(linspace(.3+Inc,.3+2*Inc,500),linspace(.025,.025+Inc,500));
zMat=zeros(size(xMat));
X=xMat;
Y=yMat;
m=1

while m <1000
   xTemp=xMat.^2-yMat.^(2)+X;
   yTemp=2*xMat.*yMat+Y;
   z=sqrt(xTemp.^2+yTemp.^2);
   xMat=xTemp;
   yMat=yTemp;
   z(find(zMat>0))=0;
   zMat(find(z>=4))=m;
   m=m+1;
   imagesc((log(zMat)))
   pause(.0000000000005)
end

%% good settings
%[xMat,yMat]=meshgrid(linspace(.3,.3125,500),linspace(.025,.0175,500));
%Inc=.005
%[xMat,yMat]=meshgrid(linspace(.3,.3+Inc,500),linspace(.025,.025+Inc,500));
%
%
%
%