四阶显式Adams法求方程组C语言,第五讲第4章线性多步法续.doc

以下是使用四阶显式Adams法求解一阶常微分方程组的C语言代码：

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

#define N 2 // 定义方程组维数
#define h 0.1 // 步长
#define T 1.0 // 终止时间

double f[N](double t, double y[]); // 定义函数f

int main() {
    double t, y[N], k1[N], k2[N], k3[N], k4[N];
    int i;

    // 初始化y和t
    y[0] = 0.0;
    y[1] = 1.0;
    t = 0.0;

    // 输出初始值
    printf("%f %f %f\n", t, y[0], y[1]);

    // 迭代求解
    while (t < T) {
        // 计算k1
        for (i = 0; i < N; i++) {
            k1[i] = h * f[i](t, y);
        }
        // 计算k2
        for (i = 0; i < N; i++) {
            k2[i] = h * f[i](t + h / 2.0, y + k1 / 2.0);
        }
        // 计算k3
        for (i = 0; i < N; i++) {
            k3[i] = h * f[i](t + h / 2.0, y + k2 / 2.0);
        }
        // 计算k4
        for (i = 0; i < N; i++) {
            k4[i] = h * f[i](t + h, y + k3);
        }
        // 计算下一个节点的y值
        for (i = 0; i < N; i++) {
            y[i] += (k1[i] + 2.0 * k2[i] + 2.0 * k3[i] + k4[i]) / 6.0;
        }
        // 更新t
        t += h;
        // 输出结果
        printf("%f %f %f\n", t, y[0], y[1]);
    }

    return 0;
}

// 定义函数f
double f[N](double t, double y[]) {
    double f[N];
    f[0] = y[1];
    f[1] = -y[0];
    return f;
}

注意：代码中的函数f是一个数组，可以根据实际情况调整。