goose - the square root of (g*(mass of body+mass of center))
cond - array of initial conditions + optional extras
cond(7)=time derivative of cond(1)
cond(8)=time derivative of cond(5)
cond(9)=mass factor for scaling output
Usage of array 'setp' --- storage for saved quantities from
setting up.  It is effectively equivalenced to the following:
b(3,2),tsv,a,e,anom0,motpi,mu2,secc,cecc,quan2,quan3,quan4,quan5,
1      7   8 9  10    11    12  13   14   15    16    17    18

kind - indicates type of setup,
0 => only want coordinates
1 => also want partials w.r.t. initial conditions
2 => also need sinc,cinc,sper,cper in extended setp
3 => also set up for moving periapse and semimajor axis
y - array for output coordinates
ry,ry2,ry3 - output radial distance, square, and cube
dy - output array of partials

Compute elliptic/hyperbolic quantities for time t from initial epoch.
Corresponding setp array must already have been set up.
nv - operation indicator.
1 do position and velocity
>1 do partials as well (up to dy(1-6,nv-1) )
0 do just position
-1 do just velocity, in y(1-3)
-2 do just acceleration, in y(1-3)