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

C,Mc,I0,Wcp,fip,dfp,fit,dft,Tp,Tt,M0:real;

begin clrscr;

assign(fe,'kurs-13v.rez');rewrite(fe);writeln(fe);

writeln(fe,' ':15,'Определение параметров вращательного движения',

' тела');

writeln(fe);

writeln(fe, ' ':40,' Лабоцкий Д.В.');

writeln(fe);

writeln(fe, ' ':30,'Вариант 13');

writeln('Введите исходные данные');

write('Момент инерции тела равен I0= ');readln(I0);

write('Коэффициент для движущего момента равен М0= ');readln(M0);

write('Момент сопротивления равен Мc= ');readln(Mc);

write('Уголразгона fip= ');readln(fip);

write('Количество интервалов разбиения n= ');readln(n);

writeln('Исходные данные занесены в файл результатов');

writeln(fe);

writeln(fe,' ':25,'Исходныеданные');

writeln(fe);

writeln(fe,' ':10,'Момент инерции тела равен I0= ',I0:5:2,' кг/м2');

writeln(fe,' ':10,'Коэффициент движущего момента М0= ',M0:5:2,' нм');

writeln(fe,' ':10,'Момент сопротивления Мc= ',Mc:5:2,' нм');

writeln(fe,' ':10,'Угол разгона fip= ',fip:5:2,' рад');

writeln(fe,' ':10,'Количество интервалов разбиения n= ',n:2);

dfp:=fip/n;

fi[1]:=0;W[1]:=0;t[1]:=0;

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

b[1]:=(Md[1]-Mc)/I0;

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