Курсовая работа: Численные методы при решении задач

f (&y0[i * n], k1, xi); // Вычислим f(xi, yi) = k1 / h

// И для каждого дифференциального уравнения системы проделываем

// операции вычисления k1, а также подготовки в ya аргумента для

// вычисления k2

for (j = 0; j < n; j++)

{

k1[j] *= h; // Вычислим наконец-то k1

ya[j] = y0[i*n+j] + k1[j] / 2.;

// И один из аргументов для функции

} // вычисления k2

f (ya, k2, xi + (h / 2.)); // Вычислим f(xi,yi) = k2 / h

for (j = 0; j < n; j++)

{ // Вычислим наконец-то k2

k2[j] *= h;

ya[j] = y0[i*n+j] + k2[j] / 2.; // И один из аргументов для функции

} // вычисления k3

f (ya, k3, xi + h / 2.); // Вычислим f(xi,yi) = k3 / h

for (j = 0; j < n; j++)

{

k3[j] *= h; // Вычислим наконец-то k3

ya[j] = y0[i*n+j] + k3[j]; // И один из аргументов для функции

} // вычисления k4

f (ya, k4, xi + h); // Вычислим f(xi,yi) = k4 / h

for (j = 0; j < n; j++) k4[j] *= h; // Вычислим наконец-то k4

// Надо вычислить приращение каждой функции из n

for (j = 0; j < n; j++) // Вычисляем следующее значение

// функции

// Y[i+1] = Yi + ...

y0[(i+1)*n+j] = y0[i*n+j] + (k1[j] + 2. * k2[j] + 2 * k3[j] + k4[j]) / 6.;