Solve a system of 3, 2nd order Boundary value problems

Hello, I want to solve a system of 3 boundary value equations. All three of them are 2nd order. Now I have created my dydx by converting them into 6 ODE's. Since they are huge equations (almost 30 lines) I can give you a gist of them. Baiscally its d2y1/dx +dy1/dx+dy2/dx+dy3/dx+y1+y2+y3+x = const d2y2/dx +dy1/dx+dy2/dx+dy3/dx+y1+y2+y3+x = const d2y3/dx +dy1/dx+dy2/dx+dy3/dx+y1+y2+y3+x = const

I know the boundary values of these. They are: y1(0) = 0.6; y1(20) = 0 y2(0) = 0.3; y2(20) = 0 y3(0) = 0.005; y3(20) = 0

So how do I write the bcfun and the solinit for this problem. I wrote the odefun by converting the above equations into a set of 6, 1st order ODE's.

If you can help me with the syntax ( just a small part of it) that will great

Thanks, Nitin Nair

[Prakash Dhanapal]

see matlab demo in help. It will make you clear. Put bvp4c or bvp5c in that help. you may see the syntax in there.

[mosi nezhadhasan]

hi I answer this in question.