高手帮忙看一下

sbstoneface

clear
D0=2.0;
Q=84100;
T=573.15;
R=8.314;
D=D0*exp(-Q/(R*T));%扩散系数
space=10;
h=0.1;%空间步长
time=1200000;
t=100;%时间步长
n=(space/h)+1;
m=(time/t)+1;
c=zeros(n,m);
c(1,1:m)=1.4;%边界条件
c(n,1:m)=0;
c(2:n-1,1)=0;%初始条件
r=D*t/h^2;%傅立叶数
for j=2:m-1
   for i=2:n-1
   c(i,j+1)=c(i,j-1)+2*r*(c(i-1,j)-2*c(i,j)+c(i+1,j)); %中心差商法
end
   end
   w=linspace(0,time,m);
w1=c(51,:);
plot(w,w1,'r')
这是我编的用差分法求解渗碳问题,可画出的图象不符实际,麻烦高手指点一下(已知扩散方程为
dc/dt=Dd^2c/dx^2,c为浓度,t为扩散时间,x为距离)