matlab pdebc,pde求解的MATLAB语言求助

如下网址上的MATLAB语言，

http://www.matlabsky.com/thread-981-1-1.html

我想按照他的代码写在M文件中，运行一下，但是一直运行不出来。

第四行，他就说我没有定义

??? Input argument "ua" is undefined.

Error in ==> pde at 4

pa=[0;ua(2)];

以下是他的代码。求助大神帮我改一下。谢谢

%% 目标PDE函数

%%syms x t u

function [c,f,s]=pdefun (x,t,u,du);

c=[1;1];

%%dudx=diff(u,'x');

f=[0.024*du(1);0.17*du(2)];

temp=u(1)-u(2);

s=[-1;1].*(exp(5.73*temp)-exp(-11.46*temp));

%% 边界条件函数

function [pa,qa,pb,qb]=pdebc(xa,ua,xb,ub,t)

%a表示下边界，b表示上边界

pa=[0;ua(2)];

qa=[1;0];

pb=[ub(1)-1;0];

qb=[0;1];

%% 初值条件函数

function u0=pdeic(x)

u0=[1;0];

%% 主调函数

clc

x=0:0.05:1;

t=0:0.05:2;

m=0;

sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);

figure('numbertitle','off','name','PDE Demo——by Matlabsky')

subplot(211)

surf(x,t,sol(:,:,1))

title('The Solution of u_1')

xlabel('X')

ylabel('T')

zlabel('U')

subplot(212)

surf(x,t,sol(:,:,2))

title('The Solution of u_2')

xlabel('X')

ylabel('T')

zlabel('U')