1、这样,前一次分割得到的函数值在分半以后仍可被利用,且易于编程。2.计算公式梯形公式复化辛普森公式复化科特斯公式 龙贝格求积公式其对应的公式为:T2n=1/2(Tn+Hn) (梯形公式)Sn=4/(4-1)T2n-1(4-1)Tn (辛普森公式)Cn=42/(42-1)S2n-1/(42-1)Sn (柯特斯公式)Rn=43/(43-1)C2n-1/(43-1)Cn (龙贝格求积法公式)3.算法描述3.1龙贝格算法基本描述先算出T0(0),从而计算出T0(1),以此类推,直到计算出|T0(0)-Tn-1(0)|e即可利用加速推算公式推算出结果。 3.2龙贝格算法程序包步骤 1.输入积分上限 2输入
2、程序下限3输入区间等分数4输入要求的函数5.算法程序表示#includemath.h#define A(x)(sin(x)/x) /宏定义若干常用函数A,B,C,D,E,G#define B(x)(cos(x*x+2*x+1)#define C(x)(atan(sqrt(x*x+1)#define D(x)(sqrt(exp(x)+sin(2*x)#define E(x)(x*x*x+3*x*x+5)#define G(x)(log10(x)/pow(2,x)double t20,s20,c20,r20;/定义全局数组double dh,fan,a,b,m; /定义全局变量int jj=0;c
3、har hs;double F(double x) /用switch调用若干被积函数 switch(hs) case A:fan=A(x);break;Bfan=B(x);Cfan=C(x);Dfan=D(x);Efan=E(x);Gfan=G(x); default :printf(输入错误!); return(fan);/返回被积函数值double H(int i) /求和函数并返回和SUM int j; double zh,SUM=0.0; /定义求和变量SUM并赋初值 for(j=1;j=pow(2,i-1);j+) zh=(a+(2*j-1)*(b-a)/pow(2,i); SUM=
4、SUM+F(zh); /调用F(x)函数 SUM=(b-a)*SUM/pow(2,i); return(SUM);double Txing(int k) /梯形公式 if(k=0) dh=tjj=(b-a)/2)*(F(a)+F(b); /分半次数为零时T形公式求积 else dh=0.5*Txing(k-1)+H(k); /Txing函数递归调用循环输出并返回dh t+jj=dh; m=pow(2,jj); printf(T%0.0lf=%0.7lft,m,tjj);/输出并返回dh return(dh);double Simpson(int k) /辛普森公式 int i,j; Txing
5、(k); /调用梯形公式 for(i=0;i=k-1;i+) si=(4.0*ti+1-ti)/3.0; /递推辛普森公式 printf(n for(j=0;j+) /循环输出 m=pow(2,j);S%0.0lf=%0.7lf,m,sj);t return(sk-1); /返回最后一个值sk-1double Cotes(int k) /科特斯公式 Simpson(k);=k-2; ci=(16.0*si+1-si)/15.0; /递推科特斯公式j+) /循环输出C%0.0lf=%0.7lft,m,cj); return(ck-2); /返回最后一个值ck-2double Romberg(in
6、t k) /隆贝格公式 /调用科特斯公式 Cotes(k);=k-3; ri=(64.0*ci+1-ci)/63.0; /递推隆贝格公式 m=pow(2,j);R%0.0lf=%0.7lft,m,rj); return(rk-3); /返回最后一个值rk-3main() int k; char y;请从以下公式中选择积分函数:A: sin(x)/(x)tB:cos(x2+2x+1)tC:atan(sqrt(x2+1)nnD:sqrt(ex+sin(2x)tE:x3+3*x2+5tG:log10(x)/pow(2,x)nn请选择函数F(x)(大写): scanf(%c,&hs); getchar
7、();请输入您选用的求积公式第一个字母(大写): y=getchar(); switch(y) /根据输入的大写字母判断所选求积公式 case T请输入分半次数: scanf(%dk); printf(请输入积分区间a,b且(ab): scanf(%lf,%lfa,&b); if(a=0) a=0.0000000001; /输入当积分区间为零时默认为一极小数 printf(-nnTxing=%0.7fn,Txing(k); /输出计算结果小数点后保留7位S请输入分半次数大于等于1的正整数: scanf(-nSimpson=%0.7fn,Simpson(k);请输入分半次数大于等于2的正整数: /输入当积分区间为零时默认为一极小数 Cotes=%0.7fn,Cotes(k);R请输入分半次数大于等于3的正整数:Romberg=%0.7fn,Romberg(k); default:公式输入有误,请查证!6.算法结果截图
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1