数值分析-常微分数值解的求法C语言.docVIP

本科生课程设计报告 实习课程 数值分析 学院名称 管理科学学院 专业名称 信息与计算科学 学生姓名 学生学号 指导教师 实验地点 实验成绩 二〇 一 六 年 六 月 二〇 一 六 年 六 摘要 常微分方程数值解法计算数学的一个分支.是解常微分方程各类定解问题的数值方法.现有的解析方法只能用于求解一些特殊类型的定解问题，实用上许多很有价值的常微分方程的解不能用初等函数来表示，常常需要求其数值解.所谓数值解，是指在求解区间内一系列离散点处给出真解的近似值.这就促成了数值方法的产生与发展.? 实验内容和要求 常微分方程初值问题 有精确解。 要求：分别取步长h = 0.1，0.01，0.001，采用改进的Euler方法、4阶经典龙格-库塔R－K方法和4阶Adams预测-校正方法计算初值问题。以表格形式列出10个等距节点上的计算值和精确值，并比较他们的计算精确度。其中多步法需要的初值由经典R-K法提供。运算时，取足以表示计算精度的有效位。 算法说明 2.1函数及变量说明 表1 函数及变量定义 函数/变量名 类 型 说 明 f double 要求的微分方程 p_f double 微分方程原函数 Euler，Euler_Pro K_R,Adams double