Source Codes zum Feigenbaumdiagramm

program GrafischeIteration_Gerade;

uses
  crt,
  graph;

const
  xl=0; xr=10; yu=0; yo=10;

var
  xmax, ymax : word;
  mx, my     : real;
  i, xb, yb,
  xbo, ybo   : integer;
  anz        : longint;
  ca, xa, x0,
  x, y, c    : real;
  auswahl,
  wahl, ok   : char;

procedure eingabe;
begin
  repeat
    clrscr;
    writeln('Welche Funktion der Art ');
    writeln;
    writeln('f(x) = m*x+c    soll untersucht werden ? ');
    writeln;
    writeln;
    write('Bitte Koeffizient m eingeben : ');readln(ca);
    write('Bitte Koeffizient c eingeben : ');readln(c);
    writeln;
    write('Anfangswert der Iteration: ');readln(x0);
    write('Anzahl der Iterationen: ');readln(anz);
    mx:=639/(xr-xl);
    my:=479/(yo-yu);
    ymax:=479;
    xmax:=639;
    writeln;
    write('Ist die Eingabe OK (j/n) ? : ');
    ok:=readkey;
  until ok='j';
end;

function fvon(x,m,c:real):real;
begin
  fvon:=m*x+c;
end;

procedure grafik_ein;
var
  gt,gm:integer;
begin
  gt:=detect;
  initgraph(gt,gm,'c:\bp\bgi');
end;

procedure Koordinatensystem_zeichnen;
begin
  xbo:=round(-xl*mx);
  ybo:=round(yo*my);
  line(0,ybo,xmax,ybo);
  line(xbo,0,xbo,ymax);
end;

procedure Schaubild_zeichnen;
begin
  for xb:=0 to xmax do
    begin
    x:=xl+xb/mx; {x-links plus xb im Ma~stab}
    y:=fvon(x,ca,c);
    yb:=round((yo-y)*my);
     {Ursprung minus (da Achse andersherum) y-Wert im Maßstab}
    putpixel(xb,yb,4);
    x:=xl + xb/mx;
    y:=x;              {1. winkelhalbierende }
    yb:=round((yo-y)*my);
    putpixel(xb,yb,3);
    end;
  xa:=x0;moveto(round(x0*mx)+xbo,ybo);
  for i:=1 to anz do
    begin
    xa:=ca*xa+c;
    yb:=ybo;
    lineto(round(x0*mx)+xbo,ybo-round(xa*my));
    lineto(xbo+round(xa*mx),ybo-round(xa*my));
    x0:=xa;
    end;
end;

begin
  eingabe;
  grafik_ein;
  Koordinatensystem_zeichnen;
  schaubild_zeichnen;
  readln;
end.