MATLAB：一维非齐次抛物型抛物型方程的向前Euler有限差分数值求解

计算算例为：

\begin{cases}\frac{\partial u}{\partial t}-\frac{\partial^2u}{\partial x^2}=x\mathrm e^t-6x,\quad0<x<1,\quad0<t\leq1,\\\\u(x,0)=x^3+x,\quad0\leqslant x\leqslant1,\\\\u(0,t)=0,\quad u(1,t)=1+\mathrm e^t,\quad0<t\leqslant1.\end{cases}

解析解：u(x,t)=x(x^2+e^t).

采用向前Euler数值计算上述偏微分方程，数值解和解析解对比结果如下：

误差曲线如下：