热传导偏微分方程差分,使用Crank-Nicloson格式,附MATLAB实例