Boundary-value problems (BVPs) for seismic rays generally have multipl
e solutions. In practical applications the number of solutions can be
large. The algorithm presented below solves a one-parameter family of
BVPs and makes it easy to obtain all the solutions of a BVP.