using System;

​

public class Program {

    static void Main(string[] args) {

        double x0 = 0;  // Valor inicial de x

        double y0 = 1;  // Valor inicial de y

        double h = 0.1; // Tamaño del paso

        double xf = 1;  // Valor final de x

        int n = (int)((xf - x0) / h); // Número de pasos

        double[] x = new double[n+1]; // Vector de valores de x

        double[] y = new double[n+1]; // Vector de valores de y

​

        // Valores iniciales

        x[0] = x0;

        y[0] = y0;

​

        // Método de Runge-Kutta de orden 4

        for (int i = 0; i < n; i++) {

            double k1 = h * (5 * x[i] - 2 * y[i]);

            double k2 = h * (5 * (x[i] + h / 2) - 2 * (y[i] + k1 / 2));

            double k3 = h * (5 * (x[i] + h / 2) - 2 * (y[i] + k2 / 2));

            double k4 = h * (5 * (x[i] + h) - 2 * (y[i] + k3));

            y[i+1] = y[i] + (k1 + 2 * k2 + 2 * k3 + k4) / 6;

            x[i+1] = x[i] + h;

        }

​

        // Imprimir los resultados

        Console.WriteLine(" x      y(x)");

        Console.WriteLine("------------");

        for (int i = 0; i <= n; i++) {

            Console.WriteLine("{0:F2}   {1:F4}", x[i], y[i]);

        }

    }

}

​