数值作业:龙贝格算法计算积分C语言实现

最新推荐文章于 2021-12-10 19:49:56 发布

WUST_陈迪洋

最新推荐文章于 2021-12-10 19:49:56 发布

阅读量1.5w

点赞数 4

分类专栏：数值分析文章标签：龙贝格算法计算积分算法 Romberg算法 c语言

本文链接：https://blog.csdn.net/Chen_dSir/article/details/71270316

版权

这篇博客介绍了如何使用C语言实现龙贝格算法来计算定积分，以解决数值计算问题。文中提到的例子是计算f(x) = sin(x)/x的积分，并提及在实现过程中遇到的一个关于中文空格的调试难题，强调了编程实践和调试的重要性。

摘要由CSDN通过智能技术生成

根据Romberg算法计算定积分，和变步长的Simpson算法的输入都一样．算法基本分析：输入a,b(积分上下限),n为积分区间等分数，eps为计算精度，我这里1e-7,代表0乘以10的负7次方．本题目取的例子为数值书137面的例子２，f(x)= sin(x)/x，下面给出代码：

/********************************************  
    > File Name: Dragon.c  
    > Author:chendiyang  
    > School:WUST_CST_1501班  
    > Myblog:www.chendsir.com  
    > Mail:1441353519@qq.com   
    > Created Time: 2017年05月6日 星期六 12时33分10秒  
 ************************************************************************/    
#include <stdio.h>
#include <math.h>
#define N 20
#define MAX 10     //数组存的最大行数
#define a 0.0000001    //积分下限
#define b 1.0      //积分上限
#define eps 1e-7 //精度

double f(double x)//所求积分公式
{
    return sin(x) / x;
}

double computeT(double aa, double bb, long int n)//复化梯形公式
{
    int i;
    double sum, h = (bb - aa) / n;
    for (i = 1; i < n; i++)
        sum &