Рунге-Кутт әдісінің программа листингі

uses crt;

label 1;

var x0,y0,x,y,h,a,b,c:real;

n,i:integer;

function f(x,y:real):real;

begin

f:=x+sin(y/2.8);

begin

textbackground(1);

clrscr;

writeln('****************************************************');

writeln('* Данной программой определяетсяся *');

writeln('* численное решение дифференциального уравнения *');

writeln('* методом Рунге-Кутта *');

writeln('* *');

writeln('****************************************************');

writeln;

window(15,5,65,15);

textbackground(0);

clrscr;

writeln('Введите исходные данные x0, y0, h, n');

readln(x0,y0,h,n);

x:=x0;

i:=0;

1: a:=h*f(x,y);

x:=x+h/2;

y:=y0+a/2;

b:=h*f(x,y);

y:=y0+b/2;

c:=h*f(x,y);

x:=x+h/2;

y:=y0+c;

y:=y0+(a+2*(b+c)+h*f(x,y))/6;

writeln(x:7:4,' ',y:7:4);

if i>n then exit else

i:=i+1;

y0:=y;

end.
*************** ****************** *****************