Курсовая работа: Моделирование движения на плоскости

Md[i]:=M0+ln(fi[i]+1)+sqrt(fi[i]);

int[i]:=(((Md[i]-Mc)+(Md[i-1]-Mc))*dfp)/2;

W[i]:=sqrt((2/I0)*(I0*sqr(W[i-1])/2+int[i]));

Wcp:=(W[i]+W[i-1])/2;

t[i]:=t[i-1]+(fi[i]-fi[i-1])/Wcp;

b[i]:=(W[i]-W[i-1])/(t[i]-t[i-1]) end;

Tp:=t[n+1];

write(fe,' ':10);

for i:=1 to 57 do

write(fe,'_');writeln(fe);

writeln(fe,' ':10,'I',' ':10,'I',' ':10,'I',' ':16,'I',' ':16,'I');

writeln(fe,' ':8,' I fi[i] I W[i] I b[i] I ',

' t[i] I ');

writeln(fe,' ':10,'I',' ':10,'I',' ':10,'I',' ':16,'I',' ':16,'I');

write(fe,' ':10);

for i:=1 to 57 do write(fe,'-');writeln(fe);

for i:=1 to n+1 do

writeln(fe,' ':7,i:2,' I',' ',fi[i]:7:3,' I',' ',W[i]:7:3,' I',' ',

b[i]:7:3,' I',' ',t[i]:7:3,' I');

writeln('Произведен расчет параметров разгона');

fit:=(I0*sqr(W[n+1]))/(2*Mc);

b[n+1]:=-Mc/I0;

dft:=fit/n;

for i:=n+2 to (2*n+1) do begin

fi[i]:=fi[i-1]+dft;

W[i]:=sqrt((2/I0)*((I0*sqr(W[i-1])/2)-(Mc*(fi[i]-fi[i-1]))));

Wcp:=(W[i]+W[i-1])/2;

t[i]:=t[i-1]+(fi[i]-fi[i-1])/Wcp;

b[i]:=(W[i]-W[i-1])/(t[i]-t[i-1]) end;