Hi,

I am tring to solve two coupled pde's with piecewise/Heaviside conditions. I am puzzled with the maple error...Rhop, Rhos are constants n given some values..

> pde_sys:=[diff(c(m,t),t)=(2/m)*diff(c(m,t),m)+diff(c(m,t),m,m)-(convert(piecewise(d(m,t)<Rhop,(Rhop-d(m,t))*d(m,t)*c(m,t)),Heaviside)+(Rhos-d(m,t))*c(m,t)*convert(piecewise(d(m,t)>Rhop,Rhop,d(m,t)),Heaviside)), diff(d(m,t),t)=convert(piecewise(d(m,t)<Rhop,(Rhop-d(m,t))*d(m,t)*c(m,t)),Heaviside...

Hello all,

I am trying to solve a pde system with piecewise conditions numerically but unsucessful. Here's what i am doing:

> restart;

> sys:=[diff(c(x,t),t)=v(x,t)*diff(c(x,t),x$2),diff(v(x,t),t)=piecewise(c(x,t)>v(x,t),c(x,t)*diff(v(x,t),x$2),0)];

> bc:={c(0,t)=2,c(1,t)=1,c(x,0)=0,v(0,t)=1,v(1,t)=2,v(x,0)=0};

> sol:=pdsolve(sys,bc,type=numeric,time=t);

I am getting the following error.......

Error, (in pdsolve/numeric/process_PDEs...