Two example applications are provided:
The calculations call functions from the IAU SOFA library.
The EPH ephemerides are with respect to ICRS axes as implemented in JPL DE405.
The adoption of VSOP2010 for EPH was for consistency with DE405, the more recent VSOP2013 having instead been fitted to INPOP10a (IMCCE, Observatoire de Paris). Users need to be aware that although the agreement between the different numerically integrated ephemerides is generally very good—for geocentric directions, milliarcsecond level—this falls off by 2-3 orders of magnitude for the outer planets.
The user guide is eph.pdf and the files are in eph.zip