一维热传导方程的Legendre-Galerkin 方法
~~~~~~~~~ 对下面方程空间方向Legendre-Galerkin method 离散,时间方向采用Crank-Nicolson格式离散。
{ u t = u x x + f x ∈ ( 0 , π ) , t > 0 u ( x , 0 ) = u 0 ( x ) 0 ≤ x ≤ π u ( 0 , t ) = 0 , u ( π , t ) = 0 t > 0 \begin{cases} &u_t=u_{xx}+f &x\in(0,\pi),t>0\\ &u(x,0)=u_0(x) &0\leq x\leq \pi \\ &u(0,t)=0,u(\pi,t)=0 & t>0 \end{cases} ⎩⎪⎨⎪⎧ut=uxx