Using BVP4C for overdetermined system

[ktsharp]

I have an overdetermined nonlinear system of ODEs:

W' = f(c)
c'' = f(W,W',c)

and boundary conditions
W(0)=a,W(L)=-a
c(0)=c(L)-b

I can split up the equations into three first order ODEs, and solve numerically with Matlab. I would like to use bvp4c, but I believe I have too many boundary conditions. Is this correct? If so, what alternative approach should I take to finding a solution. And if not, what am I doing wrong? I have already implemented the shooting method which is very sensitive on ICs for convergence, so looking for alternative method.

Regards,

 

[ktsharp]

I've realized that I haven't declared any parameters, and thus introduce an extra dimension. This was the problem!