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:

unit Unit1;

interface

uses
  Windows, math, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs,
  ExtCtrls, StdCtrls;

type
  TForm1 = class(TForm)
    Timer1: TTimer;
    Button1: TButton;
    procedure Timer1Timer(Sender: TObject);
    procedure Button1Click(Sender: TObject);
  private
    { Private-Deklarationen }
  public
    { Public-Deklarationen }
  end;

var
  Form1: TForm1;
  geschw,speed,x,y:array[1..2] of single;
  alpha,xdelta,ydelta,kraft,s:single;
implementation

{$R *.DFM}
procedure sammeln;
begin
  xdelta:=x[1]-x[2];  //die zwei seitenlängen
  ydelta:=y[1]-y[2];  //zur berechnung von s
  s:=sqrt(sqr(xdelta)+sqr(ydelta));  //abstand punkt->gravitation
  if ydelta<>0 then
    alpha:=arcsin(xdelta/s) else
    alpha:=arccos(xdelta/s);  //<--ignorieren^^'
  speed[1]:=speed[1]-sin(alpha);  //geschwindigkeit und horizontaler richtung
  if ydelta>0 then       //nötig damit die kräfte nicht ins negativ umschlagen
  speed[2]:=speed[2]-cos(alpha) //vertikale geschwindigkeit
  else
  speed[2]:=speed[2]+cos(alpha); //vertikale geschwindigkeit
  kraft:=1/(s*s);  //wirkende kraft
  geschw[1]:=speed[1]*kraft; //keine funktion
  geschw[2]:=speed[2]*kraft; //keine funktion
  x[2]:=x[2]-speed[1]*kraft; //letztendliche lageänderung
  y[2]:=y[2]-speed[2]*kraft; //letztendliche lageänderung
end;

procedure zeichnen;
var i:byte;
begin
  for i:=1 to 2 do
  Form1.Canvas.Ellipse(trunc(x[i]*10+100)-5,trunc(y[i]*10+100)-5,trunc(x[i]*10+100)+5,trunc(y[i]*10+100)+5);
end;

procedure TForm1.Timer1Timer(Sender: TObject);
begin
sammeln;
zeichnen;
end;

procedure TForm1.Button1Click(Sender: TObject);
begin
  speed[1]:=0;  //horizontale geschwindigkeit
  speed[2]:=50; //vertikale anfangsgeschwindigkeit
  x[1]:=17.5;
  x[2]:=32.3;
  y[1]:=1.5;
  y[2]:=15.3;
  timer1.Enabled:=true;
end;

end.					

1:
2:
3:
4:
5:

  speed[1]:=speed[1]-sin(alpha)*kraft;
  if ydelta>0 then
  speed[2]:=speed[2]-cos(alpha)*kraft
  else
  speed[2]:=speed[2]+cos(alpha)*kraft;					

1:
2:
3:

  x[2]:=x[2]-speed[1]*kraft; //letztendliche lageänderung   
  y[2]:=y[2]-speed[2]*kraft; //letztendliche lageänderung   
end;					

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

  kraft:=1/(s*s);
  speed[1]:=speed[1]-sin(alpha)*kraft;  //hier kommt die änderung hin
  if ydelta>0 then
  speed[2]:=speed[2]-cos(alpha)*kraft
  else
  speed[2]:=speed[2]+cos(alpha)*kraft;
  x[2]:=x[2]-speed[1];  //und hier fehlt
  y[2]:=y[2]-speed[2];  //auch das "*kraft"					

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:

type
TMassepunkt = packed record
  x: Extended; //X in Meter
  y: Extended; //Y in Meter
  m: Extended; //Masse in kg
end;

var
  pts: array [1..2] of TMassepunkt;
  r: Extended;  //Abstand
  F: Extended;  //Gravitationskraft
  dx, dy: Extended;

begin
  dx := pts[2].x - pts[1].x;
  dy := pts[2].y - pts[1].y;

  r := sqrt(dx * dx + dy * dy);
  
  F := Grav * pts[1].m * pts[2].m / r;  //Für r = 0 prüfen

  pts[1].x := pts[1].x + F / 2 / pts[1].m * (dx / r); //For Time Scaling: * (t * t);
  pts[1].y := pts[1].y + F / 2 / pts[1].m * (dy / r); //For Time Scaling: * (t * t);

  pts[2].x := pts[2].x + F / 2 / pts[2].m * (dx / r); //For Time Scaling: * (t * t);
  pts[2].y := pts[2].y + F / 2 / pts[2].m * (dy / r); //For Time Scaling: * (t * t);
end;					

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

  kraft:=1/(s*s);   
  speed[1]:=speed[1]-sin(alpha)*kraft;  //hier kommt die änderung hin   
  if ydelta>0 then   
  speed[2]:=speed[2]-cos(alpha)*kraft   
  else   
  speed[2]:=speed[2]+cos(alpha)*kraft;   
  x[2]:=x[2]-speed[1];  //und hier fehlt   
  y[2]:=y[2]-speed[2];  //auch das "*kraft"					

1:

speed[2]:=speed[2]+cos(alpha)*kraft;					

1:

x[2]:=x[2]-speed[1];					

1:

F := Grav * pts[1].m * pts[2].m / r;  //Für r = 0 prüfen