用c语言写一个高斯求积公式

以下是用C语言实现的高斯求积公式代码：

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

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

double gauss(double a, double b) {
    double x1, x2, x3, w1, w2, w3, sum;
    x1 = -0.77459667;
    x2 = 0;
    x3 = 0.77459667;
    w1 = w2 = w3 = 5.0 / 9.0;
    sum = w1 * f((b - a) * x1 / 2.0 + (a + b) / 2.0)
        + w2 * f((b - a) * x2 / 2.0 + (a + b) / 2.0)
        + w3 * f((b - a) * x3 / 2.0 + (a + b) / 2.0);
    sum *= (b - a) / 2.0;
    return sum; // 返回积分结果
}

int main() {
    double a, b, result;
    printf("请输入积分区间[a, b]: ");
    scanf("%lf %lf", &a, &b);
    result = gauss(a, b);
    printf("积分结果为: %lf", result);
    return 0;
}

其中，`f(x)` 函数定义了被积函数，`gauss(a, b)` 函数通过高斯求积公式计算积分结果，`main()` 函数接受用户输入的积分区间，调用 `gauss(a, b)` 函数计算积分结果，并输出到屏幕上。