In this study, we present a novel numerical scheme for the approximate solutions of linear as well as non-linear ordinary differential equations of fractional order with boundary conditions. This method combines Cosine and Sine (CAS) wavelets together with Green function, called Green-CAS method. The method simplifies the existing CAS wavelet method and does not require conventional operational matrices of integration for certain cases. Quasilinearization technique is used to transform non-linear fractional differential equations to linear equations and then Green-CAS method is applied. Furthermore, the proposed method has also been analyzed for convergence, particularly in the context of error analysis. Sufficient conditions for the existence of unique solutions are established for the boundary value problem under consideration. Moreover, to elaborate the effectiveness and accuracy of the proposed method, results of essential numerical applications have also been documented in graphical as well as tabular form.