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.

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.