We consider the numerical treatment of two-point boundary value problems of fractional differential equations with Caputo's derivatives of order $\gamma$ where $1<\gamma <2.$ We introduce the fractional Green functions in order to prove the existence and uniqueness of the solutions to the problem. Numerically, we introduce the simple shooting method for linear and nonlinear problems. The results obtained are accurate and competes well with results by other authors using different approaches.