[b]我最近在解一个带源项的扩散方程，源项是非线性的是一个2维问题，我采用ADI交替方向方法计算，每个方向上采用追赶法求解得到的3对角矩阵，但计算结果显示极不稳定 一会儿正，一会儿负，最后全为NaN了。希望哪位大侠给点建议，有兴趣也可以共同探讨。我的QQ:361228548 希望得到大家的帮助。 谢谢！！！！[/b][em2][em2][em2][em2]