Is there a Mathematica package/routine to solve O.D.E's, numerically,
in a boundary value problem?
(i.e. conditions of the form :
y[x1] == a1
y'[x2] == a2
...
where x1 != x2)
I've tried NDSolve but get an error message :
NDSolve::ndinn: Initial condition y[-1000.] is not a number.
David Erwin