In the study of bivariate processes, one often encounters expressions involving repeated combinations of bivariate continuum operations such as multiple bivariate convolutions, marginal convolutions, tail integration, partial differentiation and multiplication by bivariate polynomials. In many cases numerical computation of such results is quite tedious and laborious. In this paper, the bivariate Laguerre transform is developed which provides a systematic numerical tool for evaluating such bivariate continuum operations. The formalism is an extension of the univariate Laguerre transform developed by Keilson and Nunn (1979), Keilson et al. (1981) and Keilson and Sumita (1981), using the product orthonormal basis generated from Laguerre functions. The power of the procedure is proven through numerical exploration of bivariate processes arising from correlated cumulative shock models.
