jd0=2442457
dt1=0.5E0
dt2=0.9E0
peph = CalcephBin.open('example1.dat');
disp('heliocentric coordinates of Mars')
PV1 = peph.compute(jd0, dt1, 4, 11)
disp('heliocentric coordinates of Mars')
PV2 = peph.compute(jd0, dt2, 4, 11)
peph.close();
jd0=2442457
dt1=0.5E0
dt2=0.9E0
peph = CalcephBin.open('example1.dat');
disp('heliocentric coordinates of Mars')
PV1 = peph.compute(jd0, dt1, 4, 11)
disp('heliocentric coordinates of Mars')
PV2 = peph.compute(jd0, dt2, 4, 11)
peph.close();