matlab du/dt=d(du)/dx^2 x属于（0,1）,t属于(0,T]u(0,t)=u(1,t)=0u(x,

matlab du/dt=d(du)/dx^2 x属于（0,1）,t属于(0,T]u(0,t)=u(1,t)=0u(x,0)=0若在u(1/2,0)处有一个扰动1/(2^10),取h=1/16,r=1/2分别用古典显格式和richardson格式计算8层1)打印出第8层上各结点处的计算值2）预测继续算下去计算值的变化趋势3）分析上述趋势产生的原因要用matlab编程望告知具体程序~55555自己做不来啊………= =|||都没有人会做？分数已经到上限了