ADI隐式交替法三种解法及误差分析（一般的教材上只说第一种）理论部分参看孙志忠：偏微分方程数值解法注意：1. 最好不要直接看程序，中间很多公式很烦人的（一定要小心），我写了两天，终于写对了。2. 中间：例如r*(u(i-1,m1,k)+u(i+1,m1,k)形式写成分形式：r*u(i-1,m1,k)+r*u(i+1,m1,k)后面会出错，我也不是很清楚为什么，可能由于舍入误差，或者大数吃掉小数的影响。3. 下面有三个程序4. 具体理论看书，先仔细看书（孙志忠：偏微分方程数值解法）或者网上搜一些理论。Matlab程序：1.function u u0 p e x y t=ADI1(h1,h2,m1,m2,n)%ADI解二维抛物线型偏微分方程(P-R交替隐式，截断)%此程序用的是追赶法解线性方程组%h1为空间步长，h2为时间步长%m1,m2分别为x方向，y方向网格数，n为时间网格数%p为精确解，u为数值解,e为误差%定义u0(i,j,k)=u(i,j,k+1/2),因为矩阵中，i,j,k必须全为整数x=(0
