Calculation of refraction angles in the case of representation of the Euler equation in a spherical coordinate system

In this paper, we propose a numerical-analytical method for the approximate determination of the coordinates of the points of the light curve, when the refractive index n is a function of the polar coordinates r, Ѳ