A numerical inverse Laplace transform method is established using Bernoulli polynomials operational matrix of integration. The efficiency of the method is demonstrated through some standard nonlinear differential equations: Duffing equation. Van der Pol equation. Blasius equation and jerk equation. https://yeswellneses.shop/product-category/facial-serums/