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

r1[i]:=dC[i];

C[i]:=cPR[i]+r1[i]*(dX/2);

END;

Difur;

FOR i:=1 TO n DO BEGIN

r2[i]:=dC[i];

C[i]:=cPr[i]+r2[i]*(dX/2);

END;

Difur;

FOR i:=1 TO n DO BEGIN

r3[i]:=dC[i];

C[i]:=cPR[I]+r3[i]*dX;

END;

Difur;

FOR i:=1 TO n DO r4[i]:=dC[i];

FOR i:=1 TO n DO rSR[i]:=((r1[i]+r2[i])*(r2[i]+r3[i])*(r3[i]+r4[i]))/6;

END;

PROCEDURE STROKA;

BEGIN

WRITE(f2,'|',x:4:1,'|',c[1]:7:3,'|',c[2]:7:3,'|',c[3]:7:3,'|');

WRITE(f2,sum:3:0,'|',dc[1]:7:3,'|',dc[2]:7:3,'|',dc[3]:7:3,'|');

WRITELN(f2);

END;

PROCEDURE RUN;

BEGIN

WRITE('Step 3: Calculating data and writting results to file : out.rez');

X:=Xn;

dX:=0.05;

REPEAT