int j;
double jd0=2442457;
double dt1=0.5E0;
t_calcephbin *peph;
double PVAJ[12];
double P[3];
/* open the ephemeris file */
peph = calceph_open("example1.dat");
if (peph)
{
/* compute only the heliocentric position of Mars in km */
calceph_compute_order(peph, jd0, dt1,
NAIFID_MARS_BARYCENTER,
NAIFID_SUN,
CALCEPH_USE_NAIFID+CALCEPH_UNIT_KM
+CALCEPH_UNIT_SEC,
0, P);
for(j=0; j<3; j++) printf("%23.16E\n", P[j]);
/* compute positions, velocities, accelerations and jerks
of Mars in km and seconds */
calceph_compute_order(peph, jd0, dt1,
NAIFID_MARS_BARYCENTER,
NAIFID_SUN,
CALCEPH_USE_NAIFID+CALCEPH_UNIT_KM
+CALCEPH_UNIT_SEC,
3, PVAJ);
for(j=0; j<12; j++) printf("%23.16E\n", PVAJ[j]);
/* close the ephemeris file */
calceph_close(peph);
}