数值积分 (辛普森公式辛普森自适应法则) UVA

最新推荐文章于 2023-09-23 10:50:34 发布

BryantMark

最新推荐文章于 2023-09-23 10:50:34 发布

阅读量1.4k

点赞数

分类专栏： acm 数值积分数学

本文链接：https://blog.csdn.net/fanzhengchen234/article/details/8841671

版权

该博客详细介绍了如何解决UVA 1356编程问题，重点关注数值积分的方法，特别是辛普森公式及其自适应法则的应用。通过对算法的深入探讨，读者将理解如何在实际问题中有效地实施这些数学工具。

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

double F(double x) {
}
double simpson(double l, double r) {
	double m = (l + r) * 0.5;
	return (F(l) + 4 * F(m) + F(r)) * (r - l) / 6.0;
}
double asr(double l, double r, double eps, double A) {
	double m = (l + r) * 0.5;
	double L = simpson(l, m), R = simpson(m, r);
	if (fabs(L + R - A) <= 15 * eps)
		return L + R + (L + R - A) / 15.0;
	return asr(l, m, eps * 0.5, L) + asr(m, r, eps * 0.5, R);
}
double asr(double l, double r, double eps) {
	return 2.0 * asr(l, r, eps, simpson(l, r));
}
int dblcmp(double x) {
	if (x < -eps)
		return -1;
	return x > eps;
}

UVA 1356

#include<cstdio>
#include<cstring>
#include<queue>
#include<vector>
#include<iostream>
#include<string>
#include<sstream>
#include<cctype>
#include<set>
#include<sta