《8_龙贝格算法》由会员分享,可在线阅读,更多相关《8_龙贝格算法(3页珍藏版)》请在金锄头文库上搜索。
实验八龙贝格算法代码:#include#includeusingnamespacestd;intmain()intk=1;doublea,b,c,d,h,x,s,T2,S2,C2,R2;coutvv请输入a和b的值;cinab;coutc;coutvvendl;h=b-a;if(a=0)T0=(1+sin(b)/b)*h/2;elseT0=(sin(a)/a+sin(b)/b)*h/2;loop2:s=0;x=a+h/2;loop1:s=s+sin(x)/x;x+=h;if(xvb)gotoloop1;T1=T0/2+s*h/2;S1=T1+(T1-T0)/3;if(k=1)h=h/2;k+;T0=T1;S0=S1;gotoloop2;elseC1=S1+(S1-S0)/15;if(k=2)C0=C1;k+;h=h/2;T0=T1;S0=S1;gotoloop2;elseR1=C1+(C1-C0)/63;if(k=3)R0=R1;C0=C1;h=h/2;k+;T0=T1;S0=S1;gotoloop2;elsed=fabs(R1-R0);if(c=d)R0=R1;C0=C1;h=h/2;k+;T0=T1;S0=S1;gotoloop2;elsecout得出的值是vR1vendl;return0;流程图:运行界面