Runge-Kutta法（龙格－库塔法）的C语言实现

#include <stdio.h>
#include <math.h>

double f(double x, double y)
{
    double result = x/y;
    return result;
}

int main()
{
    double a = 2.0;
    double b = 2.6;
    double h = 0.2;
    double n0 = (b-a)/h + 1;
    int n = n0;
    printf("n = %d\n", n);
    
    double x[n], y[n];
    y[0] = 1;
    
    for (int i = 0; i < n; i++)//initialize x[n]
    {       
        x[i] = a + h * i;
    }
    
    double k1,  k2, k3, k4;
    for (int i = 0; i < n-1; i++)//calculate y[n]
    {
        k1 = f(x[i], y[i]);
        k2 = f(x[i]+h/2, y[i]+h*k1/2);
        k3 = f(x[i]+h/2, y[i]+h*k2/2);
        k4 = f(x[i]+h, y[i]+h*k3);
        y[i+1] = y[i] + h/6*(k1 + 2*k2 + 2*k3 + k4);
    }

    for (int i = 0; i < n; i++)//print
    {       
        printf("y(%lf) = %lf\n", x[i], y[i]);
    }
    
    return 0;
}