1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:

procedure ....;
var
Infizierte, Geheilte, Gesunde, Kranke, Kontakthaeufigkeit,
Infektionsrate, Heilungsrate, Immune : Real;
begin

Gesunde:= 1000;
Kranke:=40;
Immune:= 0;
Kontakthaeufigkeit:=0.2;
Infektionsrate:=0.01;
Heilungsrate:=0.1;
x:=0;  // gibt die Tage an, vom 0 - Nullten TAag bis zum 60 Tag praktisch.
   REPEAT
   Infizierte := Gesunde*Kranke*Kontakthaeufigkeit*Infektionsrate;
   Geheilte := Kranke * Heilungsrate;
  // Gesunde.neu:= Gesunde.alt + dt*(-Infizierte)                           
  //Kranke.neu:= Kranke.alt + dt*(Infizierte-Geheilte)
  // Immune.neu:= Immune.alt + dt*(Geheilte)


  ^^ Das ist der Teil den ich nicht hinkriege?!

   //  y:=x;     <-- hier muss dann irgendwie die funktion zum Pixel 
                     übertragen werden. 
    xb:=round(x0+x*k);
    yb:=round(y0-y*k);    // hier Y- Wert ändern
    Bild.Canvas.Pixels[xb,yb]:=farbe;
    x:=x+0.001; // damit es ein Pixel weiter rückt
  UNTIL x > 60 // 60 Tage


end;					

1:
2:
3:
4:
5:

Gesunde:= round(Gesunde.NEU - x*Infizierte);
Geheilte := round(Geheilte.NEU * Heilungsrate);
Kranke:= round(Kranke.NEU + x*(Infizierte-Geheilte));
Immune:= round(Immune.NEU + x*(Geheilte));
Infizierte := round(Gesunde.NEU*Kranke*Kontakthaeufigkeit*Infektionsrate);					

1:
2:
3:
4:

Immune:= 0;
Kontakthaeufigkeit:=0.2;
Infektionsrate:=0.01;
Heilungsrate:=0.1;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:
37:
38:
39:
40:
41:
42:
43:
44:
45:
46:
47:
48:
49:
50:
51:
52:
53:
54:
55:
56:
57:
58:
59:
60:
61:
62:
63:
64:
65:
66:
67:
68:
69:
70:
71:
72:
73:
74:
75:
76:
77:
78:
79:

type
 TData = record
  S, I, R: Extended;
 end;

function f(u: TData): TData;
const
 beta = ...
 mu = ...
 lambda = ...
 gamma = ...
begin
 with u do
 begin
  Result.S := beta -mu*S -lambda*S*I;           // susceptible
  Result.I :=      -mu*I +lambda*S*I -gamma*I;  // infected
  Result.R :=      -mu*R             +gamma*I;  // resistant
 end;
end;

function Mult(u: TData; factor: Extended): TData;
begin
 with u do
 begin
  Result.S := S*factor;
  Result.I := I*factor;
  Result.R := R*factor;
 end;
end;

function Add(u, v: TData): TData;
begin
 Result.S := u.S+v.S;
 Result.I := u.I+v.I;
 Result.R := u.R+v.R;
end;

function RKM4(u: TData; dt: Extended): TData;
var
 v, du: TData;
begin
 du := f(u);
 Result := du;
 v := Add(u, Mult(du, dt/2));
 du := f(v);
 Result := Add(Result, Mult(du, 2));
 v := Add(u, Mult(du, dt/2));
 du := f(v);
 Result := Add(Result, Mult(du, 2));
 v := Add(u, Mult(du, dt));
 du := f(v);
 Result := Add(Result, du);
 Result := Add(u, Mult(Result, dt/6));
end;

procedure TForm1.Button1Click(Sender: TObject);
const
 dt = 0.01; // Zeitschritt
var
 t: Extended;
 u: TData;
begin
 // Anfangswerte
 t := 0;
 u.S := ...
 u.I := ...
 u.R := ...
 while t < 200 do
 begin
  u := RKM(u, dt);
  with Image1.Canvas do
  begin
   Pixels[Round(t), Round(u.S)+100] := clBlack;
   Pixels[Round(t), Round(u.I)+101] := clRed;
   Pixels[Round(t), Round(u.R)+100] := clBlue;
  end;
  t := t + dt;
 end;
end;					

1:
2:
3:
4:
5:
6:

function RKM1(u: TData; dt: Extended): TData;  
var du: TData;
begin  
 du := Mult(f(u), dt);  
 Result := Add(u, du);  
end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:
37:
38:
39:
40:
41:
42:
43:
44:
45:
46:
47:
48:
49:
50:
51:
52:
53:
54:
55:
56:
57:
58:
59:
60:
61:
62:
63:
64:
65:
66:
67:
68:
69:
70:
71:
72:
73:
74:
75:
76:
77:
78:
79:
80:
81:
82:
83:
84:
85:
86:
87:

type
 TData = record
  S, I, R: Extended;
 end;

function f(u: TData): TData;
const
 beta = 0; // Zuwachs
 mu = 0; // Todesrate
 lambda = 0.0004; // = Kontakthäufigkeit*Infektion_Pro_Kontakt
 gamma = 0.1; // = Heilungsrate
