This paper concerns linear and nonlinear n^th order boundary value problems that contain a finite number of random variables in the boundary conditions or in the differential equation. The results extend methods previously known for corresponding initial value problems. Numerically implementable procedures are given for the determination of the joint density of the solution at an arbitrary point. The possible use of Liouvilleapos;s equation to reduce a random boundary value problem to a random initial value problem is also indicated.