Как сделать модель движения заряженной частицы в однородном магнитном поле? И как сделать эту модель визуально видимой?

uses graphABC;
procedure pole;
var d,i:integer;
begin
d:=windowheight div 5;
setpencolor(clBlue);
setpenwidth(2);
for i:=1 to 4 do
 begin
  line(0,d*i,windowwidth,d*i);
  line(0,d*i,10,d*i-5);
  line(0,d*i,10,d*i+5);
 end;
end;
procedure zaryad(x,y:integer);
begin
setpenwidth(1);
setpencolor(clRed);
setbrushcolor(clRed);
circle(x,y,10);
setpencolor(clWhite);
setbrushcolor(clWhite);
rectangle(x-8,y-3,x+8,y+3);
rectangle(x-3,y-8,x+3,y+8);
end;
var r,h,yc,x,y:integer;
    xc,u:real;
begin
setwindowsize(600,400);
centerwindow;
r:=windowheight div 2-50;
h:=r div 2;
pole;
yc:=windowheight div 2;
xc:=h;
u:=-pi;
zaryad(h,yc);
lockdrawing;
 
repeat
clearwindow;
pole;
x:=round(xc+h*cos(u));
y:=round(yc+r*sin(u));
zaryad(x,y);
sleep(5);
xc:=xc+0.1;
u:=u+0.01;
redraw;
until x>windowwidth;
end.