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();