1-I got the results by using BVP4c but when I tried to change or manipulate the boundary conditions for different configurations, I got meaningless results.
2-I have difficulty evaluating when using the ODE45 solver, considering that I do not use time as a variable.
I revised the coding many times, and I got negative values of parameters which is wrong process.