求助matlab pdepe解偏微分方程的问题

kevintheunique

程序如下

clc

m=0;

x=linspace(0, 1.25e-5,1000);

t= linspace(0,1000,1000);

sol=pdepe(m,@pdefun,@pdeic,@pdebc,x,t);

u1=sol(:,:,1);

u2=sol(:,:,2);

figure(1)

surf(x,t,u1)

title('u1(x,t)')

xlabel('Distance x')

ylabel('Time t')

zlabel('u1')

figure(2)

surf(x,t,u2)

title('u2(x,t)')

xlabel('Distance x')

ylabel('Time t')

zlabel('u2')

a=3.6e7;

b=4;

cc=2.6e-7;

d=1.6e18;

e=4e-13;

ff=1e10;

g=0.3;

function [c,f,s]=pdefun(x,t,u,dudx)

c=[1;1];

f=[cc*du1dx-b+cc*e*(d-u2);cc*e*(d-u2)];

s=[a-ff*u2-e*b*(d-u2)*u(1);e*b*(d-u2)*u1 -ff*u2];

function u0=pdeic(x)

u0=[0;0];

function [p1,q1,p2,q2]=pdebc(x1,u1,x2,u2,t)

p1=[b*u1;0];

q1=[c;0];

p2=[(b-g)*u1;0];

q2=[c;1];

运行时报错说

??? Error using ==> fevalUndefined function or method 'pdeic' for input arguments of type 'double'.
Error in ==> pdepe at 244temp = feval(ic,xmesh(1),varargin{:});错误应该是出在初值函数pdeic那里 但我试过改过数据类型，也不对

我看文献上都是这么写的 结果自己跑就有问题

查了好久都不知道错在哪里  求高手赐教

kevintheunique

没人用过这个么