In this paper we determined the variation of the Keplerian orbital elements of a GPS satellite due to the direct and indirect action of solar radiation pressure. For this study, we created a soft program to determine the Keplerian elements. This soft uses the initial conditions of position and speed of a GPS satellite and solves the Laplace integration problem using Runge - Kutta algorithm of 4