$\displaystyle \frac{\partial u}{\partial t}=\alpha(x,t)\cdot \frac{\partial^2 u}{\partial x^2}+b(x,t)$

$u(x,0)=f(x)$初始条件

$u_x(0,t)=0$第二类边界条件

$u_x(1,t)=0$第二类边界条件

Maple 中有我的第一种类型 BC 的代码我真的不知道，我应该在这段代码中更改什么来解决我的第二种类型 BC 的问题我应该在循环#Boundary cond中更改什么？我必须在哪里使用中心有限差分作为边界条件？显式方案中会发生什么变化（对于 2nd type BC）？我的第一类 BC 的一维热方程的显式方法代码如下：

如果有人帮助我，我将不胜感激。