龙贝格求积公式 c语言实现 数值积分

龙贝格求积公式 c语言实现 数值积分

标签：计算方法实验

---

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

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

int main(){
    double a = 0, b = 1, eps = 0.00001;  //上下限，精度
    double s, x, s1, s2, c1, c2, r1, r2, h = b - a, t1 = h * (f(a) + f(b)) /