The shooting method is a widely used numerical technique for solving boundary value problems (BVPs). The basic idea involves transforming the BVP into an equivalent initial value problem (IVP), which is easier to solve using standard numerical methods like Runge-Kutta. The process starts with an initial guess for the missing initial conditions at one boundary. The system of differential equations is then solved as an IVP over the domain, using these guessed conditions. At the other boundary, the calculated solution is compared with the prescribed boundary conditions. If the boundary conditions are not satisfied, the initial guess is adjusted iteratively. Techniques such as the secant method or Newton’s method can be used to systematically refine the guesses. This process is repeated until the solution meets the boundary conditions within a specified tolerance, yielding an approximate solution to the BVP.In this study we considered as a number of examples to illustrate the shooting method and the solution obtained by this method are compaired with exact solution. Finally we observe the convergency, consistency and stability of the of the result with the varation of the mesh size ℎ.