laplace_orbit_fit() implements the Laplace method for orbit determination from 3 distinct azimuth/elevation observations of a body. The underlying theory is described in Bate, White and Mueller. Note that the observations must be closely spaced in or