I have a system of spatial ODEs to solve... Actually a DAE system, but here's the issue: The equations are vaild over a specific domain, x = 0..L The equations are only bound at one point (thier "initial point") but at either 0 or L f1(0)=0 f2(0)=100 f3(L)=0 f4(L)=100 (also an algebraic expression that links all of the functions) Essentially, those functions bound at L are moving "backwards" with respect to those bound at 0. My question is, is this an initial value problem, or a boundary value problem? I started in MAPLE, and it decided (automatically) that it was a BVP. In MATLAB I am attempting to use bvp4c to come up with a solution, but is it really even a BVP if it only has one boundary condition, essentially an initial conditition, but defined at different spatial coordinates for different functions? Your insight is appreciated.