42 std::vector<double>& rCurrentYValues,
43 std::vector<double>& rNextYValues)
53 std::vector<double> k1(num_equations);
54 std::vector<double>& dy = rNextYValues;
59 for (
unsigned i=0; i<num_equations; i++)
61 k1[i] = timeStep*dy[i];
62 k1[i] = k1[i]/2.0 + rCurrentYValues[i];
67 for (
unsigned i=0; i<num_equations; i++)
69 rNextYValues[i] = rCurrentYValues[i] + timeStep*dy[i];