复化Simpson公式求定积分及其误差和调用函数次数的C语言代码

下面是复化Simpson公式求定积分及其误差和调用函数次数的C语言代码：

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

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

double Simpson(double a, double b, double (*f)(double), int n) { // 复化Simpson公式
    double h = (b - a) / n;
    double sum = 0;
    for (int i = 1; i < n; i += 2) {
        sum += f(a + i * h) * 4;
    }
    for (int i = 2; i < n - 1; i += 2) {
        sum += f(a + i * h) * 2;
    }
    return (f(a) + f(b) + sum) * h / 3;
}

int main() {
    double a, b, eps;
    int n = 2;
    printf("请输入积分区间[a, b]: ");
    scanf("%lf%lf", &a, &b);
    printf("请输入误差限eps: ");
    scanf("%lf", &eps);
    double I1, I2 = 0;
    do {
        I1 = I2;
        I2 = Simpson(a, b, f, n);
        n *= 2;
    } while (fabs(I2 - I1) >= eps);
    printf("积分值: %lf\n", I2);
    printf("误差: %lf\n", fabs(I2 - I1) / 15);
    printf("调用函数次数: %d\n", (n - 1) * 2);
    return 0;
}

其中，`f` 函数为被积函数，`Simpson` 函数为复化Simpson公式，`main` 函数为主函数，实现输入积分区间和误差限，输出积分值、误差和调用函数次数。