I suggest you study some references of numerically solving such differential equations (numerical calculations or numerical analysis). The next step is compiling the required steps in software.

Generally, numerical solution of such a system, needs some generic steps:

1. Initial conditions and boundary conditions

2. Convergence criterion

3. Generating next solution using recursive methods and increments

Dear Ma Ka

I think, you can never solve these set of equations, using Maple pre-defined command syntax. You should first make a digestion of the available solution strategies as applied to the equation, then some remaining parts of the analysis could be expected to be performed by Maple. I repeat again that, these hardly nonlinear system of coupled ODEs with two-point boundary-value problem setup, require much human thoughts and mind digestion, before being fed into a computer ODE solver.

Hello

Thank you very much for your valuable tips.

hi, first of all, it is better Mathematica or Matlab( i prefer Matlab).

second, you should use explicit numerical methods like rung_kutta 4th, rkf, adams_boshforth, or predictor _ corrector method (adams_boshford_moulton).

ode15s,ode23s can be helpful in solving nonlinear stiff odes sys.

Hello

Thank you very much for your comment.

You are right. It is better to solve by MATLAB, But I did not write much program with it. Also, I do not have much time.

Dear Ma Ka

I suggest you to start your excursion into Matlab, by learning it from a class fellow or friend who is conversant with Maple or Matlab languages.

Learning Matlab by self study is a bit difficult, however, interactive videos available online, are as well a good option.