JPL site http://ssd.jpl.nasa.gov/horizons.cgi
I do not now remember why I wrote this for. I think it was an initial attempt in Matlab, because in the final report I used Mathematica. But here is the listing. It seems to be do something. I should make an animation of this.
nma_project2_EMA550_v3.m
nma_project2_EMA550_driver