1) 令n?0。 2)
k?1。
3) 求解三对角线性方程组(3.3a)得到差分解
?un?1/21,kn?1/2n?1/2,u2,,?,uNk?1,k?。
4) 若k?N,则k增加1,转步骤3)。否则转5)。 5) 令j?1。
6) 求解三对角线性方程组(3.3b)得到差分解
?unj,1n,unj,2,?,uj,N?1?。
7) 若j?N,则j增加1,转步骤6)。否则转8)。 8) 若n?T/?,则n增加1,转步骤2)。否则结束。
LOD格式的基本想法是,由第n层计算n?1/2层时,对
k?1,?,N?1,依次固定y?yk,然后计算y?yk这条直线上各个
网点上的近似值;因为这时y不变,所以原来的二维微分方程退化为关于x的一维微分方程。接着,当由第n层计算n?1/2层时,则依次固定x?xj。LOD格式可以直接推广到任意维抛物方程。LOD格式对任意网比都是稳定的,截断误差阶和收敛阶是O(?2?h2)。