mayatnk_L2

program mayatnik_L2;
  label M;
  
   var
      i:integer;
      l1,l2,om,g,ddx1,dx1,x1,ddfi,dfi,fi,dt,t:real;
      ddx2,dx2,x2,m1,m2,b,pi,la1,ta1:real;
      dds,ds,s,x0,c2,la2,ta2,E,J1,ar,mu,h,lot:real;
  BEGIN
      E:=20000000000.0; J1:=0.00016;
      m1:=1000.0; l1:=4.0; om:=2.0; pi:=3.14;
       g:=9.81; b:=1.0; h:=0.0;
       ddx1:=0.0; dx1:=0.0; x1:=0.0;
       ddx2:=0.0; dx2:=0.0; x2:=0.0;
       ddfi:=0.0; dfi:=0.0; fi:=0.0;
       dds:=0.0; ds:=0.0; s:=0.0;
        t:=0.0; dt:=0.001; i:=0; 
        lot:=0.5; mu:=0.5;
        l2:=l1*lot;
        m2:=m1*mu;
       mu:=1.0; 
        la1:=sqrt(3.0*E*J1/m1/exp(3*ln(l1))); ta1:=2.0*pi/la1;
        la2:=sqrt(g/(lot*l1)); 
        ar:=0.0;   
       M:
       i:=i+1; t:=t+dt;
        x0:=t*exp(-b*t)*Sin(2.0*pi*om*t);
    ddx1:=sqr(la1)*(x0-x1)/(1.0+mu)-mu/(1+mu)*(g*Cos(fi)+sqr(dfi)*l2)*Sin(fi);
    ddfi:=-(g*Sin(fi)+ddx1*Cos(fi))/l2;
    dx1:=dx1+ddx1*dt;
    dfi:=dfi+ddfi*dt;
    x1:=x1+dx1*dt;
    fi:=fi+dfi*dt;
    x2:=l2*Sin(fi);
  
       if abs(x0-x1)>ar then ar:=abs(x0-x1);
   if i div 100*100=i then begin
   { writeln(t:10:3,x0:10:5,ddx1:10:5,ddx2:10:5,ddfi:10:5);
   writeln(t:10:3,x0:10:5,dx1:10:5,dx2:10:5,dfi:10:5); }
    writeln(t:10:3,x0:20:10,x1:20:10,x2:20:10,fi:20:10);
    writeln(ar:10:5); end;
    if t<4.0 then goto M;
    {writeln(x1:10:5);}
  END.