Переписать код с Delphi на Python

var
  Form1: TForm1;
const
   h=0.01;
   imax=50000;
 
implementation
 
{$R *.lfm}
 
{ TForm1 }
 
var
  s,g,ep :real;
  sx,sy,sz, sxdr,sydr,szdr,sxp,syp,szp,enosx,enosy,enosz:real;
  i:Integer;
  sstart, send ,sstep:real;
 
  k1x,k1y,k1z,k2x,k2y,k2z,k3x,k3y,k3z,k4x,k4y,k4z:real;
 
  dvozmx,dvozmy,dvozmz:real;
  x,y,z,vozmx_x,vozmx_y,vozmx_z,vozmy_x,vozmy_y,vozmy_z,vozmz_x,vozmz_y,vozmz_z:real;
  nl1x, nl1y, nl1z, lsum1x, lsum1y, lsum1z: REAL;
  l1x, l1y, l1z, t1, dote1e2, dote2e3, dote1e3 : REAL;
 
 
function fx(x,y,z,t:real):Real;
begin
 result:=2*s*x + y - g*z
end;
function fy(x,y,z,t:real):Real;
begin
 result:= -x
end;
function fz(x,y,z,t:real):Real;
begin
 result:= (x- (8.592 * z - 22.0 * z *z + 14.408 * z *z*z))/ep;
end;
 
function fxc(x,y,z,t:real):Real;
begin
 result:=(2*s*x + y - g*z )/(-x)
end;
function fyc(x,y,z,t:real):Real;
begin
 result:= -x /(-x)
end;
function fzc(x,y,z,t:real):Real;
begin
 result:= ((x- (8.592 * z - 22.0 * z *z + 14.408 * z *z*z))/ep)/(-x);
end;
 
 
procedure rungenorm(var sx,sy,sz:real);
begin
k1x:=h*fx(sx,sy,sz,i*h);
k1y:=h*fy(sx,sy,sz,i*h);
k1z:=h*fz(sx,sy,sz,i*h);
k2x:=h*fx(sx+k1x/2.0,sy+k1y/2.0,sz+k1z/2.0,i*h);
k2y:=h*fy(sx+k1x/2.0,sy+k1y/2.0,sz+k1z/2.0,i*h);
k2z:=h*fz(sx+k1x/2.0,sy+k1y/2.0,sz+k1z/2.0,i*h);
k3x:=h*fx(sx+k2x/2.0,sy+k2y/2.0,sz+k2z/2.0,i*h);
k3y:=h*fy(sx+k2x/2.0,sy+k2y/2.0,sz+k2z/2.0,i*h);
k3z:=h*fz(sx+k2x/2.0,sy+k2y/2.0,sz+k2z/2.0,i*h);
k4x:=h*fx(sx+k3x,sy+k3y,sz+k3z,i*h);
k4y:=h*fy(sx+k3x,sy+k3y,sz+k3z,i*h);
k4z:=h*fz(sx+k3x,sy+k3y,sz+k3z,i*h);
sx:=sx+(k1x+2.0*k2x+2.0*k3x+k4x)/6.0;
sy:=sy+(k1y+2.0*k2y+2.0*k3y+k4y)/6.0;
sz:=sz+(k1z+2.0*k2z+2.0*k3z+k4z)/6.0;
end;
 
 
procedure rungeeno;
begin
  k1x:=-syp*fxc(sx,sy,sz,i*h);
  k1y:=-syp*fyc(sx,sy,sz,i*h);
  k1z:=-syp*fzc(sx,sy,sz,i*h);
  k2x:=-syp*fxc(sx+k1x/2.0,sy+k1y/2.0,sz+k1z/2.0,i*h);
  k2y:=-syp*fyc(sx+k1x/2.0,sy+k1y/2.0,sz+k1z/2.0,i*h);
  k2z:=-syp*fzc(sx+k1x/2.0,sy+k1y/2.0,sz+k1z/2.0,i*h);
  k3x:=-syp*fxc(sx+k2x/2.0,sy+k2y/2.0,sz+k2z/2.0,i*h);
  k3y:=-syp*fyc(sx+k2x/2.0,sy+k2y/2.0,sz+k2z/2.0,i*h);
  k3z:=-syp*fzc(sx+k2x/2.0,sy+k2y/2.0,sz+k2z/2.0,i*h);
  k4x:=-syp*fxc(sx+k3x,sy+k3y,sz+k3z,i*h);
  k4y:=-syp*fyc(sx+k3x,sy+k3y,sz+k3z,i*h);
  k4z:=-syp*fzc(sx+k3x,sy+k3y,sz+k3z,i*h);
  enosx:=sx+(k1x+2.0*k2x+2.0*k3x+k4x)/6.0;
  enosy:=sy+(k1y+2.0*k2y+2.0*k3y+k4y)/6.0;
  enosz:=sz+(k1z+2.0*k2z+2.0*k3z+k4z)/6.0;
