We have calculated the astronomical almanac 1994 and simulated the trajectory of a satellite orbit considering all perturbative forces with various initial conditions. In this work, Gauss Jackson multistep integration method has been used to calculate our basic equation of motion with high numerical accuracy. It has beer. found that our results agree well with the Astronomical Almanac Data distributed by JPL of NASA and the orbit simulations have been carried out with fast speed, stability and excellent round-off error accumulation, comparing with other numerical methods. In order to be carried out our works on almanac and orbit calculations easily by anyone who uses a personal computer, we have made a computer program on graphical user interface to provide various menus for detail works selected by a mouse.