#stability

隐式扩散与托马斯算法 — 用一次三对角求解换来无条件稳定

r=0.6时显式扩散60步内就冲到1e12。隐式法只做一次托马斯三对角扫描，同一时间步仍保持有限，即便把r推到2也不会炸掉。