Нарисовать дуанты циклотрона и начертить траекторию движения частицы

Uses Graph;
Label 1;
Var q,m,v,vx,vy,vz,x,y,z,r,rd,a,bz,ex,e0,w,t,dt,mx,fx,fy:Real;
gd,gm,xe,ye:Integer;
Begin
Rd:=1; a:=5e-2; bz:=1; q:=2*1,6e-19; m:=4*1,7e-27;
Mx:=200;
W:=q*bz/m; x:=0; y:=0; yx:=0; vy:=0; dt:=2e-10;
Gd:=0;  InitGraph(gd,gm,’C:\tp\bgi’);
SetColor(3); Circle(315,225,Round(rd*mx));
Line(315-Round(a*mx), 225-Round(rd*mx),
315-Round(a*mx), 225+Round(rd*mx));
Line(315+Round(a*mx),
225-Round(rd*mx), 315+Round(a*mx),
225+Round(rd*mx);
Repeat
If (x<-a/2) or (x>a/2) then e0:=0 else e0:=3e6;
Ex;+e0*cos(w*t); fx:=q*(ex+vy*bz); fy:=q*vx*bz;
Vx:=vx+fx*dt/m; vy:=vy+fy*dt/m;
X:=x+vx*dt; y:=y+vy*dt; t:=t+dt;
Xe:=314+Round(x*mx); ye:=225-Round(y*mx);
If (xe<0) or (xe>630) or (ye<0) or (ye>450) then
Goto 1
PutPixel(xe,ye,15);
R:=sqrt(x*x+y*y); If (r>=rd) and (Abs(x)<a/2) then
Bz:=0
Until t>1e-5;
1:writeln(‘v=’,sqrt(vx*vx+vy*vy),’m/s’);
Readln; CloseGraph;
End.