An expression for the first-passage density of Brownian motion to a curved boundary is expanded as a series of multiple integrals. Bounds are given for the error due to truncation of the series when the boundary is wholly concave or wholly convex. Extensions to the Brownian bridge and to continuous Gauss–Markov processes are given. The series provides a practical method for calculating the probability that a sample path crosses the boundary in a specified time-interval to a high degree of accuracy. A numerical example is given.