from calcephpy import *

def printcoord(PV,name):
    print('{0} :\n{1}\n'.format(name,PV))

jd0=2442457
dt1=0.5E0
dt2=0.9E0

peph = CalcephBin.open("example1.dat")

PV1 = peph.compute(jd0, dt1, 4, 11)
printcoord(PV1,"heliocentric coordinates of Mars")

PV2 = peph.compute(jd0, dt2, 4, 11)
printcoord(PV2,"heliocentric coordinates of Mars")

peph.close()