原创力文档- 1、本文档共9页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
数值分析实习报告--第1页
数值分析实习报告
姓名：***
学号：***
班级：***
数值分析实习报告--第1页
数值分析实习报告--第2页
前欧拉法
1.题目取h=0.1,利用Euler公式求解
dy/dx=y-2*x/y(0=x=1)
y(0)=1
2.思路
利用左矩形公式得到的公式yn+1=yn+hf(xn,yn)进行迭代，在进行迭代的过程使用for循
环，让n从1取到10，每迭代一次输出一个x与y，同时x加0.1再进行下一次迭代，一直
到循环结束。
3.程序
clear;y=1,x=0,%初始化
forn=1:10
y=1.1*y-0.2*x/y,x=x+0.1,
end
4.运行结果
y=1x=0
y=1.1000x=0.1000
y=1.1918x=0.2000
y=1.2774x=0.3000
y=1.3582x=0.4000
y=1.4351x=0.5000
y=1.5090x=0.6000
y=1.5803x=0.7000
y=1.6498x=0.8000
y=1.7178x=0.9000
y=1.7848x=1.0000
数值分析实习报告--第2页
数值分析实习报告--第3页
后欧拉法
1.题目取h=0.1,利用Euler公式求解
dy/dx=y-2*x/y(0=x=1)
y(0)=1
2.思路
利用右矩形公式得到的公式yn+1=yn+hf(xn+1,yn+1)进行迭代，在进行迭代的过程使用for
循环，让n从1取到10，每迭代一次输出一个x与y，同时x加0.1再进行下一次迭代，一
直到循环结束。
3.程序
clear;y0=1,x0=0,%初始化
forn=1:10
yp=y0+0.1*(y0-2*x0/y0);
x=x0+0.1;
x0=x
yp=y0+0.1*(yp-2*x/yp)
y0=yp
end
4.运行结果
y=1x=0
y=1.0918x=0.1000
y=1.1763x=0.2000
y=1.2546x=0.3000
y=1.3278x=0.4000
y=1.3964x=0.5000
y=1.4609x=0.6000
文档评论（0）