[算法][求积分][复合辛普森公式]
1 //这里f()为被积函数,输入a,b为积分上下限, 2 //eps为计算精度[这里要注意假收敛,一般设小一点好] 3 #include <iostream> 4 #include <cmath> 5 #define eps 1e-6 6 using namespace std; 7 double f(double x){ 8 return sqrt(1+cos(x)*cos(x)); 9 }//被积函数 10 double Sn(double a,double b,double n){ 11 double h=(a+b)/(2*n); 12 double sum=0; 13 for(int k=1;k<=n;k++){ 14 sum+=2*f(a+(2*k-1)*h); 15 sum+=f(a+2*k*h); 16 } 17 sum=(sum*2+f(a)-f(b))*h/3; 18 return sum; 19 }//辛普森公式 20 double Solve(double a,double b){ 21 int i=1; 22 double S1n=Sn(a,b,i); 23 double S2n=Sn(a,b,2*i); 24 while(fabs(S1n-S2n)>eps){ 25 cout<<"n = "<<i<<" S1 = "<<S1n<<" S2 = "<<S2n<<'\n'; 26 i++; 27 S1n=S2n; 28 S2n=Sn(a,b,2*i); 29 } 30 return S2n; 31 } 32 int main(){ 33 cout<<Solve(0,48)<<'\n'; 34 }
分类:
算法_算法与数据结构
标签:
复合辛普森
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· go语言实现终端里的倒计时
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· 周边上新:园子的第一款马克杯温暖上架
· Open-Sora 2.0 重磅开源!
· 分享 3 个 .NET 开源的文件压缩处理库,助力快速实现文件压缩解压功能!
· Ollama——大语言模型本地部署的极速利器
· DeepSeek如何颠覆传统软件测试?测试工程师会被淘汰吗?