In this paper, a Chebyshev-collocation spectral method is developed for Volterra integral equations (VIEs) of second kind with weakly singular kernel. We first change the equation into an equivalent VIE so that the solution of the new equation possesses better regularity. The integral term in the resulting VIE is approximated by Gauss quadrature formulas using the Chebyshev collocation points. The convergence analysis of this method is based on the Lebesgue constant for the Lagrange interpolation polynomials, approximation theory for orthogonal polynomials, and the operator theory. The spectral rate of convergence for the proposed method is established in the L∞-norm and weighted L2-norm. Numerical results are presented to demonstrate the effectiveness of the proposed method.

We study inversion of the spherical Radon transform with centres on a sphere (the data acquisition set). Such inversions are essential in various image reconstruction problems arising in medical, radar and sonar imaging. In the case of radially incomplete data, we show that the spherical Radon transform can be uniquely inverted recovering the image function in spherical shells. Our result is valid when the support of the image function is inside the data acquisition sphere, outside that sphere, as well as on both sides of the sphere. Furthermore, in addition to the uniqueness result, our method of proof provides reconstruction formulas for all those cases. We present a robust computational algorithm and demonstrate its accuracy and efficiency on several numerical examples.

A Legendre-collocation method is proposed to solve the nonlinear Volterra integral equations of the second kind. We provide a rigorous error analysis for the proposed method, which indicate that the numerical errors in L2-norm and L∞-norm will decay exponentially provided that the kernel function is sufficiently smooth. Numerical results are presented, which confirm the theoretical prediction of the exponential rate of convergence.

A new computationally simple, speedy and accurate method is proposed to construct first-passage-time probability density functions for Gauss–Markov processes through time-dependent boundaries, both for fixed and for random initial states. Some applications to Brownian motion and to the Brownian bridge are then provided together with a comparison with some computational results by Durbin and by Daniels. Various closed-form results are also obtained for classes of boundaries that are intimately related to certain symmetries of the processes considered.

In this paper we study a Volterra integral equation of the second kind, including two arbitrary continuous functions, in order to determine first-passage-time probability density functions through time-dependent boundaries for time-non-homogeneous one-dimensional diffusion processes with natural boundaries. These results generalize those which were obtained for time-homogeneous diffusion processes by Giorno et al. [3], and for some particular classes of time-non-homogeneous diffusion processes by Gutiérrez et al. [4], [5].

The first-crossing-time problem through two time-dependent boundaries for one-dimensional diffusion processes is considered. The first-crossing p.d.f.'s from below and from above are proved to satisfy a new system of Volterra integral equations of the second kind involving two arbitrary continuous functions. By a suitable choice of such functions a system of continuous-kernel integral equations is obtained and an efficient algorithm for its solution is provided. Finally, conditions on the drift and infinitesimal variance of the diffusion process are given such that the system of integral equations reduces to a non-singular single integral equation for the first-crossing-time p.d.f.

The first-passage-time p.d.f. through a time-dependent boundary for one-dimensional diffusion processes is proved to satisfy a new Volterra integral equation of the second kind involving two arbitrary continuous functions. Use of this equation is made to prove that for the Wiener and the Ornstein–Uhlenbeck processes the singularity of the kernel can be removed by a suitable choice of these functions. A simple and efficient numerical procedure for the solution of the integral equation is provided and its convergence is briefly discussed. Use of this equation is finally made to obtain closed-form expressions for first-passage-time p.d.f.'s in the case of various time-dependent boundaries.

The distribution of total claims payable by an insurer is considered when the frequency of claims is a mixed Poisson random variable. It is shown how in many cases the total claims density can be evaluated numerically using simple recursive formulae (discrete or continuous).

Mixed Poisson distributions often have desirable properties for modelling claim frequencies. For example, they often have thick tails which make them useful for long-tailed data. Also, they may be interpreted as having arisen from a stochastic process. Mixing distributions considered include the inverse Gaussian, beta, uniform, non-central chi-squared, and the generalized inverse Gaussian as well as other more general distributions.

It is also shown how these results may be used to derive computational formulae for the total claims density when the frequency distribution is either from the Neyman class of contagious distributions, or a class of negative binomial mixtures. Also, a computational formula is derived for the probability distribution of the number in the system for the M/G/1 queue with bulk arrivals.