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.