使用SciPy绘制含时扩散方程解析解的多时刻图像

我想要实现的功能是：绘制含时扩散方程（Time Dependent Diffusion Equation）解析解在不同时刻的图像。具体逻辑很明确——固定某一个时刻的t值，只改变自变量x，然后绘制对应的因变量y的曲线。

目前我已经写了部分代码，如下：

import numpy as np
from scipy.sparse import diags
import scipy as sp
import scipy.sparse
from scipy import special
import matplotlib.pyplot as plt

def Canalytical(intervals, D=1):
    points = 1000
    x=np.linspace(0, 1, intervals+1)
    t=1
    c=np.ones([intervals+1])
    sm = 0
    pos = 0
    for xi in x...

（注：代码里的for xi in x...还未完成，后续可以根据具体的解析解形式补充计算逻辑——比如如果是有限边界的级数解，就可以在循环里累加级数项；如果是无限域的高斯型解，直接调用scipy.special里的误差函数就能快速完成计算）

内容的提问来源于stack exchange，提问作者toumperlekis