For numerical methods in optimal control problems whose dynamics are represented by ordinary differential equations, theories are of indirect and gradient method are presented with variational approach. Theories developed are applied to the problem of determination of optimal trajectory injecting to a parabolic orbit. Two-point boundart-value problem resulted from application of indirect method is solved using initial value method and neighboring optimal methods.