[求助]差分法求解微分方程(万分感谢)
 

对下列定解问题
du/dt=(4/(pi)^2)d(du)/dx^2,当0<x<4,t>0时；
u(x,0)=sin(pi/4*x)(1+2cos(pi/4*x))，当0<=x<=4时；
u(0,t)=u(4,t)=0,当t>0时；
取h=0.2,r=1,分别用古典显格式，古典隐格式和Crank-Nicolson格式计算t=0.4时的近似解，并与精确解u(x,t)=
e^(-t)sin(pi/2*x)+e^(-t/4)sin(pi/4*x)比较。
小弟初学MATLAB,老师布置了上题让我们写上面题目的程序,请热心人帮忙解答,谢谢大家了