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): maple-shoot_RK4.mw