We investigate through this research the numerical inversion technique for the Laplace transforms cooperated by the integration Boubaker polynomials operational matrix. The efficiency of the presented approach is demonstrated by solving some differential equations. Also, this technique is combined with the standard Laplace Homotopy Perturbation Method. The numerical results highlight that there is a very good agreement between the estimated solutions with exact solutions.