This study is concerned with the development and application of a quasi-Newton descent algorithm for solving finite dimensional optimization problems subject to nonlinear equality and inequality constraints by the use of gradient projection and constraint restoration. The algorithm is applied to the positioning of a geostationary satellite which is interpreted as a problem of optimal TV-impulse orbital rendezvous.