吧里关心微分方程的人不多，但这个例子确实很独特，姑且提一句。我们知道（好吧，相当多的人不知道，想着自己写经典龙格库塔的人时不时就会冒出来），NDSolve默认的常微分方程求解器（细节不明——甚至可以说是NDSolve中最不透明的部分——似乎一般是LSODA方法，有时似乎也会调用StiffnessSwitching方法）性能是非常好的，在NDSolve求解出错的时候，常微分方程求解器永远应该放在最后检查。而当常微分方程求解器出现性能问题时，“使用MaxStepSize限制步长”又可以解决掉绝大多数问题，求解方法导致的问题可以说是少之又少，而上个月，我居然发现了一个。完整讨论请参见SE帖子
《Boundary condition for heat equation in polar coordinates deduced with L'Hôpital's rule fails for method of lines, but works well for FDM》（编号305852）
简单地说，我发现了一类似乎只能用特定的隐式方法求解的常微分方程组，而默认的常微分方程求解器
并不总是能正确地处理它。