x=0:0.05:3.0;
t=0:0.05:2;
m=0;
sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);
u1=sol(:,:,1);
function u0=pdeic(x)
u0=1
end
u2=sol(:,:,2);
figure(1);
surf(x,t,u1)
title('u1(x,t)')
xlabel(' 位置 x')
ylabel(' 时间 t')
zlabel(' 电子密度(10^1^6 cm^-3)')
figure(2);
surf(x,t,u2)
title('u2(x,t)')
xlabel('Distance x')
ylabel('Time t')
figure(3);
surf(x,t,u2./u1)
title('u(x,t)')
xlabel('Distance x')
figure(4);
plot(x,u1)