Курсовая работа: Решение систем дифференциальных уравнений методом Рунге - Кутты 4 порядка

Difur;

sum:=C[1]+C[2]+C[3];

STROKA;

p:=p+dp; END;

FOR i:=1 TO n DO Cpr[i]:=C[i];

RK_4;

X:=X+dX;

UNTIL(X>Xk);

WRITELN(' - done.');

END;

PROCEDURE INIT;

BEGIN

ClrScr;

WRITELN('Smith-04: v1.0 (c) 1998 by Mike Smith [email protected] ');

WRITELN;

WRITELN;

WRITE('Step 1: Read data from file : in.dat');

ASSIGN(f1,'in.dat');

RESET(f1);

READLN(f1,C[1],C[2],C[3]);

READLN(f1,k1,k2,k3,k4);

READLN(f1,Xn,Xk,dp,n,eps,p);

WRITELN(' - done.');

ASSIGN(f2,'out.rez');

REWRITE(f2);

WRITE('Step 2: Write header to file : out.rez');

WRITELN(f2,'');

WRITELN(f2,'| t,c| Ca,% | Cb,%| Cc,% | SUM | dCa | dCb | dCc |');

WRITELN(f2,'=');