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:

procedure Tmain.punktberechnen;
var temp:string;
fEl_nachoben: boolean;
vV1,vH1: extended; //temporär
//Weitere Variablen:
{
v0: Gesamtgeschwindigkeit
vH0,vV0: vertikale, horizontale Geschwindigkeit
B,E: Feldstärken (B=magn.; E=el.)
Q,M, usw sollte klar sein..
k: Kräftemultiplier (damit man was sieht ;-)
}
begin
 //MAGNETFELD
  v0 := sqrt(Power(vH0,2)+Power(vV0,2));  
  fL := B * Q * v0;
  fL := fL * k;
   //Alten "Vektor" berechnen
    vH1 := vH0; //vH1: temporäre horizointale Geschwindigkeit
    vV1 := vV0; //vV1: temporäre vertikale Geschwindigkeit
   //"Vektor" um 90° drehen
    vH1 := vV0;
    vV1 := -vH0;
   //Stärke der temporären, "gedrehten" Geschwindigkeit
    vH1 := vH1 * fL;
    vV1 := vV1 * fL;
   //Zur Globalen Geschwindigkeit hinzufügen...
    vH0 := vH0 + vH1;
    vV0 := vV0 + vV1;

 //ELEKTRISCHES FELD
  fEL := Q * E;
  fEL := fEL * k;
  vV0 := vV0 + fEl //Zur Globalen Geschwindigkeit hinzufügen...

 //GRAVITATIONSFELD
  ... //An der liegts nicht, schon getestet..

//POSITION
 eY := eY - vV0; //Da die Geschwindigkeit bei jedem Durchlauf
 eX := eX + vH0; //addiert wird, ersetzt sie die Beschleunigung

//KOLLISION
 ...

//TEXT OUTPUT
 ...

//ZEICHNEN
 ...                    

1:
2:

    vH1 := vH0; //was soll das sein
    vV1 := vV0; //was soll das sein                    

1:
2:

    vH1 := vV0;
    vV1 := -vH0;                    

1:
2:

    vH1 := vH1 * fL;//versteh nicht was das soll
    vV1 := vV1 * fL;//soll das die Zentrifugalkraft werden                    

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:

type TVec=record
x:extended;
y:extended;
z:extended;
end;

implementation

function skalarprodukt(vec1,vec2:TVec):extended;
begin
   result:=vec1.x*vec2.x+vec1.y*vec2.y+vec1.z*vec2.z;
end;

function kreuzprodukt(vec1,vec2:TVec):TVec;
begin
   result.x:=vec1.y*vec2.z-vec1.z*vec2.y;
   result.y:=-vec1.x*vec2.z+vec1.z*vec2.x;
   result.z:=vec1.x*vec2.y-vec1.y*vec2.x;
end;

//F=Q*(E + v x B) Q=Ladung E=el.Feld B=magn Feld v=Geschwindigkeit
//result=Kraftvektor
function Lorentzforce(Q:integer;E,B,v:TVec):TVec;
var BForce:TVec;
begin
   BForce:=kreuzprodukt(v,B);
   
   result.x:=Q*(E.x+BForce.x);
   result.y:=Q*(E.y+BForce.y);
   result.z:=Q*(E.z+BForce.z);
end;

//dp/dt=F Newton =>m*deltav=F*deltatime
//Pneu=Palt+v*deltatime
function newpos(istpos,istv,Force:TVec):TVec;
begin
   //das musste aber noch selber ausfüllen
   //hier könntest du zum beispiel variable schrittweiten zur genaueren Analyse etc. 
   //einbauen
end;                    

1:
2:

Fx := Q*(vX*B);
Fy := Q*(E - vY*B);