A novel numerical convolution of rational Laplace transform inversion

A novel numerical convolution method is proposed as an alternative to conventional inversion by Heaviside´s expansion of the rational Laplace transforms. The approach is based on divided differences of inverse Laplace transforms. This new convolution method can be easily extended to the case when the Laplace transforms contain repeated poles. Numerical examples are given to demonstrate their application in educational learning and numerical tool development