%6.b. julia sets:
function M=julia(zMax, c, N)
m=500;
[Re,Im]=meshgrid(linspace(-zMax, zMax, m),linspace(-zMax, zMax, m));
Z=Re+1i*Im;

M=zeros(m);
for a=1:m 
    for b=1:m
M(a,b)=escapeVelocity(Z(a,b),c,N);
    end
end

imagesc(-zMax:zMax, -zMax:zMax, atan(0.1*M));
%furán mûködik a skálázás!!!
xlabel('Re(Z)');
ylabel('Im(Z)');
axis xy;
end

%julia (1, -.297491+i*0.641051, 250);
%julia (.5, -.297491+i*0.641051, 250);