Finite-time ruin probabilities using bivariate Laguerre series

In this paper, we revisit the finite-time ruin probability in the classical compound Poisson risk model. Traditional general solutions to finite-time ruin problems are usually expressed in terms of infinite sums involving the convolutions related to the claim size distribution and their integrals, which can typically be evaluated only in special cases where the claims follow exponential or (more generally) mixed Erlang distribution. We propose to tackle the partial integro-differential equation satisfied by the finite-time ruin probability and develop a new approach to obtain a solution in terms of bivariate Laguerre series as a function of the initial surplus level and the time horizon for a large class of light-tailed claim distributions. To illustrate the versatility and accuracy of our proposed method which is easy to implement, numerical examples are provided for claim amount distributions such as generalized inverse Gaussian, Weibull and truncated normal where closed-form convolutions are not available in the literature.