Romberg（龙贝格求积）

最新推荐文章于 2023-04-22 16:18:45 发布

原创最新推荐文章于 2023-04-22 16:18:45 发布 · 2.2k 阅读

1 ·

CC 4.0 BY-SA版权

文章标签：

#C语言 #龙贝格

计算方法专栏收录该内容

13 篇文章

订阅专栏

本文展示了一个使用C语言实现的简单数值积分程序。该程序通过复合梯形法则近似计算了一个特定函数在指定区间上的积分值，并展示了如何通过逐步减小区间宽度来提高积分精度。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include <stdio.h>
float f(float x)
{
    return(exp(x)/(4+x*x));
}
main()
{
    float a=0,b=1;//积分上下界
    float h=b-a,T1,T2,s,x;//h是区间长度，T1是分为n等分的近似值，T2是分为2n等分的近似值
    int i;//次数
    T1=h/2*(f(a)+f(b));//公式
    for(i=0;i<4;i++)
    {
        printf("h=%f,T=%f\n",h,T1);
        s=0;
        x=a+h/2;
        while(x<b){
            s+=f(x);
            x+=h;
        }
        T2=T1/2+h/2*s;//公式
        h/=2;//将区间逐次分半
        T1=T2;//赋值后继续循环
    }
    printf("h=%f,T=%f\n",h,T1);
    return;
}