变步长梯形求积公式 c语言实现 数值积分

变步长梯形求积公式 c语言实现 数值积分

标签：计算方法实验

---

/*
    本实验用变步长梯形求积公式求x / (4 + x^2)在[0, 1]的定积分。
*/
#include <stdio.h>
#include <math.h>

double f(double x){
    return x / (4 + x * x);
}

int main(){
    double a = 0, b = 1, eps = 0.00001;  //上下限，精度
    double s, x, h = b - a, t1, t2 = h * (f(a) + f(b)) / 2, temp;

    do{
        s =