% Dieses Programm visualisiert eine Julia-Menge 
% (vgl. Lehrbuch, Kap. 9.3).
%
% Die Festlegung des Bildausschnitts erfolgt durch
% Angabe der folgenden Parameter:
% ux = kleinste x-Koordinate
% ox = größte x-Koordinate
% nx = Anzahl der Gitterpunkte auf der x-Achse
% uy = kleinste y-Koordinate
% oy = größte y-Koordinate
% ny = Anzahl der Gitterpunkte auf der y-Achse
%
% Der Parameter kmax legt die Anzahl der Iterationen fest,
% der Parameter c hat dieselbe Bedeutung wie im Buch.

ux = -1.8;
ox = 1.8;
nx = 400;
x = linspace(ux,ox,nx);
uy = -1.15;
oy = 1.15;
ny = 255;
y = linspace(uy,oy,ny);
[xx,yy] = meshgrid(x,y);
kmax = 80;
w = xx+sqrt(-1)*yy;
c = -1.25
% c = -0.75
% c = 0.35+0.35*sqrt(-1)
% c = -0.03+0.655*sqrt(-1)
% c = -0.12+0.74*sqrt(-1)
% c = 0.365-0.3*sqrt(-1)
for k = 1:kmax
  k
  w = w.^2+c;
end
w = min(100,abs(w));
colormap(gray)
image([ux ox], [uy oy], w)
axis equal