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

writeln(' h=',(h):8:5,'; e2max=',e2max:16:8,'; e4max=',e4max:16:8,' ',n:8);

writeln(outfile);

writeln(outfile,' h=',h:6:5,'; e2max=',e2max:10:8,'; e4max=',e4max:10:8,'.');

close(outfile);

readkey;

end.

Функция вычисления точного решения

function clearsolve(t:real):real;

begin

clearsolve:=exp(-t)*(cos(t)+sin(t)+t*sin(t));

end;

Процедура вычисления правых частей системы уравнений в нормальной форме Коши

procedure right(t:real;var x,f:vector_n);

begin

f[1]:=x[2];

f[2]:=2*exp(-t)*cos(t)-2*x[1]-2*x[2];

end;

Процедура RK2

procedure RK2(t:real;h:real;var x:vector_n);

var h4, h23: real;

f1, f2, xr: vector_n;

begin

h4:=0.25*h;

h23:=0.66666667*h;

right(t,x,f1);

for j:=1 to 2 do xr[j]:=x[j]+h23*f1[j];

right(t+h23,xr,f2);

for j:=1 to 2 do x[j]:=x[j]+h4*(f1[j]+3*f2[j]);

end;

Процедура RK4

procedure RK4(t:real;h:real;var x:vector_n);

var h2,h6:real;

f,fs,xr:vector_n;