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://www.roneverhart.com/Wiley-X-CAG-1-Gloves/