麻烦大家帮我看看这个程序该怎么调用
大家好,我编了一个M文件解偏微分方程组,但不知道怎么调用,M文件如下:function pdex4
m=0;
x=;
t=;
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=pdex4pde(x,t,u,DuDx)
c=; f=[-0.5668*u1;0];F=5.0787*u2-u1;s=;
%.......................................................
function u0=pdex4ic(x);
u0=;
%.......................................................
function=pdex4bc(xl,cl,t)
pl=;ql=;
我最后想求的是x=1.0时各个t所对应的u1,u2值。本人对MATLAB不太懂,哪位高手能不能指点我一下?太感谢了
[ 本帖最后由 eight 于 2007-9-22 22:50 编辑 ]
回复 #1 星月女孩 的帖子
代码放到一个m文件里,存成pdex4.m.如果程序没有错误,那么在命令行输入pdex4应该就可以看到结果了。
谢谢!我试试看
页:
[1]