复化辛普森求积分c语言简单程序

复化辛普森法（Composite Simpson's Rule）是一种数值积分方法，用于近似计算函数在某一区间上的定积分。在C语言中，实现这个算法需要一些数学计算和循环结构。以下是一个简单的C语言程序示例，用于计算给定函数在一个给定区间的定积分为三次近似：

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

// 复化辛普森公式函数
double simpson(double (*func)(double), double a, double b, int n) {
    double h = (b - a) / n;
    double sum = func(a) + func(b); // 边界项
    for (int i = 1; i < n; i += 2) {
        sum += 4 * func(a + i * h); // 奇数位置的函数值
    }
    for (int i = 2; i < n - 1; i += 2) { // 偶数位置的函数值
        sum += 2 * func(a + i * h);
    }
    return h * sum / 3;
}

// 示例函数，你需要根据实际需求替换为你的函数
double example_func(double x) {
    // 这里只是一个例子，替换为你的函数
    return sin(x);
}

int main() {
    double a = 0.0;
    double b = M_PI; // 例如，考虑0到π的正弦函数积分
    int n = 10; // 分段数量，越多精度越高

    double integral_approx = simpson(example_func, a, b, n);
    printf("复化辛普森法估算的积分值: %lf\n", integral_approx);

    return 0;
}