Рунге-Кутт!!!

[Ёжик86]

Составить программу для решения задачи Коши для дифференциального уравнения первого порядка по методу Рунге-Кутта
y'=-(y*y)/(x*x); y(2)=2; x на промежутке [2;3]
это все что дано в условии

[Fog]

Метод Рунге-Кутты решения дифференциального уравнения первого порядка y' = F(x,y) является методом четвертого порядка, в котором вычисления производятся по формуле:
yk+1 = yk +(k1 +2k2 +2k3 +k4 )/6, где
k1 = Fk h = F(xk , yk )h
k2 = F(xk +h/2, yk +k1 /2)h
k3 = F(xk +h/2, yk +k2 /2)h
k4 = F(xk +h, yk +k3 )h,
k = 0, ..., n-1
h = (xf -x0 )/n
Приведенный ниже код решает дифференциальное уравнение на отрезке от x до x1. Начальная точка имеет кординаты (x,y). Расчет происходит с постоянным шагом h=(x1-x)/n, где n - количество точек, заданное пользователем.
...

i : Integer;
  h, y1, k1, k2, k3: Double;
begin
  h := (x1-x)/n;
  y1 := y;
  i := 1;
  repeat
    k1 := h*F(x, y);
    x := x+h/2;
    y := y1+k1/2;
    k2 := F(x, y)*h;
    y := y1+k2/2;
    k3 := F(x, y)*h;
    x := x+h/2;
    y := y1+k3;
    y := y1+(k1+2*k2+2*k3+F(x, y)*h)/6;
    //Вывод решения в таблицу
    ListBox1.Items.Add(FormatFloat('y(0.###) = ', x) +
       FormatFloat('0.###', y));
    //Вывод решения на график
    Series1.AddXY(x, y);
    y1 := y;
    i := i+1;
  until  not (i<=n);
end;

[Fog]

Прогу качай тут: http://depositfiles.com/ru/files/8551856
Еще посмотри тут http://plagiata.net.ru - пробей по поиску метод Рунге Кутт

[Ёжик86]

мне на паскале надо(((
или я вообще ничего не понимаю(((