[算法][求积分][复合辛普森公式]

 

 

复制代码
 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 }
复制代码

 

posted @   beautifulzzzz  阅读(2912)  评论(0编辑  收藏  举报
编辑推荐:
· go语言实现终端里的倒计时
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
阅读排行:
· 周边上新:园子的第一款马克杯温暖上架
· Open-Sora 2.0 重磅开源!
· 分享 3 个 .NET 开源的文件压缩处理库,助力快速实现文件压缩解压功能!
· Ollama——大语言模型本地部署的极速利器
· DeepSeek如何颠覆传统软件测试?测试工程师会被淘汰吗?
点击右上角即可分享
微信分享提示