求解非线性扩散方程，该方程可以线性化，如 Richtmyer & Morton [1] 中的一般非线性扩散方程所示。 该方法是将 pde 线性化并应用 Crank-Nicolson 隐式有限差分格式以数值方式求解方程。 Matlab 运行命令--------------------------
类型： IsoFreeSurfaceSolver
解决 pde： -------------------
\frac{\partial h}{\partial t}=\frac{1}{12}\frac{\partial^2 h^4}{\partial x^2}
pde 可应用于在重力作用下在水平基底上扩散的等温粘性流体流动 - Huppert [2]。 请注意，PDE 已无量纲化。
初始条件： t=0: h = (1 - x^2)_{+} + 10^-6（有预湿膜） 考虑到 x = 0
