function Y=julia(X)
x=X(1); y=X(2);
z=x+1i*y;
c=1i-0.08375;
f=0.5*(z-c)^2;
Y(1)=real(f);
Y(2)=imag(f);