sovamatlab测试代码

Ⅰ Matlab 这段代码怎么成功运行

functionhtt

m=0;
x=linspace(0,10,20)*1e-2;
t=linspace(0,15,16);

sol=pdepe(m,@pdefun,@icbun,@bcfun,x,t);

figure
mesh(x,t,sol(:,:,1))%温度与时间和空间位置的关系图

figure
%画1、2、4、6、8、15s时刻温度分布图
%1s时刻，（因为本题sol第一行为0时刻）

plot(x,sol(2,:,1))

holdon

plot(x,sol(3,:,1))
plot(x,sol(5,:,1))
plot(x,sol(7,:,1))
plot(x,sol(9,:,1))
plot(x,sol(16,:,1))

%pdefun.m%偏微分方程(一维动态传热)
function[c,f,s]=pdefun(x,t,u,dx)
c=1/2e-4;f=dx;s=0;

%icbun.m%偏微分方程初始条件(一维动态传热)
functionu0=icbun(x)
u0=20;

%bcfun.m%偏微分方程边界条件(一维动态传热)
function[pl,ql,pr,qr]=bcfun(xl,ul,xr,ur,t)
pl=ul-120;ql=0;pr=ur-20;qr=0;