Investigation of the step variable numerical integration methods is the main purpose of this paper. Also it is shown how the one-dimensional numerical integration method of solving an ordinary differential equation can be employed to solve a three-dimensional differential equation numerically. In order to integrate the orbit, the step variable methods of Runge-Kutta-Fehlberg and Adams are employed and their accuracies are investigated.