C#龙格库塔解微分方程

Question

请问大家，C#解龙格库塔该怎么写，刚学一点C#基本语句，只会定义k的四行，不知道完整的该怎么写

Answer 1

float dydx(float x, float y)
{
    return ((x - y) / 2);
}
float rungeKutta(float x0, float y0, float x, float h)
{
    // Count number of iterations using step size or
    // step height h
    int n = (int)((x - x0) / h);

    float k1, k2, k3, k4, k5;

    // Iterate for number of iterations
    float y = y0;
    for (int i = 1; i <= n; i++)
    {
        // Apply Runge Kutta Formulas to find
        // next value of y
        k1 = h * dydx(x0, y);
        k2 = h * dydx(x0 + 0.5f * h, y + 0.5f * k1);
        k3 = h * dydx(x0 + 0.5f * h, y + 0.5f * k2);
        k4 = h * dydx(x0 + h, y + k3);

        // Update next value of y
        y = y + (1.0f / 6.0f) * (k1 + 2 * k2 + 2 * k3 + k4); ;

        // Update next value of x
        x0 = x0 + h;
    }

    return y;
}

//调用：
rungeKutta(1, 2, 3, 4);

/*函数引用自：
https://www.geeksforgeeks.org/runge-kutta-4th-order-method-solve-differential-equation/
*/

题主可以试一下。