Fraktale - Regelmäßigkeit im Chaos

6. Ein Pascalprogramm für die Mandelbrotmenge

program fraktal;
uses graph,crt;
const x1      :real=-2;         {kleinster x-Wert}
      y1      :real=-1;         {kleinster y-Wert}
      wt      :real=3;          {x-Schrittweite}
      ht      :real=2;          {y-Schrittweite}
      maxx    :word=500;        {x-Breite der Zahlenebene in Pixel}
      maxy    :word=400;        {y-Länge der Zahlenebene in Pixel}
var dx,dy,x,y                 :real;
    i,j,col                   :word;{col=Farbennr;i,j=Laufvariablen}
    a,b                       :integer;
function mandeltest(x,y:real;w:word):word;
var re,im        :real;
    re2,im2      :real;
    i            :word;
begin
  re:=x;
  im:=y;
  for i:=0 to 2+w do begin
    re2:=re*re;
    im2:=im*im;
    if re2+im2>256 then begin
      mandeltest:=i+1;
      break;
    end;
    im:=2*re*im+y;{bei Juliamengen fuer y imaginaeren Teil von c einsetzen}
    re:=re2-im2+x;{bei Juliamengen fuer x realen Teil von c einsetzen}
  end;
  mandeltest:=i+1;
end;
begin
  a:=detect;
  b:=0;
  initgraph(a,b,'');
  dx:=wt/maxx;
  dy:=ht/maxy;
  for i:=0 to maxx do
    for j:=0 to maxy do begin
      x:=x1+dx*i;
      y:=y1+dy*j;
      col:=mandeltest(x,y,300);
      if not(col>=300) then putpixel(i,j,col);
    end;
  readkey;
end.