数值计算方法实验三

最新推荐文章于 2021-12-23 18:43:20 发布

小wal

最新推荐文章于 2021-12-23 18:43:20 发布

阅读量1.7k

点赞数 3

分类专栏：数值计算方法文章标签：数值计算方法龙贝格算法变步长梯形算法

本文链接：https://blog.csdn.net/qq_41106517/article/details/94720408

版权

本实验中，我们应用变步长梯形法和龙贝格算法来求解积分问题，以达到10^-4的精度。对于变步长梯形法，用户需输入精度ε和区间[a, b]，程序将逐步计算并显示每一步的结果。而对于龙贝格算法，我们需要绘制其实现框图，并由用户输入精度和区间。计算过程中，将依次打印梯形序列、辛普森序列、柯特斯序列和龙贝格序列，最后输出满足精度要求的积分近似值。" 132956481,20014998,DOM操作：append与appendChild的区别,"['前端开发', 'javascript', 'HTML']

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

1.用变步长梯形法求下列积分,使精确度达到10-4。

要求：

(1) 精度ε和区间端点a，b的信息从键盘终端输入；

(2) 打印输出每一步的计算结果。

源代码：
#include<iostream>
#include<cmath>
using namespace std;
int main(){
	double x,a,s,b,e,h,T1,T2;
	cout<<"请输入区间端点a,b和精度e:"<<endl; 
	cin>>a>>b>>e;
	cout<<endl;
	h = b-a;
	T1 = (h/2)*(4/(1+pow(a,2))+4/(1+(pow(b,2))));
	int count=0;
	cout<<"K    "<<"  T1   "<<"     T2  "<<endl;
	do{
		if(count!=0){
			h  = h/2;
			T1 = T2;
		}
		count++;
		s = 0;
		x = a + h/2;
		do{
			s = s + (4/(1+pow(x,2)));
			x = x + h;
		}while(x<b);
		T2 = (T1/2)+h*s/2;
		printf("%d:  %.6lf   %.6lf\n",count,T1,T2); 
	}while(abs(T2-T1)>=e);
	cout<<endl;
	printf("最后结果为：%.4lf",T2);
	return 0;
}

2. 用龙贝格算法计算下列积分值，使精确度达到10-4。