I have a system which consists of the following coupled second order ODEBVP
d2a1/dx2 = C1*a1*a2 - C2*(1-a2^2);
d2a2/dx2 = C3*a1*a2 - C4*(1-a2^2);
The boundary conditions are
At x = 0, a1 = 1, da2/dx = 0;
At x = 1, a1 = C5, da2/dx = 0;
a1, a2 - Variables, C1-C5 - Constants
I tried using bvp4c but I am not able to solve it. Can someone please suggest me how can I solve my ststem and get plots for a1 and a2 wrt x. I have also tried using finite difference for solving the system but was not able to solve.