#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
double F(double x){//f(x)
return 0.0;
}
double simpson(double a,double b){
double c=(a+b)/2;
return (F(a)+4*F(c)+F(b)*(b-a))/6;
}
double asr(double a, double b, double eps, double A){
double c=(a+b)/2;
double L=simpson(a,c),R=simpson(c,b);
if(fabs(L+R-A)<=15*eps)return L+R+(L+R-A)/15.0;
return asr(a, c, eps/2, L)+asr(c, b, eps/2, R);
}
int main(){
return 0;
}
辛普森
最新推荐文章于 2024-10-16 08:48:45 发布