Question: single PDE solver error

Dear all,

I'm a new maple user trying to solve a single PDE in maple using finite volume.

For simplicity I'm using three volumes. When I'm trying to apply a simple flux limiter, up-wind scheme, (relation between the wall and the node values u[j+0.5](t)=u[j](t)), I'm solver works well and produces accurate results.

However, when I try to a more complicated flux limiter, the solver doesn't initialize and I'm getting an error:
"Error, (in dsolve/numeric/DAE/make_proc) specified dependent variables [u[0](t), u[.5](t), u[1](t), u[1.5](t), u[2](t), u[2.5](t), u[3](t), u[3.5](t), u[4](t)] do not agree with input system {u[0], u[1], u[2], u[3], u[4], u[1/2], u[1.5], u[2.5], u[3.5]}, differences: {u[1/2]}, {u[.5]}".

I would greatly appreciate it if somebody helps me on in this issue.

error.mw

Please Wait...