
axis([-1.5 1.5 -1.5 1.5]);
axis square
hold on;
set(gca,'Fontsize',12)
xlabel('\it x','Fontsize',12)
ylabel('\it y','Fontsize',12)
title('A simple sketch of the Julia set of Figure 11.15','Fontsize',12)

n=500;
a=linspace(-1.5,1.5,n);
b=linspace(-1.5,1.5,n);

for i=1:n
for j=1:n
    
x=a(i);
y=b(j);

for k=1:20
f1=x^2-y^2-.75;
f2=2*x*y;
x=f1;
y=f2;
end
if sqrt(x^2+y^2)<2
plot(a(i),b(j),'erasemode','none');
drawnow;
end
end %end for b
end %end for a
