[求助]关于pdepe使用的问题

weiz5022

function dd
clear reset;
global De b Cli v0 r1 r0 L0 k Imax km Cmin tt pi;
De=2.0e-06;b=1;Cli=6e+00;V0=1e-07;r1=3.3e-01;r0=2e-02;
L0=2.5e+02;k=8e-03;Imax=1e-05;km=2.0e-02;Cmin=2.0e-03;tt=8.64e+05;
rn1=100;tn2=100;
m=1;r=linspace(r0,r1,rn1);t=linspace(0,tt,tn2);
sol=pdepe(m,@pdefun,@pdeSNic,@pdeSNbc,r,t);
u=sol(:,:,1);
figure;plot(r,u(end,:));
xlabel('\fontsize{12}the distance from the center');
ylabel('\fontsize{12}盐浓度');
nt=length(t);iok=2:nt;
for j=iok;[uout(iok),I(j)]=pdeval(m,r,u(j,:),r0);end;I;
figure;C=u(:,2);
f3=Imax*(C-Cmin)./(km+C-Cmin);
plot(t,f3);
xlabel('\fontsize{12}时间');ylabel('\fontsize{12}表面盐通量');
f4=2*pi*r0*L0*f3*tt/nt;
T1=sum(f4);
LL(1)=0;
for h=1:tn2-1;LL(h+1)=LL(h)+k*tt/nt;LL(h+1);end;LL;
LN=2*pi*r0*LL'.*f3*tt/nt;
T2=sum(LN)
T=T1+T2
function [c,f,s]=pdefun(r,t,u,DuDr)
global De b r0 v0;
c=1;f=De*DuDr+r0*v0*u/(b*r);s=0;
function u0=pdeSNic(r)
global Cli;
u0=Cli;
function [pl,ql,pr,qr]=pdeSNbc(rl,ul,rr,ur,t)
global De b Cli v0 Imax km Cmin;
pl=v0*ul-Imax*(ul-Cmin)./(km+ul-Cmin);ql=De*b;
pr=ur-Cli;qr=0;
相关错误为
??? Error using ==> run
Attempted to access pL(1); index out of bounds because numel(pL)=0.
不知道问题出在哪里了？查了好久了，找不到缘由，能不能帮帮我啊，非常感谢！