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

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);