Question: Solving ODE - shooting method, Runge-Kutta method of 4th order

Hi guys,

I have got serious problem with solving this system of ODE, where psi is equal to 10*sqrt(da) and uB,E,kL,uc are constants:

firstly I have to find missing initial condition using shooting method and calculate cA(z=2) using Runge-Kutta 4th order then. And plot following concentration profile for cA and dA. So far, I have dealt with first diff eq 2nd order dividing by two diff eq of 1st order using mentioned constants so I got these ones:


I understand I need a condition da(z=0) that I should obtain using shooting method, but I do not know how to do it in spite of I understand it in theoretical way. And same problem I have with RK4th order. Anybody here with a hint, please? 

There is my maple file (unfinished version of mine):

Please Wait...