%6.a. julia sets: escape velocity
function n=escapeVelocity(z0, c, N)
%zn=zn-1^2+c
z=zeros(1, N);
z(1, 1)=z0;
n=0;
if abs((z0^2)+c)<=2
  a=2;
  nagyobb=0;
  while a<=N  && nagyobb==0;
    z(1,a)=(z(1,a-1)^2)+c;
    if abs(z(1,a))>2
      nagyobb=1;
      n=a;
    end
    a=a+1;
  end
  if nagyobb==0
      n=N;
  end
end
end