How to solve second order coupled ODEs with bvp4c?

I want to solve the following ODEs:

With boundary conditions

I'm aware that this problem needs to be solved with bvp4c or bvp5c, although I have great difficulty solving this because the equations are coupled. Is it possible to solve this problem with bvp4/5c, and how do I couple this?

Answer 1

[Jan]

function dy = YourODE(t, y)
k  = 17;
dy = [y(3); ...
      y(4);
      k * y(1) * y(2); ...
      k * y(1) * y(2)];
end

If you call the coordinates x, y or c_a, c_b does not matter. I do not see, why a coupling is a problem. Does this help you already?