欢迎来到天天文库
浏览记录
ID:14295419
大小:73.93 KB
页数:5页
时间:2018-07-27
《数值分析龙贝格算法》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
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(x6、/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运行结果
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运行结果
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运行结果
此文档下载收益归作者所有