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

int main() {
    double a, b, eps;
    int n = 2, cnt = 0;
    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, &cnt);
        n *= 2;
    } while (fabs(I2 - I1) >= eps);
    printf("积分值: %lf\n", I2);
    printf("误差: %lf\n", fabs(I2 - I1) / 15);
    printf("函数调用次数: %d\n", cnt);
    return 0;
}

其中，`f` 函数为被积函数，`Simpson` 函数为复化Simpson公式，`main` 函数为主函数，实现输入积分区间和误差限，输出积分值、误差和函数调用次数。主要修改是增加了一个 `cnt` 变量，用于记录函数调用次数，并且在 `Simpson` 函数中通过指针方式将 `cnt` 传递进去，每次调用函数时增加 `cnt` 的值。