您所在位置：主题：[求助]差分法求解微分方程
king01965
[专家分：0] 发布于 2009-06-13 09:12:00
对下列定解问题
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)比较。
麻烦各位高手帮忙弄下程序万分感谢!
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)比较。
麻烦各位高手帮忙弄下程序万分感谢!