热 拉普拉斯方程的n维求解器的实现如下: 在空间上使用二阶方案,并在时间上使用后向欧拉方案。 为了解决所得的线性系统,使用了共轭梯度法。 例子: y(t = 0)= 10sin(x)的初始条件,扩散系数为0.1