i don't find where is the error plz help me

1 次查看（过去 30 天）

显示更早的评论

mouna 2011-5-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/6607-i-don-t-find-where-is-the-error-plz-help-me

关闭： Walter Roberson 2016-1-22

在 MATLAB Online 中打开

hi,

this is my program but i don't find where is the error plz help me it's urgent. thanks for your help in advance.

function pdex
m = 0;
z = linspace(0,100,100);
t = linspace(0,5,5);
sol = pdepe(m,@pdexpde,@pdexic,@pdexbc,z,t);
C = sol(:,:,1);
figure;
plot(z,C(end,:));
% ----------------------------------------------------------
function [g,f,s] = pdexpde(z,t,C,DCDz);
rhos=0.55;
rhof=0.385;
epsilon=0.64;
u=0.05;
keq=23.53;
A=1+(((1-epsilon)/epsilon)*keq*(rhos/rhof));
B=u/epsilon;
g=A;
f=0;
s=(-B)*DCDz;
% ----------------------------------------------------------
function C0 = pdexic(z);
a=0.0015;
C0=a;
% ----------------------------------------------------------
function [pl,ql,pr,qr] = pdexbc(zl,Cl,zr,Cr,t)
pl = 0;                        
ql = 0;                             
pr = 0;                           
qr = 1;