复化梯形公式.docVIP

复化梯形公式 #includeiostream.h #includemath.h double function1(double x)//被积函数 { double s; s=x/(4+x*x); return s; } double function2(double x)//被积函数 { double s; s=sqrt(x); return s; } double ReiterationOfTrapezoidal(double a,double b,double n,double f(double x))//复化梯形公式 { double h,fa,fb,xk,t; h=(b-a)/n;//步长 fa=f(a); fb=f(b); double s=0.0; for(int k=1;kn;k++) { xk=a+k*h; s=s+f(xk); } t=h*(fa+fb)/2+h*s;//复化梯形公式 return t; } double ReiterationOfSimpson(double a,double b,double n,double f(double x))//复化Simpson公式 { double h,fa,fb,xk,xj; h=(b-a)/n; fa=f(a); fb=f(b); double s1=0.0; double s2=0.0; for(int k=1;kn;k++) { xk=a+k*h; s1=s1+f(xk); } for(int j=0;jn;j++) { xj=a+(j+0.5)*h; s2=s2+f(xj); } double sn;//和 sn=h/6*(fa+fb+2*s1+4*s2);//复化Simpson公式 return sn; } main() { double a,b,Result,n; /*第一小题*/ cout请输入第一小题积分下限:endl; cina; cout请输入第一小题积分上限:endl; cinb; cout请输入第一小题分割区间数n:endl; cinn; cout第一小题复化梯形公式计算结果:; Result=ReiterationOfTrapezoidal(a, b, n,function1); coutResultendl; cout第一小题复化Simpson公式计算结果:; Result=ReiterationOfSimpson(a, b, n,function1); coutResultendl; /*第二小题*/ cout请输入第二小题积分下限:endl; cina; cout请输入第二小题积分上限:endl; cinb; cout请输入第二小题分割区间数n:endl; cinn; cout第二小题复化梯形公式计算结果:; Result=ReiterationOfTrapezoidal(a, b, n,function2); coutResultendl; cout第二小题复化Simpson公式计算结果:; Result=ReiterationOfSimpson(a, b, n,function2); coutResultendl; return 0; } 本文来自CSDN博客，转载请标明出处：/liuhongyi666/archive/2008/12/22/3584482.aspx #include conio.h /* 此头函数请不要删除 */ #include stdio.h float N,A,B,X; float F(float x); float Ti() //复化梯形公式 { float i,n,h,t,g=0; n=N; h=(B-A)/n; for(i=1;in;i++) {X=A+i*h; g=F(X)+g; } t=h*(F(A)+2*g+F(B))/2; return(t); } float Si() //复化辛普生公式 {float i,n,h,t,p=0,g=0; n