end;
 
 
procedure TForm1.Button1Click(Sender: TObject);
begin
s:= floatspinedit1.value;
g:= floatspinedit2.value;
ep:= floatspinedit3.value;
sx:= floatspinedit4.value;
sy:= floatspinedit5.value;
sz:= floatspinedit6.value;
sxdr:= floatspinedit10.value;
sydr:= floatspinedit11.value;
szdr:= floatspinedit12.value;
 
sstart:= floatspinedit7.value; диапазон для дерева
send:= floatspinedit9.value; конечное значение
sstep:= floatspinedit8.value;
xt.clear;
xy.clear;
xtdr.clear;
xycut.clear;
 
for i:=0 to imax do
begin if (abs(sx)>1000) or (abs(sy)>1000) or(abs(sz)>1000)  then break;
   rungenorm(sx,sy,sz);  end;
for i:=imax to imax+10000 do
 
begin
if (abs(sx)>1000) or (abs(sy)>1000) or(abs(sz)>1000)  then break;
if (abs(sxdr)>1000) or (abs(sydr)>1000) or(abs(szdr)>1000)  then break;
sxp:=sx; syp:=sy; szp:=sy;
rungenorm(sx,sy,sz);
rungenorm(sxdr,sydr,szdr);
   if (syp<0) and (sy>0) {and (sxp<0)} then   begin
   rungeeno; xycut.AddXY(enosx,enosy);  end;
xt.AddXY(i*h,sx);
xtdr.AddXY(i*h,sxdr);
xy.AddXY(sx,sy);
end;
 
 
end;
 
procedure TForm1.Button2Click(Sender: TObject); (дерево)
begin
s:= floatspinedit1.value;
g:= floatspinedit2.value;
ep:= floatspinedit3.value;
sx:= floatspinedit4.value;
sy:= floatspinedit5.value;
sz:= floatspinedit6.value;
sstart:= floatspinedit7.value;
send:= floatspinedit9.value;
sstep:= floatspinedit8.value;
tree_s_enosx.clear;
s:=sstart-sstep;
while s<send do begin
s:=s+sstep;
sx:= floatspinedit4.value;
sy:= floatspinedit5.value;
sz:= floatspinedit6.value;
 
for i:=0 to imax do
begin if (abs(sx)>1000) or (abs(sy)>1000) or(abs(sz)>1000)  then break;
    rungenorm(sx,sy,sz);  end;
for i:=imax to imax+10000 do
begin
if (abs(sx)>1000) or (abs(sy)>1000) or(abs(sz)>1000)  then break;
sxp:=sx; syp:=sy; szp:=sy;
rungenorm(sx,sy,sz);
   if (syp<0) and (sy>0) {and (sxp<0)} then   begin
    //rungeeno;  tree_s_enosx.AddXY(s,enosx); end;
    tree_s_enosx.AddXY(s,sx); end;
end;
end;
end;