复化Gauss公式求定积分的C语言代码

以下是用C语言实现复化Gauss公式求定积分的代码示例：

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

double f(double x) {
    // 定义被积函数
    return sin(x);
}

double Gauss(double a, double b) {
    // 定义Gauss公式参数
    double x[5] = {-0.9061798459, -0.5384693101, 0, 0.5384693101, 0.9061798459};
    double w[5] = {0.2369268851, 0.4786286705, 0.5688888889, 0.4786286705, 0.2369268851};

    double t = (b - a) / 2;
    double s = (a + b) / 2;
    double sum = 0;

    for (int i = 0; i < 5; i++) {
        double u = t * x[i];
        sum += w[i] * f(s + u);
    }

    return t * sum;
}

double GaussIntegral(double a, double b, double eps) {
    // 复化Gauss公式求定积分
    double h = b - a;
    double t = (a + b) / 2;
    double S = Gauss(a, b);
    double T = Gauss(a, t) + Gauss(t, b);

    if (fabs(T - S) < eps) {
        return T;
    } else {
        return GaussIntegral(a, t, eps / 2) + GaussIntegral(t, b, eps / 2);
    }
}

int main() {
    double a = 0, b = M_PI;
    double eps = 1e-6;
    double result = GaussIntegral(a, b, eps);
    printf("The integral of sin(x) from %.2f to %.2f is %.6f\n", a, b, result);
    return 0;
}

在上述代码中，`f()` 函数定义了被积函数，`Gauss()` 函数定义了Gauss公式的参数和计算方法，`GaussIntegral()` 函数则实现了复化Gauss公式求定积分的递归算法。在 `main()` 函数中，我们指定被积函数为 `sin(x)`，积分区间为 `[0, pi]`，精度为 `1e-6`，并输出计算结果。