Рунге-Кутта второго порядка c++

[Sylar220]

Нужно написать формулу Рунге Кутта 2го порядка вместо 4го, которая в коде.


for (int i = 1; i <= n; i++)
{
k1 = Y(i)+h/2*runge(X[i - 1], Y[i - 1]);
cout << "k1 = " << k1 << " " << endl;
k2 = h*runge(X[i - 1] + h / 2, Y[i - 1] + k1 / 2);
cout << "k2 = " << k2 << " " << endl;
k3 = h*runge(X[i - 1] + h / 2, Y[i - 1] + k2 / 2);
cout << "k3 = " << k3 << " " << endl;
k4 = h*runge(X[i - 1] + h, Y[i - 1] + k3);
cout << "k4 = " << k4 << " " << endl;
dy = (k1 + 2 * k2 + 2 * k3 + k4) * (1.00 / 6.00);
cout << "dy = " << dy << " " << endl;
Y[i] = Y[i - 1] + dy;
cout << "Y[" << i << "]= " << Y[i] << " " << endl;
}