Using the expression of the orbit of an ellipse of small eccentricity \(e\) in polar coordinates, show that the variation of the planet's angular velocity as a function of eccentricity is equal to \(\frac{\Delta \omega}{\omega} \simeq 4 e\).