实验8欧拉法_改进欧拉法_龙格-库塔法_线性多步法

《实验8欧拉法_改进欧拉法_龙格-库塔法_线性多步法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、西华大学计算机系学生上机实践报告西华数学与计算机学院上机实践报告课程名称：计算方法A年级：上机实践成绩：指导教师：严常龙姓名上机实践名称：解常微分方程初值问题学号：上机实践日期：2013.12.20上机实践编号：1上机实践时间：10:00-11:40一、目的1．通过本实验加深对欧拉法、改进欧拉法、龙格-库塔法、线性多步法的构造过程的理解；2．能对上述四种方法提出正确的算法描述编程实现，观察计算结果的改善情况。二、内容与设计思想自选常微分方程的初值问题，分别用欧拉法、改进欧拉法求解。分别用以上两种方法求解常微分方程初值问题：三、使用环境操作系统：WindowsXP软件环境：MicrosoftVi

2、sualC++四、核心代码及调试过程#include#include#definef(y)(y*y+1)#definem0.0//初值为0#defineh0.1//步长为0.1#definen14//迭代次数为14#definea0.0//定义区间长度#defined1.4voidgjol();//改进欧拉法voidol();//欧拉法main(){ol();printf("");gjol();}voidgjol(){第3页共3页西华大学计算机系学生上机实践报告inti;floaty[n+1];y[0]=m;//赋初值printf("改进欧拉法");f

3、or(i=0;i

4、f(y[0]);}else{y[i+1]=y[i-1]+2*h*f(y[i]);}printf("y[%d]=%f",i,y[i]);}第3页共3页西华大学计算机系学生上机实践报告printf("欧拉法结果为：");printf("y[%d]=%f",i,y[i]);}五、总结本次试验利用C语言实现了用欧拉法、改进欧拉法求常微分方程的初值问题,通过本实验加深了对欧拉法、改进欧拉法求常微分方程的初值问题的理解，对C语言在计算方法中的应用有了进一步的掌握。六、附录参考：《数值计算方法与算法》（第二版）第3页共3页