论文部分内容阅读
对扩散方程提出了精度为O(τ3+h2)的差分格式,首先对空间变量中心差分格式离散,所得到常微分方程组利用指数函数的Pade[ 2/1]逼近,得到空间二阶时间三阶精度的两层绝对稳定的隐式差分格式,并讨论了稳定性.数值结果与Crank-Nicholson格式进行比较,数值结果表明,该格式不但有效地解决初始边界条件间断问题,而且适合于大时间步长问题.