Modeling of the radial and transverse components of the albedo function

In this work, we shall confine ourselves to solve analytically the integrals called , for the two components of albedo radiation pressure on an Earth’s satellite. When the Earth’s albedo is variable, as far as we know, this case has not been dealt with analytically. We shall solve these two integrals when the satellite’ horizon is illuminated and when the sun lies on the equator. This approach will facilitate the evaluation of the mentioned two equations. We will also compare our results with previous works obtained numerically.