谁有关于定积分的算法阿,谢谢,急用 谁有关于定积分的算法阿,谢谢,急用 解决方案 » 免费领取超大流量手机卡,每月29元包185G流量+100分钟通话, 中国电信官方发货 #include <iostream.h>#include <iomanip.h>#include<math.h>double f(double x);void main(){ int n=1; //初值 double a=0,b=1; double h,Tn,T2n,In,I2n; const double eps=1e-8; h=b-a; T2n=I2n=h*(f(a)+f(b))/2; In=0; while(fabs(I2n-In)>=eps) //求积分 { Tn=T2n; In=I2n; double sigma=0.0; for(int k=0;k<n; K++) //求变步长梯形的和部分 { double x=a+(K+0.5)*h; sigma+=f(x); } T2n=(Tn+h*sigma)/2.0; //变步长梯形 I2n=(4*T2n-Tn)/3.0 //辛普生公式 n*=2; //划分 h/=2; } cout<<"the integral of f(x) from" <<a<<"to" <<b<<"is \n" <<setiosflags(ios::fixed) <<setprecision(8) <<setw(10)<<I2n<<endl; //输出结果} double f(double x) { return exp(x)/(1+x*x); } #include <iostream.h>#include <iomanip.h>#include<math.h>double f(double x);void main(){ int n=1; //初值 double a=0,b=1; double h,Tn,T2n,In,I2n; const double eps=1e-8; h=b-a; T2n=I2n=h*(f(a)+f(b))/2; In=0; while(fabs(I2n-In)>=eps) //求积分 { Tn=T2n; In=I2n; double sigma=0.0; for(int k=0;k<n; K++) //求变步长梯形的和部分 { double x=a+(K+0.5)*h; sigma+=f(x); } T2n=(Tn+h*sigma)/2.0; //变步长梯形 I2n=(4*T2n-Tn)/3.0 //辛普生公式 n*=2; //划分 h/=2; } cout<<"the integral of f(x) from" <<a<<"to" <<b<<"is \n" <<setiosflags(ios::fixed) <<setprecision(8) <<setw(10)<<I2n<<endl; //输出结果} double f(double x) { return exp(x)/(1+x*x); } java 读取oracle库数值型“空值”字段为何变为0.0,寻求帮助! URL相对转绝对的问题 J2ME/J2SE/J2SE,在编程技术有什么区别? 运行JSP要设置哪些 为什么没有打印出异常呢 小妹刚学java 照书上抄的代码怎么会有NullPointerException呢 java GUI设计中,多按钮响应的问题 郁闷 散分~~~ 如何得到jbutton的key event? 有人知道这个用java怎么写嘛 。。。我是一个萌新 applet显示问题 一直没有解决好的问题:如果判断一个数据库的连接是否正常
#include <iomanip.h>
#include<math.h>double f(double x);void main(){
int n=1; //初值
double a=0,b=1;
double h,Tn,T2n,In,I2n;
const double eps=1e-8;
h=b-a;
T2n=I2n=h*(f(a)+f(b))/2;
In=0; while(fabs(I2n-In)>=eps) //求积分
{
Tn=T2n;
In=I2n;
double sigma=0.0;
for(int k=0;k<n; K++) //求变步长梯形的和部分
{
double x=a+(K+0.5)*h;
sigma+=f(x);
}
T2n=(Tn+h*sigma)/2.0; //变步长梯形
I2n=(4*T2n-Tn)/3.0 //辛普生公式 n*=2; //划分
h/=2;
}
cout<<"the integral of f(x) from" <<a<<"to" <<b<<"is \n"
<<setiosflags(ios::fixed)
<<setprecision(8)
<<setw(10)<<I2n<<endl; //输出结果
} double f(double x)
{
return exp(x)/(1+x*x);
}
#include <iomanip.h>
#include<math.h>double f(double x);void main(){
int n=1; //初值
double a=0,b=1;
double h,Tn,T2n,In,I2n;
const double eps=1e-8;
h=b-a;
T2n=I2n=h*(f(a)+f(b))/2;
In=0; while(fabs(I2n-In)>=eps) //求积分
{
Tn=T2n;
In=I2n;
double sigma=0.0;
for(int k=0;k<n; K++) //求变步长梯形的和部分
{
double x=a+(K+0.5)*h;
sigma+=f(x);
}
T2n=(Tn+h*sigma)/2.0; //变步长梯形
I2n=(4*T2n-Tn)/3.0 //辛普生公式 n*=2; //划分
h/=2;
}
cout<<"the integral of f(x) from" <<a<<"to" <<b<<"is \n"
<<setiosflags(ios::fixed)
<<setprecision(8)
<<setw(10)<<I2n<<endl; //输出结果
} double f(double x)
{
return exp(x)/(1+x*x);
}