发布时间：
本文归纳常见常微分方程的解析解解法以及基于Python的微分方程数值解算例实现。
考虑常微分方程的解析解法，我们一般可以将其归纳为如下几类：
这类微分方程可以变形成如下形式： 两边同时积分即可解出函数，难点主要在于不定积分，是最简单的微分方程。
某些方程看似不可分离变量，但是经过换元之后，其实还是可分离变量的，不要被这种方程迷惑。
形如 的方程叫做一阶线性微分方程，若 为0，则方程齐次，否则称为非齐次。
解法： (直接套公式)
伯努利方程 形如 的方程称为伯努利方程，这种方程可以通过以下步骤化为一阶线性微分方程： 令 , 方程两边同时乘以 ，得到 即 . 这就将伯努利方程归结为可以套公式的一阶线性微分方程。
形如 的方程称为二阶常系数微分方程，若 ，则方程称为齐次的，反之称为非齐次的。以下默认方程是非齐次的。
求解此类方程分两步：
原方程的解 = 齐次通解 + 非齐次特解
首先假设 .用特征方程法，写出对应的特征方程并且求解： 解的情况分为以下三种：
情况一：方程有两个不同的实数解
假设两个实数解分别是 , 此时方程的通解是 情况二：方程有一个二重解 假设该解等于 ，此时方程的通解是 情况三：方程有一对共轭复解 假设这对解是 , 此时方程的通解是
对于 和特征根的情况，对特解的情况做如下归纳：
形如 的方程叫做高阶常系数微分方程，若 ，则方程是齐次的，否则是非齐次的。下面默认方程是非齐次的。
求解此类方程分两步：
原方程的解 = 齐次通解 + 非齐次特解
考虑带有第三类边界条件的二阶常系数微分方程边值问题
问题一:两点边值问题的解析解
由于此方程是非齐次的,故 求解此类方程分两步：
原方程的解 = 齐次通解 + 非齐次特解
首先假设 . 用特征方程法，写出对应的特征方程 求解得到两个不同的实数特征根: .
此时方程的齐次通解是
由于 . 所以非齐次特解形式为 将上式代入控制方程有 求解得: , 即非齐次特解为 .
原方程的解 = 齐次通解 + 非齐次特解
于是,原方程的全解为 因为该问题给出的是第三类边界条件,故需要求解的导函数 且有 将以上各式代入边界条件 解此方程组可得: .
综上所述,原两点边值问题的解为
对一般的二阶微分方程边值问题 假定其解存在唯一, 为求解的近似值, 类似于前面的做法,
考虑带有第三类边界条件的二阶常系数微分方程边值问题
问题二:有限差分方法算出其数值解及误差 对于 原问题 ,取步长 h= ,用 有限差分 求其 近似解 ,并将结果与 精确解y(x)=-x-1 进行比较.
因为
先以将区间划分为5份为例，求出数值解
结果:
是不是解出数值解就完事了呢？当然不是。我们可以将问题的差分格式解与问题的真解进行比较，以得到解的可靠性。通过数学计算我们得到问题的真解为 ，现用范数来衡量误差的大小：
结果:
接下来绘图比较 时数值解与真解的差距：
结果:
将区间划分为 份, 即 时.
结果:
绘图比较 时数值解与真解的差距：
最后，我们还可以从数学的角度分析所采用的差分格式的一些性质。因为差分格式的误差为 ， 所以理论上来说网格每加密一倍，与真解的误差大致会缩小到原来的 . 下讨论网格加密时的变化：
结果：
我只会前两种方法，即欧拉法和龙哥库塔法。后面的不是很熟悉
计算方法又称数值分析。是为各种数学问题的数值解答研究提供最有效的算法，计算方法主要内容包括函数逼近论、数值微分、数值积分、误差分析等，常用方法有迭代法、差分法、插值法、有限元素法等，现代计算方法要求适应电子计算机的特点。
误差与原则误差种类模型误差、观测误差、截断误差和舍入误差，法则加减运算近似数加减时，把其中小数位数较多的数四舍五入，使其比小数位数最少的数多一位小数，计算保留的小数位数与原近似数最小数位数最少者相同。
乘除运算近似数乘除时，各因子保留位数应比小数位数最少的数多一位小数，计算保留的小数位数与原近似数最小数位数最少者位数至多少一位，乘方与开方运算近似数乘方与开方时，计算保留的小数位数与原近似数位数相同，对数运算近似数对数时，计算保留的小数位数与原近似数位数相同，注意避免两个相近的数相减，避免除数绝对值远远小于被除数绝对值的除法，避免大数吃掉小数，计算讲效率，尽可能减少运算。
计算方法的特点
插值方法Lagrange插值线性插值、抛物线插值，Newton插值，分段插值，Hermite插值，分段三次Hermite插值，三次样条插值，最小二乘法直线拟合与多项式拟合，数值积分机械求积法梯形公式、中矩形公式、Simpson公式，Newton-Cotes求积法，复化求积法复化梯形公式、复化Simpson公式、复化Cotes公式，Romberg求积法，Guass求积法，数值微分求积法。
常微分方程的数值解法尤拉方法尤拉法、隐式尤拉法、二步尤拉法，改进尤拉方法，龙格-库塔方法，线性多步法亚当姆斯方法， 方程求根的数值解法二分法，迭代法，埃特金法，牛顿法牛顿下山法，近似牛顿法简化牛顿法、弦截法抛物线法，线性方程组的解法高斯消去法顺序消去法、列主元消去法、全主元消去法，矩阵三角分解法，追赶法平方根法，范数，简单迭代法Jacobi迭代法，Gauss-Seidel迭代法。
随机环境中经济增长模型研究广义生产函数假设下的经济增长模型分析考虑市场预期的供求关系模型基于Matlab的离散事件模拟用风险预算进行资产配置有向图上的PAR贯序模拟系统单圈图的一般Randic指标的极值问题模糊数学在公平评奖问题中的应用模糊矩阵在环境评估中的初步应用模糊评判在电脑中的初步应用数学家的数学思想Riemann积分定义的网收敛表述微积分思想在不等式证明中的应用用有限的尺度标量无限的过程-略论极限ε语言在微积分及现代数学中的位置及意义微积分思想在几何问题中的应用齐次平衡法求KdV-Burgers方程的Backlund变换Painleve分析法判定MKdV-Burgers方程的可积性直接法求KdV-Burgers方程的对称及精确解行波求解KdV-Burgers方程因子有向图的矩阵刻划简单图上的lit-only sigma-game半正则图及其线图的特征多项式与谱分数有向图的代数表示WWW网络的拓扑分析作者合作网络等的拓扑分析古诺模型价格歧视用数学软件做计算微分方程的计算器用数学软件做矩阵计算的计算器弹簧－质点系统的反问题用线性代数理论做隐含语义搜索对矩阵若当标准型理论中变换阵求法的探讨对矩阵分解理论的探讨对矩阵不等式理论的探讨(1)对矩阵不等式理论的探讨(2)函数连续性概念及其在现代数学理论中的延伸从有限维空间到无限维空间Banach空间中脉冲泛函微分方程解的存在性高阶脉冲微分方程的振动性具有积分边界条件的分数阶微分方程解的存在唯一性分数阶微分方程的正则摄动一个形态形成模型的摄动解一个免疫系统常微分方程模型的渐近解前列腺肿瘤连续性激素抑制治疗的数学模型前列腺肿瘤间歇性激素抑制治疗的数学模型病毒动力学数学模型肿瘤浸润数学模型耗散热方程初边值问题解的正则性耗散波方程初边值问题解的正则性耗散Schrodinger方程初边值问题解的正则性非线性发展方程解得稳定性消费需求的鲁棒调节生产函数的计量分析企业的成本形态分析的研究分数阶Logistic方程的数值计算分数阶捕食与被捕食模型的数值计算AIDS传播模型的全局性分析HIV感染模型的全局性分析风险度量方法的比较及其应用具有区间值损益的未定权益定价分析模糊规划及其在金融分析中的应用长依赖型金融市场股票价格与长相依性分数布朗运动下的外汇期权定价不确定性与资产定价加油站点的分布与出租车行业的关系
论文发表写作指导：
有一个网站叫中华论文中心，貌似有很多文章，你自己上去看下吧！
任找一本数学手册查（我是认真的！！！！）
举例说明常微分方程模型是各类数学建模竞赛中常见的模型, 并通过列举一些参考文献来说明此类模型的建模方法和求解求解技巧不仅相同. 从而得出"常微分方程在数学建模中的应用"是值得研究的.
1500字太夸张了，给你一下提示吧！ 1、运用微分方程或微分方程组，可以描述经济系统的动态运行规律。2、运用微分方程，可以分析经济系统的均衡与稳定性。3、在微分方程中加入控制变量，将经济学问题转化为最优控制问题，可以分析经济系统的最优控制策略。目前比较常用的微分方程在经济学中的应用有：（1）最早的哈罗德-多马经济增长模型、索罗模型等均属于微分方程（或转化为差分方程）模型。（2）后来的经济增长的世代交替模型等也是运用的微分方程。（3）技术扩散的巴斯模型，以及分析竞争洛克塔-瓦塔利亚模型也是微分方程模型。（4）亚瑟的路径依赖与锁定模型是随机微分方程。（5）布莱克-斯科尔斯期权定价模型，源于随机微分方程和变分法。（6）各种进化博弈模型中的复制动态方程是微分方程。
要的话请联系我邮箱（点我可见）。13 【篇名】 偏微分方程组的对称群及其在弹性力学方程组中应用 CAJ原文下载 PDF原文下载 【作者】 张鸿庆. 朝鲁. 唐立民. 【刊名】 大连理工大学学报 1997年03期 编辑部Email 《中文核心期刊要目总览》来源期刊 “中国期刊方阵”入选期刊 ASPT来源刊 CJFD收录期刊 【机构】 大连理工大学数学科学研究所. 大连理工大学工程力学研究所. 【关键词】 偏微分方程. 弹性力学. 对称群／不变向量场. 符号运算. 【聚类检索】 同类文献 引用文献 被引用文献 【摘要】 给出了非退化线性偏微分方程组及二次型泛函对称群的不变向量场的一般形式和一类特殊形式非线性偏微分方程组对称群的简化计算条件；利用以上结论及作者以往工作，借助符号运算语言MathematicaTM计算了平面弹性力学方程组一阶Lie－Bactlund对称群的不变向量场，以及应力函数对应的三维弹性力学方程组的Lie代数．为构造弹性力学方程组的一类广泛精确解及守恒律提供了必要的基础，并说明了结论对计算偏微分方程组对称群时的简化作用 【光盘号】 SCTC9706 14 【篇名】 力学中一类变系数微分方程可调参数模型解法 CAJ原文下载 PDF原文下载 【作者】 赵文福. 封营儒. 连星耀. 黎明安. 【刊名】 西安理工大学学报 1995年02期 编辑部Email CJFD收录期刊 【机构】 西安理工大学机械工程系. 【关键词】 可调参数. 变系数微分方程. 非均匀控制参数. 【聚类检索】 同类文献 引用文献 被引用文献 【摘要】 结合一种非均匀控制参数，提出了一种变系数微分方程的可调整参数模型解法，可以很方便地处理由于物理上、几何上的非均匀、非线性而导致数学上的变系数微分方程，应用这种模型可以用非常少的单元得到较满意的数值结果。 【光盘号】 SCTC9508 31 【篇名】 材料力学弯曲问题中集中量与分布量的统一处理 CAJ原文下载 PDF原文下载 【作者】 周锡勤. 张存道. 【刊名】 现代电力 1995年02期 编辑部Email CJFD收录期刊 【机构】 北京动力经济学院. 【关键词】 集中量. 分布量. 弯曲变形. 【聚类检索】 同类文献 引用文献 被引用文献 【摘要】 介绍了利用δ函数统一处理集中量与分布量的一般方法。着重讨论了这种方法在建立含集中量的杆件弯曲时的平衡微分方程的应用，从而推广了材料力学中杆件弯曲时的平衡微分方程。该方程更全面更精确地反映了杆件弯曲这一物理现象。作者把它称为梁弯曲时的广义平衡微分方程。 【光盘号】 SCTC95S5 38 【篇名】 双相材料空间中平片界面裂纹问题的超奇异积分-微分方程 CAJ原文下载 PDF原文下载 【作者】 乐金朝. 汤任基. 【刊名】 科学通报 1996年15期 编辑部Email 《中文核心期刊要目总览》来源期刊 “中国期刊方阵”入选期刊 ASPT来源刊 CJFD收录期刊 【机构】 郑州工学院道路检测与CAE技术研究中心. 上海交通大学工程力学系 郑州 450002 . 上海 200030. 【关键词】 双相材料. 平片界面裂纹. 超奇异积分-微分方程. 【聚类检索】 同类文献 引用文献 被引用文献 【摘要】 <正> 随着复合材料的广泛应用,界面断裂力学成为国际断裂界的前沿研究课题,该领域的研究工作引起了国内外力学家、金属物理学家及材料科学家的广泛关注,并取得了许多新进展。据作者所知,目前的工作主要是研究二维问题,由于数学和力学等方面的困难,三维界面断裂力学方面的研究工作报道较少。本文利用双相材料空间在集中力作用下的弹性力学基本解,使用边界元法,在有限部积分的意义下将任意形状的平片界面裂纹问题归结为一组以裂纹面上的位移间断为未知函数的超奇异积分-微分方程。此组方程对于进一步开展三维界面断裂力学问题的研究具有重要意义。 【光盘号】 SCTA96S4 39 【篇名】 常微分方程的不变式在量子力学中的应用 CAJ原文下载 PDF原文下载 【作者】 杨进. 【刊名】 大学物理 1998年08期 编辑部Email 《中文核心期刊要目总览》来源期刊 CJFD收录期刊 【机构】 成都气象学院基础科学系. 【关键词】 常微分方程. 不变式. 库仑场. 【聚类检索】 同类文献 引用文献 被引用文献 【摘要】 利用常微分方程的不变式，非常方便地求解了一些量子力学问题． 【光盘号】 SCTA9809 40 【篇名】 保守力系的变形拉格朗日方程及其应用 CAJ原文下载 PDF原文下载 【作者】 梁志强. 【刊名】 泰安师专学报 2000年06期 编辑部Email CJFD收录期刊 【机构】 泰安师专物理系!山东泰安271000. 【关键词】 Lagrandge方程. 轨道微分方程. 轨道方程. 【聚类检索】 同类文献 引用文献 被引用文献 【摘要】 从保守力系的拉格朗日方程出发 ,导出一种用于求解保守系统轨道微分方程的变形拉格朗日方程。并将其应用于有心力问题及抛体问题 ,导出了有心力问题的轨道微分方程Binet公式及抛体轨道方程。保守力系的变形拉格朗日方程提供了求解运动物体轨道方程的新方法 ,同时也丰富了分析力学的教学内容。 【光盘号】 SOCI0105
在一篇数学 教育 论文中，题目是论文的要件之首，它不同于一般 文章 的题目，我们要重视题目的重要性。以下是我为大家精心准备的数学教育论文题目，欢迎阅读!数学教育论文题目(一) 1、浅谈中学数学中的反证法 2、数学选择题的利和弊 3、浅谈计算机辅助数学教学 4、数学研究性学习 5、谈发展数学思维的 学习 方法 6、关于整系数多项式有理根的几个定理及求解方法 7、数学教学中课堂提问的误区与对策 8、中学数学教学中的创造性思维的培养 9、浅谈数学教学中的“问题情境” 0、市场经济中的蛛网模型 11、中学数学教学设计前期分析的研究 12、数学课堂差异教学 13、浅谈线性变换的对角化问题 14、圆锥曲线的性质及推广应用 15、经济问题中的概率统计模型及应用 数学教育论文题目(二) 1、二阶变系数齐次微分方程的求解问题 2、一种函数方程的解法 3、微分中值定理的再讨论 4、学生数学学习的障碍研究; 5、中学数学教育中的素质教育的内涵; 6、数学中的美; 7、数学的和谐和统一----谈论数学中的美; 8、推测和猜想在数学中的应用; 9、款买房问题的决策; 10、线性回归在经济中的应用; 11、数学规划在管理中的应用; 12、初等数学解题策略; 13、浅谈数学CAI中的不足与对策; 14、数学创新教育的课堂设计; 15、中学数学教学与学生应用意识培养; 16、关于培养和提高中学生数学学习能力的探究; 17、运用多媒体培养学生 18、高等数学课件的开发 19、 广告 效益预测模型; 数学教育论文题目(三) 1、浅谈菲波纳契数列的内涵和应用价值 2、一道排列组合题的解法探讨及延伸 3、整除与竞赛 4、足彩优化 5、向量的几件法宝在几何中的应用 6、递推关系的应用 7、坐标方法在中学数学中的应用 8、小议问题情境的创设 9、数学概念探索启发式教学 10、柯西不等式的推广与应用 11、关于几个特殊不等式的几种巧妙证法及其推广应用 12、一道高考题的 反思 13、数学中的研究性学习 15、数字危机 16、数学中的化归方法 17、高斯分布的启示 18、 的变形推广及应用 19、网络优化 20、泰勒公式及其应用 猜你喜欢： 1. 数学教育教学论文参考范文 2. 关于数学专业毕业论文题目参考 3. 数学教育专业毕业论文 4. 有关数学教育的论文范文 5. 数学教育专业毕业论文参考
根据我搜集的一些网站来看，建议看看这个，要做毕业论文以及毕业设计的，推荐一个网站 ，里面的毕业设计什么的全是优秀的，因为精挑细选的，网上很少有，都是相当不错的毕业论文和毕业设计，对毕业论文的写作有很大的参考价值，希望对你有所帮助。 别的相关范文很多的，推荐一些比较好的范文写作网站，希望对你有帮助，这些精选的范文网站，里面有大量的范文，也有各种文章写作方法，注意事项，应该有适合你的，自己动手找一下，可不要照搬啊，参考一下，用自己的语言写出来那才是自己的。 如果你不是校园网的话，请在下面的网站找： 毕业论文网： 分类很细 栏目很多 毕业论文： 毕业设计： 开题报告： 实习论文： 写作指导：
要的话请联系我邮箱（点我可见）。13 【篇名】 偏微分方程组的对称群及其在弹性力学方程组中应用 CAJ原文下载 PDF原文下载 【作者】 张鸿庆. 朝鲁. 唐立民. 【刊名】 大连理工大学学报 1997年03期 编辑部Email 《中文核心期刊要目总览》来源期刊 “中国期刊方阵”入选期刊 ASPT来源刊 CJFD收录期刊 【机构】 大连理工大学数学科学研究所. 大连理工大学工程力学研究所. 【关键词】 偏微分方程. 弹性力学. 对称群／不变向量场. 符号运算. 【聚类检索】 同类文献 引用文献 被引用文献 【摘要】 给出了非退化线性偏微分方程组及二次型泛函对称群的不变向量场的一般形式和一类特殊形式非线性偏微分方程组对称群的简化计算条件；利用以上结论及作者以往工作，借助符号运算语言MathematicaTM计算了平面弹性力学方程组一阶Lie－Bactlund对称群的不变向量场，以及应力函数对应的三维弹性力学方程组的Lie代数．为构造弹性力学方程组的一类广泛精确解及守恒律提供了必要的基础，并说明了结论对计算偏微分方程组对称群时的简化作用 【光盘号】 SCTC9706 14 【篇名】 力学中一类变系数微分方程可调参数模型解法 CAJ原文下载 PDF原文下载 【作者】 赵文福. 封营儒. 连星耀. 黎明安. 【刊名】 西安理工大学学报 1995年02期 编辑部Email CJFD收录期刊 【机构】 西安理工大学机械工程系. 【关键词】 可调参数. 变系数微分方程. 非均匀控制参数. 【聚类检索】 同类文献 引用文献 被引用文献 【摘要】 结合一种非均匀控制参数，提出了一种变系数微分方程的可调整参数模型解法，可以很方便地处理由于物理上、几何上的非均匀、非线性而导致数学上的变系数微分方程，应用这种模型可以用非常少的单元得到较满意的数值结果。 【光盘号】 SCTC9508 31 【篇名】 材料力学弯曲问题中集中量与分布量的统一处理 CAJ原文下载 PDF原文下载 【作者】 周锡勤. 张存道. 【刊名】 现代电力 1995年02期 编辑部Email CJFD收录期刊 【机构】 北京动力经济学院. 【关键词】 集中量. 分布量. 弯曲变形. 【聚类检索】 同类文献 引用文献 被引用文献 【摘要】 介绍了利用δ函数统一处理集中量与分布量的一般方法。着重讨论了这种方法在建立含集中量的杆件弯曲时的平衡微分方程的应用，从而推广了材料力学中杆件弯曲时的平衡微分方程。该方程更全面更精确地反映了杆件弯曲这一物理现象。作者把它称为梁弯曲时的广义平衡微分方程。 【光盘号】 SCTC95S5 38 【篇名】 双相材料空间中平片界面裂纹问题的超奇异积分-微分方程 CAJ原文下载 PDF原文下载 【作者】 乐金朝. 汤任基. 【刊名】 科学通报 1996年15期 编辑部Email 《中文核心期刊要目总览》来源期刊 “中国期刊方阵”入选期刊 ASPT来源刊 CJFD收录期刊 【机构】 郑州工学院道路检测与CAE技术研究中心. 上海交通大学工程力学系 郑州 450002 . 上海 200030. 【关键词】 双相材料. 平片界面裂纹. 超奇异积分-微分方程. 【聚类检索】 同类文献 引用文献 被引用文献 【摘要】 <正> 随着复合材料的广泛应用,界面断裂力学成为国际断裂界的前沿研究课题,该领域的研究工作引起了国内外力学家、金属物理学家及材料科学家的广泛关注,并取得了许多新进展。据作者所知,目前的工作主要是研究二维问题,由于数学和力学等方面的困难,三维界面断裂力学方面的研究工作报道较少。本文利用双相材料空间在集中力作用下的弹性力学基本解,使用边界元法,在有限部积分的意义下将任意形状的平片界面裂纹问题归结为一组以裂纹面上的位移间断为未知函数的超奇异积分-微分方程。此组方程对于进一步开展三维界面断裂力学问题的研究具有重要意义。 【光盘号】 SCTA96S4 39 【篇名】 常微分方程的不变式在量子力学中的应用 CAJ原文下载 PDF原文下载 【作者】 杨进. 【刊名】 大学物理 1998年08期 编辑部Email 《中文核心期刊要目总览》来源期刊 CJFD收录期刊 【机构】 成都气象学院基础科学系. 【关键词】 常微分方程. 不变式. 库仑场. 【聚类检索】 同类文献 引用文献 被引用文献 【摘要】 利用常微分方程的不变式，非常方便地求解了一些量子力学问题． 【光盘号】 SCTA9809 40 【篇名】 保守力系的变形拉格朗日方程及其应用 CAJ原文下载 PDF原文下载 【作者】 梁志强. 【刊名】 泰安师专学报 2000年06期 编辑部Email CJFD收录期刊 【机构】 泰安师专物理系!山东泰安271000. 【关键词】 Lagrandge方程. 轨道微分方程. 轨道方程. 【聚类检索】 同类文献 引用文献 被引用文献 【摘要】 从保守力系的拉格朗日方程出发 ,导出一种用于求解保守系统轨道微分方程的变形拉格朗日方程。并将其应用于有心力问题及抛体问题 ,导出了有心力问题的轨道微分方程Binet公式及抛体轨道方程。保守力系的变形拉格朗日方程提供了求解运动物体轨道方程的新方法 ,同时也丰富了分析力学的教学内容。 【光盘号】 SOCI0105
古代方程发展史中国古代是一个在世界上数学领先的国家，用近代科目来分类的话，可以看出无论在算术、代数、几何和三角各方而都十分发达。现在就让我们来简单回顾一下初等数学在中国发展的历史。 (一)属于算术方面的材料 大约在3000年以前中国已经知道自然数的四则运算，这些运算只是一些结果，被保存在古代的文字和典籍中。乘除的运算规则在后来的“孙子算经”(公元三世纪)内有了详细的记载。中国古代是用筹来计数的，在我们古代人民的计数中，己利用了和我们现在相同的位率，用筹记数的方法是以纵的筹表示单位数、百位数、万位数等；用横的筹表示十位数、千位数等，在运算过程中也很明显的表现出来。“孙子算经”用十六字来表明它，“一从十横，百立千僵，千十相望，万百相当。” 和其他古代国家一样，乘法表的产生在中国也很早。乘法表中国古代叫九九，估计在2500年以前中国已有这个表，在那个时候人们便以九九来代表数学。现在我们还能看到汉代遗留下来的木简(公元前一世纪)上面写有九九的乘法口诀。 现有的史料指出，中国古代数学书“九章算术”(约公元一世纪前后)的分数运算法则是世界上最早的文献，“九章算术”的分数四则运算和现在我们所用的几乎完全一样。 古代学习算术也从量的衡量开始认识分数，“孙子算经”(公元三世纪)和“夏候阳算经”(公元六、七世纪)在论分数之前都开始讲度量衡，“夏侯阳算经”卷上在叙述度量衡后又记着：“十乘加一等，百乘加二等，千乘加三等，万乘加四等；十除退一等，百除退二等，千除退三等，万除退四等。”这种以十的方幂来表示位率无疑地也是中国最早发现的。 小数的记法，元朝(公元十三世纪)是用低一格来表示，如作1356 。在算术中还应该提出由公元三世纪“孙子算经”的物不知数题发展到宋朝秦九韶(公元1247年)的大衍求一术，这就是中国剩余定理，相同的方法欧洲在十九世纪才进行研究。 宋朝杨辉所著的书中(公元1274年)有一个1—300以内的因数表，例如297用“三因加一损一”来代表，就是说297=3×11×9，(11＝10十1叫加一，9＝10—1叫损一)。杨辉还用“连身加”这名词来说明201—300以内的质数。 (二)属于代数方面的材料 从“九章算术”卷八说明方程以后，在数值代数的领域内中国一直保持了光辉的成就。 “九章算术”方程章首先解释正负术是确切不移的，正象我们现在学习初等代数时从正负数的四则运算学起一样，负数的出现便丰富了数的内容。 我们古代的方程在公元前一世纪的时候已有多元方程组、一元二次方程及不定方程几种。一元二次方程是借用几何图形而得到证明。 不定方程的出现在二千多年前的中国是一个值得重视的课题，这比我们现在所熟知的希腊丢番图方程要早三百多年。具有x3+px2+qx=A和x3+px2=A形式的三次方程，中国在公元七世纪的唐代王孝通“缉古算经”已有记载，用“从开立方除之”而求出数字解答(可惜原解法失传了)，不难想象王孝通得到这种解法时的愉快程度，他说谁能改动他著作内的一个字可酬以千金。 十一世纪的贾宪已发明了和霍纳(1786—1837)方法相同的数字方程解法，我们也不能忘记十三世纪中国数学家秦九韶在这方面的伟大贡献。 在世界数学史上对方程的原始记载有着不同的形式，但比较起来不得不推中国天元术的简洁明了。四元术是天元术发展的必然产物。 级数是古老的东西，二千多年前的“周髀算经”和“九章算术”都谈到算术级数和几何级数。十四世纪初中国元代朱世杰的级数计算应给予很高的评价，他的有些工作欧洲在十八、九世纪的著作内才有记录。十一世纪时代，中国已有完备的二项式系数表，并且还有这表的编制方法。 历史文献揭示出在计算中有名的盈不足术是由中国传往欧洲的。 内插法的计算，中国可上溯到六世纪的刘焯，并且七世纪末的僧一行有不等间距的内插法计算。 十四世纪以前，属于代数方面许多问题的研究，中国是先进国家之一。 就是到十八，九世纪由李锐(1773—1817)，汪莱(1768—1813)到李善兰(1811—1882)，他们在这一方面的研究上也都发表了很多的名著。
人类对一元二次方程的研究经历了漫长的岁月，早在公元前2000年左右，居住在底格里斯河和幼法拉底河的古巴比伦人已经能解一些一元二次方程。而在中国，《九章算术》“勾股”章中就有一题：“今有户高多于广六尺八寸，两隅相去适一丈，问户高、广各几何？。”之后的丢番图（古代希腊数学家），欧几里德（古代希腊数学家），赵爽，张遂，杨辉对一元二次方程的贡献更大贝祖（Bezout Etienne ）法国数学家。少年时酷爱数学，主要从事方程论研究。他是最先认识到行列式价值的数学家之一。最早证明了齐次线性方程组有非零解的条件是系数行列式等于零。他在其第一篇论文《几种类型的方程》中用消元法将只含一个未知数的n次方程问题与解联立方程组问题联系起来，提供了某些n次方程的解法。他还用消元法解次数高于1的两个二元方程，并证明了关于方程次数的贝祖定理。1086～1093年，中国宋朝的沈括在《梦溪笔谈》中提出“隙积术”和“会圆术”，开始高阶等差级数的研究。 十一世纪，阿拉伯的阿尔·卡尔希第一次解出了二次方程的根。 十一世纪，阿拉伯的卡牙姆完成了一部系统研究三次方程的书《代数学》。 十一世纪，埃及的阿尔·海赛姆解决了“海赛姆”问题，即要在圆的平面上两点作两条线相交于圆周上一点，并与在该点的法线成等角。 十一世纪中叶，中国宋朝的贾宪在《黄帝九章算术细草》中，创造了开任意高次幂的“增乘开方法”，并列出了二项式定理系数表，这是现代“组合数学”的早期发现。后人所称的“杨辉三角”即指此法。 十二世纪，印度的拜斯迦罗著《立刺瓦提》一书，这是东方算术和计算方面的重要著作。 1202年，意大利的裴波那契发表《计算之书》，把印度—阿拉伯记数法介绍到西方。 1220年，意大利的裴波那契发表《几何学实习》一书，介绍了许多阿拉伯资料中没有的示例。 1247年，中国宋朝的秦九韶著《数书九章》共十八卷，推广了“增乘开方法”。书中提出的联立一次同余式的解法，比西方早五百七十余年。 1248年，中国宋朝的李治著《测圆海镜》十二卷，这是第一部系统论述“天元术”的著作。 1261年，中国宋朝的杨辉著《详解九章算法》，用“垛积术”求出几类高阶等差级数之和。 1274年，中国宋朝的杨辉发表《乘除通变本末》，叙述“九归”捷法，介绍了筹算乘除的各种运算法。 1280年，元朝《授时历》用招差法编制日月的方位表(中国 王恂、郭守敬等)。 十四世纪中叶前，中国开始应用珠算盘。 1303年，中国元朝的朱世杰著《四元玉鉴》三卷，把“天元术”推广为“四元术”。 1464年，德国的约·米勒在《论各种三角形》(1533年出版)中，系统地总结了三角学。 1494年，意大利的帕奇欧里发表《算术集成》，反映了当时所知道的关于算术、代数和三角学的知识。 1545年，意大利的卡尔达诺、费尔诺在《大法》中发表了求三次方程一般代数解的公式。 1550～1572年，意大利的邦别利出版《代数学》，其中引入了虚数，完全解决了三次方程的代数解问题。 1591年左右，德国的韦达在《美妙的代数》中首次使用字母表示数字系数的一般符号，推进了代数问题的一般讨论。 1596～1613年，德国的奥脱、皮提斯库斯完成了六个三角函数的每间隔10秒的十五位小数表。 1614年，英国的耐普尔制定了对数。 1615年，德国的开卜勒发表《酒桶的立体几何学》，研究了圆锥曲线旋转体的体积。 1635年，意大利的卡瓦列利发表《不可分连续量的几何学》，书中避免无穷小量，用不可分量制定了一种简单形式的微积分。 1637年，法国的笛卡尔出版《几何学》，提出了解析几何，把变量引进数学，成为“数学中的转折点”。 1638年，法国的费尔玛开始用微分法求极大、极小问题。 1638年，意大利的伽里略发表《关于两种新科学的数学证明的论说》，研究距离、速度和加速度之间的关系，提出了无穷集合的概念，这本书被认为是伽里略重要的科学成就。 1639年，法国的迪沙格发表了《企图研究圆锥和平面的相交所发生的事的草案》，这是近世射影几何学的早期工作。 1641年，法国的帕斯卡发现关于圆锥内接六边形的“帕斯卡定理”。 1649年，法国的帕斯卡制成帕斯卡计算器，它是近代计算机的先驱。 1654年，法国的帕斯卡、费尔玛研究了概率论的基础。 1655年，英国的瓦里斯出版《无穷算术》一书，第一次把代数学扩展到分析学。 1657年，荷兰的惠更斯发表了关于概率论的早期论文《论机会游戏的演算》。 1658年，法国的帕斯卡出版《摆线通论》，对“摆线”进行了充分的研究。 1665～1676年，牛顿(1665～1666年)先于莱布尼茨(1673～1676年)制定了微积分，莱布尼茨(1684～1686年)早于牛顿(1704～1736年)发表了微积分。 1669年，英国的牛顿、雷夫逊发明解非线性方程的牛顿—雷夫逊方法。 1670年，法国的费尔玛提出“费尔玛大定理”。 1673年，荷兰的惠更斯发表了《摆动的时钟》，其中研究了平面曲线的渐屈线和渐伸线。 1684年，德国的莱布尼茨发表了关于微分法的著作《关于极大极小以及切线的新方法》。 1686年，德国的莱布尼茨发表了关于积分法的著作。 1691年，瑞士的约·贝努利出版《微分学初步》，这促进了微积分在物理学和力学上的应用及研究。 1696年，法国的洛比达发明求不定式极限的“洛比达法则”。 1697年，瑞士的约·贝努利解决了一些变分问题，发现最速下降线和测地线。 1704年，英国的牛顿发表《三次曲线枚举》《利用无穷级数求曲线的面积和长度》《流数法》。 1711年，英国的牛顿发表《使用级数、流数等等的分析》。 1713年，瑞士的雅·贝努利出版了概率论的第一本著作《猜度术》。 1715年，英国的布·泰勒发表《增量方法及其他》。 1731年，法国的克雷洛出版《关于双重曲率的曲线的研究》，这是研究空间解析几何和微分几何的最初尝试。 1733年，英国的德·勒哈佛尔发现正态概率曲线。 1734年，英国的贝克莱发表《分析学者》，副标题是《致不信神的数学家》，攻击牛顿的《流数法》，引起所谓第二次数学危机。 1736年，英国的牛顿发表《流数法和无穷级数》。 1736年，瑞士的欧拉出版《力学、或解析地叙述运动的理论》，这是用分析方法发展牛顿的质点动力学的第一本著作。 1742年，英国的麦克劳林引进了函数的幂级数展开法。 1744年，瑞士的欧拉导出了变分法的欧拉方程，发现某些极小曲面。 1747年，法国的达朗贝尔等由弦振动的研究而开创偏微分方程论。 1748年，瑞士的欧拉出版了系统研究分析数学的《无穷分析概要》，这是欧拉的主要著作之一。 1755～1774年，瑞士的欧拉出版了《微分学》和《积分学》三卷。书中包括微分方程论和一些特殊的函数。 1760～1761年，法国的拉格朗日系统地研究了变分法及其在力学上的应用。 1767年，法国的拉格朗日发现分离代数方程实根的方法和求其近似值的方法。 1770～1771年，法国的拉格朗日把置换群用于代数方程式求解，这是群论的开始。 1772年，法国的拉格朗日给出三体问题最初的特解。 1788年，法国的拉格朗日出版了《解析力学》，把新发展的解析法应用于质点、刚体力学。 1794年，法国的勒让德出版流传很广的初等几何学课本《几何学概要》。 1794年，德国的高斯从研究测量误差，提出最小二乘法，于1809年发表。 1797年，法国的拉格朗日发表《解析函数论》，不用极限的概念而用代数方法建立微分学。 1799年，法国的蒙日创立画法几何学，在工程技术中应用颇多。 1799年，德国的高斯证明了代数学的一个基本定理：实系数代数方程必有根。 微分方程:大致与微积分同时产生 。事实上，求y′＝f(x)的原函数问题便是最简单的微分方程。I.牛顿本人已经解决了二体问题：在太阳引力作用下，一个单一的行星的运动。他把两个物体都理想化为质点，得到3个未知函数的3个二阶方程组，经简单计算证明，可化为平面问题，即两个未知函数的两个二阶微分方程组。用现在叫做“首次积分”的办法，完全解决了它的求解问题。17世纪就提出了弹性问题，这类问题导致悬链线方程、振动弦的方程等等。总之，力学、天文学、几何学等领域的许多问题都导致微分方程。在当代，甚至许多社会科学的问题亦导致微分方程，如人口发展模型、交通流模型……。因而微分方程的研究是与人类社会密切相关的。当初，数学家们把精力集中放在求微分方程的通解上，后来证明这一般不可能，于是逐步放弃了这一奢望，而转向定解问题：初值问题、边值问题、混合问题等。但是，即便是一阶常微分方程，初等解(化为积分形式)也被证明不可能，于是转向定量方法（数值计算）、定性方法，而这首先要解决解的存在性、唯一性等理论上的问题。 方程对于学过中学数学的人来说是比较熟悉的；在初等数学中就有各种各样的方程，比如线性方程、二次方程、高次方程、指数方程、对数方程、三角方程和方程组等等。这些方程都是要把研究的问题中的已知数和未知数之间的关系找出来，列出包含一个未知数或几个未知数的一个或者多个方程式，然后取求方程的解。 但是在实际工作中，常常出现一些特点和以上方程完全不同的问题。比如：物质在一定条件下的运动变化，要寻求它的运动、变化的规律；某个物体在重力作用下自由下落，要寻求下落距离随时间变化的规律；火箭在发动机推动下在空间飞行，要寻求它飞行的轨道，等等。 物质运动和它的变化规律在数学上是用函数关系来描述的，因此，这类问题就是要去寻求满足某些条件的一个或者几个未知函数。也就是说，凡是这类问题都不是简单地去求一个或者几个固定不变的数值，而是要求一个或者几个未知的函数。 解这类问题的基本思想和初等数学解方程的基本思想很相似，也是要把研究的问题中已知函数和未知函数之间的关系找出来，从列出的包含未知函数的一个或几个方程中去求得未知函数的表达式。但是无论在方程的形式、求解的具体方法、求出解的性质等方面，都和初等数学中的解方程有许多不同的地方。 在数学上，解这类方程，要用到微分和导数的知识。因此，凡是表示未知函数的导数以及自变量之间的关系的方程，就叫做微分方程。 微分方程差不多是和微积分同时先后产生的，苏格兰数学家耐普尔创立对数的时候，就讨论过微分方程的近似解。牛顿在建立微积分的同时，对简单的微分方程用级数来求解。后来瑞士数学家雅各布�6�1贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论。 常微分方程的形成与发展是和力学、天文学、物理学，以及其他科学技术的发展密切相关的。数学的其他分支的新发展，如复变函数、李群、组合拓扑学等，都对常微分方程的发展产生了深刻的影响，当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具。 牛顿研究天体力学和机械力学的时候，利用了微分方程这个工具，从理论上得到了行星运动规律。后来，法国天文学家勒维烈和英国天文学家亚当斯使用微分方程各自计算出那时尚未发现的海王星的位置。这些都使数学家更加深信微分方程在认识自然、改造自然方面的巨大力量。 微分方程的理论逐步完善的时候，利用它就可以精确地表述事物变化所遵循的基本规律，只要列出相应的微分方程，有了解方程的方法。微分方程也就成了最有生命力的数学分支。
■ 有些微分方程求不出函数解(解析解)，只能求数值解，MMA软件的函数命令 tt＝NDSolve[微分方程]，然后 ▲赋值ⅹ＝2，求出 y＝？ ▲赋值 x＝3，求出 y＝？ ··· 赋值ⅹ＝n，求出 y＝？，这些就是微分方程的数值解。虽然解不出未知函数y(ⅹ)表达式，但MMA可画出它的函数图像，很复杂的图像都能画出来。也碰到过特例，从(ⅹ0)向左图像就没了，对y(x)赋值后发现，x≤xo时，函数值y(ⅹ)变成复数了，包括( 1、ⅰ )二个维度，MMA当然无法画图了。多数工程技术出现的微分方程组，总求不出函数解析式，所以数值解的意义和作用不言而喻。■ 从数值分析来看，偏微分方程及微分方程数值解常用二种方法。① 差分法~原理是用《差商》替代微商(导数)。②有限元法~原理是泛函变分法。将微分方程边值问题→泛函求极值问题→线性代数方程求解。MMA求解数值解时在各种方法中选择最优法。