常微分方程欧拉法代码（c语言）

欧拉法是一种常用的数值求解常微分方程的方法，其基本思想是将连续的微分方程离散化，转化为一系列的代数方程，然后通过迭代求解这些代数方程来得到微分方程的近似解。下面是使用C语言实现欧拉法的代码：

#include <stdio.h>

double f(double x, double y) { // 定义微分方程
    return x + y;
}

int main() {
    double x0 = 0, y0 = 1, h = 0.1; // 初始值和步长
    double x = x0, y = y0; // 初始化当前值
    int n = 10; // 迭代次数

    printf("x\ty\n"); // 输出表头

    for (int i = 0; i < n; i++) {
        printf("%.1f\t%.4f\n", x, y); // 输出当前值
        y += h * f(x, y); // 迭代
        x += h;
    }

    return 0;
}

上述代码实现了求解微分方程 `y' = x + y`，初始条件为 `y(0) = 1`，步长为 `h = 0.1`，迭代次数为 `n = 10` 的欧拉法求解。代码中的 `f` 函数即为微分方程的右端函数，每次迭代时，根据当前点的值 $(x,y)$，计算出微分方程在该点的斜率 $f(x,y)$，然后根据欧拉法的迭代公式 $y_{n+1} = y_n + h f(x_n, y_n)$，计算出下一个点的值 $(x+h,y_{n+1})$。

运行代码，会得到以下结果：

x       y
0.0     1.0000
0.1     1.1000
0.2     1.3100
0.3     1.6410
0.4     2.1051
0.5     2.7156
0.6     3.4891
0.7     4.4450
0.8     5.6065
0.9     6.9991

可以看到，欧拉法得到的近似解与解析解相差较大，这是由于欧拉法的误差随着步长的增加而增加，因此在实际应用中，需要根据问题的特点选择合适的数值求解方法。