begin
 with u do
 begin
  Result.S := beta -mu*S -lambda*S*I;           // susceptible
  Result.I :=      -mu*I +lambda*S*I -gamma*I;  // infected
  Result.R :=      -mu*R             +gamma*I;  // resistant
 end;
end;

function Mult(u: TData; factor: Extended): TData;
begin
 with u do
 begin
  Result.S := S*factor;
  Result.I := I*factor;
  Result.R := R*factor;
 end;
end;

function Add(u, v: TData): TData;
begin
 Result.S := u.S+v.S;
 Result.I := u.I+v.I;
 Result.R := u.R+v.R;
end;

function RKM1(u: TData; dt: Extended): TData;
var du: TData;
begin
 du := Mult(f(u), dt);
 Result := Add(u, du);
end;

function RKM4(u: TData; dt: Extended): TData;
var
 v, du: TData;
begin
 du := f(u);
 Result := du;
 v := Add(u, Mult(du, dt/2));
 du := f(v);
 Result := Add(Result, Mult(du, 2));
 v := Add(u, Mult(du, dt/2));
 du := f(v);
 Result := Add(Result, Mult(du, 2));
 v := Add(u, Mult(du, dt));
 du := f(v);
 Result := Add(Result, du);
 Result := Add(u, Mult(Result, dt/6));
end;

procedure TForm1.Button1Click(Sender: TObject);
const
 ScaleT = 5;
 dt = 0.2;
var
 t: Extended;
 u: TData;
begin
 // Anfangswerte
 t := 0;
 u.S := 1000;
 u.I := 40;
 u.R := 0;
 while t < 70 do
 begin
  u := RKM4(u, dt);
  with Image1,Image1.Canvas do
  begin
   Pixels[Round(ScaleT*t), Height-1-Round(u.S/10)] := clBlack;
   Pixels[Round(ScaleT*t), Height-1-Round(u.I/10)] := clRed;
   Pixels[Round(ScaleT*t), Height-1-Round(u.R/10)] := clBlue;
  end;
  t := t + dt;
 end;
end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:

function Mult(u: TData; factor: Extended): TData;
begin
 with u do
 begin
  Result.S := S*factor;
  Result.I := I*factor;
  Result.R := R*factor;
 end;
end;

function Add(u, v: TData): TData;
begin
 Result.S := u.S+v.S;
 Result.I := u.I+v.I;
 Result.R := u.R+v.R;
end;

function RKM1(u: TData; dt: Extended): TData;
var du: TData;
begin
 du := Mult(f(u), dt);
 Result := Add(u, du);
end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:

procedure TForm2.FormActivate(Sender: TObject);
VAR dt,x :real;
    sp  : byte;
     t: Extended;
     u: TData;
begin
Tab.Cells[0,0]:='Tage';
Tab.Cells[1,0]:='Infizierte';

// Anfangswerte
 x:=1;
 sp:=1;
 dt := 0.2;
 t := 0;
 u.S := 1000;
 u.I := 40;
 u.R := 0;
 while x < 70 do
 begin
  u := RKM1(u, dt);
 Tab.Cells[0,sp]:=FloatToStrF(x,fffixed,10,1);
 [b]Tab.Cells[1,sp]:=FloatToStrF(Round(u.I),fffixed,10,1);[/b]
 t:=t+0.2;
 x:= x +1;
 sp:=sp+1;

  end;
end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:

 while t < 70 do
 begin
  u := RKM1(u, dt);
 Tab.Cells[0,sp]:=FloatToStrF(t*5,fffixed,10,1);
 Tab.Cells[1,sp]:=FloatToStrF(Round(u.I/10),fffixed,10,1);
 t:=t+0.2;

 sp:=sp+1;

  end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:

function Mult(u: TData; factor: Extended): TData;
begin
 with u do
 begin
  Result.S := S*factor;
  Result.I := I*factor;
  Result.R := R*factor;
 end;
end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:

function f(u: TData): TData;
const

a = 0.08; //Nettozuwachsrate Hasen
b = 0.2; //Gewichtsverlust der Füchse/Woche
c = 0.002; // Wahrscheinlichkeit, bei Treffen gefressen zu werden
d = 0.0004; // Beutewahrscheinlichkeit der Füchse
wp = 2000;

begin
 with u do
 begin
Result.H:= a*H*w*0.001 - c*H*F;  {w*0.001}                // Population Hasen
Result.F:= -b*F + d*H*F;                  // " Füchse
Result.w:= Wp-H;                          // Weidekapazität


 end;
end;

function Mult(u: TData; factor: Extended): TData;
begin
 with u do
 begin
  Result.H := H*factor;
  Result.F := F*factor;
 Result.w := w*factor;
 end;
end;

function Add(u, v: TData): TData;
begin
 Result.H := u.H+v.H;
 Result.F := u.F+v.F;
 Result.w := u.w+v.w;
end;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:

const
 N = 10; // 10 Zwischenschritte
 dt = 1/N;
begin
 ...
 while t < 70 do  
 begin  
  for J := 1 to N do
  begin
   t := t + dt;  
   u := RKM4(u, dt);
  end;
  // Hier t und u tabellieren
 end;					

1:
2:
3:
4:

v[0] := u[0] + du[0]*dt/2;
v[1] := u[1] + du[1]*dt/2;
v[2] := u[2] + du[2]*dt/2;
v[3] := u[3] + du[3]*dt/2;