Info
此问题已关闭。 请重新打开它进行编辑或回答。
i don't find where is the error plz help me
1 次查看(过去 30 天)
显示 更早的评论
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;
2 个评论
Walter Roberson
2011-5-2
Do we get a hint about what error you observe? An error message and an associated line?
Jan
2011-5-3
If your question is any kind or "urgent", it would be a very good idea to explain anything about the problem.
回答(0 个)
此问题已关闭。
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!