Abstract
In this research work, we discuss an approximation techniques for boundary value problems (BVPs) of differential equations having fractional order (FODE). We avoid the method from discretization of data by applying polynomials of Laguerre and developed some matrices of operational types for the obtained numerical solution. By applying the operational matrices, the given problem is converted to some algebraic equation which on evaluation gives the required numerical results. These equations are of Sylvester types and can be solved by using matlab. We present some testing examples to ensure the correctness of the considered techniques.