1、古典显式格式求解抛物型偏微分方程（一维热传导方程） 2、古典隐式格式求解抛物型偏微分方程（一维热传导方程） 3、Crank-Nicolson隐式格式求解抛物型偏微分方程 4、正方形区域Laplace方程Diriclet问题的求解 如： function [U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C) %古典显式格式求解抛物型偏微分方程 %[U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C) % %方程：u_t=C*u_xx 0 <= x <= uX,0 <= t <= uT %初值条件：u(x,0)=phi(x) %边值条件：u(0,t)=psi1(t), u(uX,t)=psi2(t)

学数值分析时，用c++编写的一个用于计算三体问题的程序。 可能有人用得着。

MATLAB 提供了两种方法解决PDE 问题，一是pdepe()函数，它可以求解一般的PDEs，据用较大的通用性，但 只支持命令行形式调用。二是PDE 工具箱，可以求解特殊PDE 问题，PDEtool 有较大的局限性，比如只能求解 二阶PDE 问题，并且不能解决偏微分方程组，但是它提供了GUI 界面，从繁杂的编程中解脱出来了，同时 可 以通过File->Save As 直接生成M 代码

MATLAB

§1、引言 §2、初值问题的数值解法--单步法 §3、龙格-库塔方法 §4、收敛性与稳定性 §5、初值问题的数值解法―多步法 §6、方程组和刚性方程 §7、习题和总结

Introduction to Finite and Spectral Element Methods Using MATLAB®Second Edition This book contains information obtained from authentic and highly regarded sources. Reasonable efforts have been made to publish reliable data and information, but the author and publisher cannot assume responsibility for the validity of all materials or the consequences of their use. The authors and publishers have attempted to trace the copyright holders of all material reproduced in this publication and apologize to copyright holders if permission to publish in this form has not been obtained. If any copyright material has not been acknowledged please write and let us know so we may rectify in any future reprint.

Adams隐式公式 注：在这些方法中，4阶的Adams预测校正方法具有方法 简洁、使用方便、易排序、高精度等优点。尤其当函数 f比较复杂时更能显出它的优越性。 k p 公式 1 2 2 3 3 4 4 5

第二章线性方程组的数值解法线性方程组对于线性方程组Ax=b,其中若系数阵A非奇异，则方程组有唯一解.本章只讨论系数矩阵为非奇异的线性方程组.计算机上求解线性方程

1、古典显式格式求解抛物型偏微分方程（一维热传导方程） 2、古典隐式格式求解抛物型偏微分方程（一维热传导方程） 3、Crank-Nicolson隐式格式求解抛物型偏微分方程 4、正方形区域Laplace方程Diriclet问题的求解 如： function [U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C) %古典显式格式求解抛物型偏微分方程 %[U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C) % %方程：u_t=C*u_xx 0 <= x <= uX,0 <= t <= uT %初值条件：u(x,0)=phi(x) %边值条件：u(0,t)=psi1(t), u(uX,t)=psi2(t)

介绍了应用最为广泛的椭圆型、双曲型、抛物型偏微分方程的数值解法，而且还详细编程实现了每种方程的多种常见数值解法。 附件使用MATLAB编程来实现这些算法。