用pdepe的方法解求解方程组,有二个方程是PDE方程,有一个就是丫普通线性的。
我一开始用pdepe解的时候,直接把这个非PDE方程,定义成c=0,f=0,s=f(u).这种形势,可是方程组运行的时候就有毛病,怎么调也不行,是不是这种方法不行啊
下面是错误提示:
??? Error using ==> funfun\private\daeic12
This DAE appears to be of index greater than 1.
Error in ==> ode15s at 394
[y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...
Error in ==> pdepe at 289
[t,y] = ode15s(@pdeodes,t,y0,opts,pde,ic,bc,m,xmesh,xi,xim,zxmp1,xzmp1,varargin{:});
>>
求了,求了求了求大神
我一开始用pdepe解的时候,直接把这个非PDE方程,定义成c=0,f=0,s=f(u).这种形势,可是方程组运行的时候就有毛病,怎么调也不行,是不是这种方法不行啊
下面是错误提示:
??? Error using ==> funfun\private\daeic12
This DAE appears to be of index greater than 1.
Error in ==> ode15s at 394
[y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...
Error in ==> pdepe at 289
[t,y] = ode15s(@pdeodes,t,y0,opts,pde,ic,bc,m,xmesh,xi,xim,zxmp1,xzmp1,varargin{:});
>>
求了,求了求了求大神
