数值分析龙贝格算法

《数值分析龙贝格算法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、《数值分析》课程设计报告龙贝格算法专业：班级：学号：姓名：1.实验目的运用龙贝格求积公式，编制龙贝格算法程序，计算积分的近似值，要求误差小于。2.程序流程图开始输入a,b,eb-a->h,h/2[f(a)+f(b)]->T1,1->k0->s,a+h/2->xs+f(x)->sx+h->xxT2T2+1/3(T2-T1)->S2k+1->k,h/2->hT2->T1,S2->S1k=1=S2+1/15(S2-S1)->C2C2->C1k=2=C2+1/63(C2-C1)->R2k=3R2->R1输出a,b,R2

2、R2-R1

3、

4、.程序代码#include"stdio.h"#include"math.h"intmain(void){floatf(float);floata,b,e,h,T1=0,T2=0,S1=0,S2=0,C1=0,C2=0,R1=0,R2=0,k,s,x;inti=0;printf("请输入积分的下限：");scanf("%f",&a);printf("请输入积分的上限：");scanf("%f",&b);printf("请输入允许误差：");scanf("%f",&e);k=1;h=b-a;T1=h*(f(a)+f(b))/2;printf("--------

5、--------------");printf("kT2S2C2R2");printf("%d%10.7f%10.7f%10.7f%10.7f",i,T1,S1,C1,R1);do{s=0;x=a+h/2;while(x

6、/63;C2=S2+(S2-S1)/15;C1=C2;k=k+1;h=h/2;T1=T2;S1=S2;}else{C2=S2+(S2-S1)/15;R2=C2+(C2-C1)/63;if(fabs(R2-R1)

7、rn0;}floatf(floatx){/*floaty=0;if(x==0.0)return1;y=(float)sin(x)/(x);returny;*/floaty;y=1/(sqrt(1-x*x));returny;}4运行结果