龙贝格积分法C语言(西安交大)(共2页)_第1页
龙贝格积分法C语言(西安交大)(共2页)_第2页
龙贝格积分法C语言(西安交大)(共2页)_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、精选优质文档-倾情为你奉上#include<stdio.h>#include<math.h>void main() int i; float a,b,fa,fb; double T,T1,S,S1,C,C1,R,R1,e;double fx1(double x);printf("请输入积分下界 a=n"); /*输入X的下界*/ scanf("%f",&a); printf("请输入积分下界 b=n"); /*输入X的上界*/ scanf("%f",&b);printf(&q

2、uot;请输入计算精度 e=n"); scanf("%lf",&e); fa=fx1(a); /*计算fx函数上下界的值*/ fb=fx1(b); printf("a=%3.7fnb=%3.7fnfa=%3.7fnfb=%3.7fne=%3.10f",a,b,fa,fb,e); /*输出上、下界值，上下界fx函数值,计算精度e*/ printf("n"); /*程序开始*/ T1=(b-a)*(fa+fb)/2; for(i=0;i+) double fx2(double i,double a,double b);

3、double Tkk; Tkk=fx2(i,a,b); /*调用fx函数*/ T=T1/2+Tkk; T1=T; S=T+(T-T1)/(pow(4,1)-1); S1=S; if(i>=1) C=S+(S-S1)/(pow(4,2)-1); C1=C; if(i>=2) R=C+(C-C1)/(pow(4,3)-1); if(abs(R-R1)<0.) /*判断是否满足计算精度，满足则退出迭代，不满足则继续*/ break; R1=R; printf("积分计算数值为n"); /*输出计算结果*/printf("%2.7fn",T);double fx1(double x) /*编写一个被积函数fx*/ double fx; fx=1/(1+x); return(fx); double fx2(double i,double a,double b) /*编写T的高次迭代函数*/ /double fx1(double z); double m,l; double Tk,Tkk,z; m=pow(2,i);Tk=0;for(l=1;l<=m;l+) z=a+(2*l-1)*(b-a)

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论