m. slade  feb. 1970           subroutine mensun
computation of approximate  solar position in rectangular
coordinates (equator and equinox of ref.epoch) for use in phase
corrections and  star occultations by moon

if k=-1, velocity only
if k=0, position only
if k=1, position and velocity
distances in a.u.
velocity in  a.u./ephemeris days