麻烦大家帮我看看这个程序该怎么调用

星月女孩

大家好，我编了一个M文件解偏微分方程组，但不知道怎么调用，M文件如下：
function pdex4
m=0;
x=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0];
t=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0];
sol=pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t);
u1=sol(:,:,1); u2=sol(:,:,2);
figuge;   surf(x,t,u1);   title('u1(x,t)');   xlabel('x');    ylabel('t');
figuge;   surf(x,t,u2);   title('u2(x,t)');   xlabel('x');    ylabel('t');
%.......................................................
function[c,f,s]=pdex4pde(x,t,u,DuDx)
c=[1;1]; f=[-0.5668*u1;0];F=5.0787*u2-u1;s=[11.6761*F;-4.3774*F];
%.......................................................
function u0=pdex4ic(x);
u0=[0;9.89128];
%.......................................................
function[pl,ql]=pdex4bc(xl,cl,t)
pl=[u1;0];ql=[0;1];
我最后想求的是x=1.0时各个t所对应的u1,u2值。本人对MATLAB不太懂，哪位高手能不能指点我一下？太感谢了

[ 本帖最后由 eight 于 2007-9-22 22:50 编辑 ]

花如月

代码放到一个m文件里，存成pdex4.m.如果程序没有错误，那么在命令行输入pdex4
应该就可以看到结果了。

星月女孩