龙贝格积分-matlab通用程序_第1页
龙贝格积分-matlab通用程序_第2页
龙贝格积分-matlab通用程序_第3页
全文预览已结束

下载本文档

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

文档简介

龙贝格积分-matlab通用程序龙贝格积分-matlab通用程序龙贝格积分-matlab通用程序资料仅供参考文件编号：2022年4月龙贝格积分-matlab通用程序版本号：A修改号：1页次：1.0审核：批准:发布日期：数值分析课程的大作业，教材《数值分析》李乃成.梅立泉clear

clc

formatlong

f=input('请输入原函数f=','s');

a=input('积分下限a=');

b=input('积分上限b=');

eps1=input('精度eps1=');

T(1)=double((b-a)/2*(limit(sym(f),findsym(sym(f)),a)+limit(sym(f),findsym(sym(f)),b)));

fork=2:4

sum1=0;

fori=1:2^(k-2)

sum1=sum1+subs(sym(f),findsym(sym(f)),(a+(2*i-1)*(b-a)/2^(k-1)));

end

T(k)=1/2*T(k-1)+(b-a)/(2^(k-1))*sum1;

end

fork=1:3

S(k)=T(k+1)+1/(4-1)*(T(k+1)-T(k));

end

fork=1:2

C(k)=S(k+1)+1/(4^2-1)*(S(k+1)-S(k));

end

R(1)=C(2)+1/(4^3-1)*(C(2)-C(1));

k=3;

while1

T(1)=T(2);

T(2)=T(3);

T(3)=T(4);

sum2=0;

fori=1:2^k

sum2=sum2+subs(sym(f),findsym(sym(f)),(a+(2*i-1)*(b-a)/2^(k+1)));

end

T(4)=1/2*T(4)+(b-a)/2^(k+1)*sum2;

S(1)=S(2);

S(2)=S(3);

S(3)=T(4)+1/(4-1)*(T(4)-T(3));

C(1)=C(2);

C(2)=S(3)+1/(4^2-1)*(S(3)-S(2));

R(2)=C(2)+1/(4^3-1)*(C(2)-C(1));

ifabs(R(2)-R(1))<eps1

人人文库> 全部分类> 应用文书 > 项目管理

温馨提示

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

最新文档

评论

0/150

提交评论