Hi everyone,

I am trying to solve the following dimensionless system of pdes:

u1t=u1x+u2-u1

u2t=3.139uxx-2u2+u3+u1

u3t=-u3x+u2-u3

the boundary conditions are u1(5,t)=0, u2(0,t)=0, u2(5,t)=0, u3=1-heaviside(t-0.1)

initial conditions are u1(x,0)=0, u2(x,0)=0, u3(x,0)=0

pdsolve returns some results but they are not stable. Since I'm fairly new to maple, I don't know if it is possible to discretize...