Estoy tratando de calcular la cantidad de oblateness que es causada por la rotación planetaria. Me imagino a la fuerza de la gravedad añadida a la fuerza centrífuga causada por la rotación del planeta como sigue:
$\hspace{4cm}$
Es decir, en el punto en cuestión, a una latitud de $\phi$, la distancia desde el eje de rotación es $r\cos(\phi)$. Por lo tanto, la fuerza centrífuga sería $\omega^2r\cos(\phi)$ en una dirección perpendicular al eje de rotación. La radial y tangencial de los componentes se $\omega^2r\cos^2(\phi)$$\omega^2r\cos(\phi)\sin(\phi)$, respectivamente.
Mi suposición es que la superficie del planeta se ajustará de modo que sea perpendicular a la efectiva $g$; es decir, la suma de la gravedad y las fuerzas centrífugas. Esto llevaría a la ecuación $$ \frac{\mathrm{d}r}{r\,\mathrm{d}\phi}=-\frac{\omega^2r\cos(\phi)\sin(\phi)}{g-\omega^2r\cos^2(\phi)} $$ Podemos hacer varias suposiciones aquí, y voy a suponer que $\omega^2r$ es pequeña en comparación a $g$. Por lo tanto, obtenemos $$ \int_{\text{eq}}^{\text{pn}}\frac{\mathrm{d}r}{r^2} =-\frac{\omega^2}{g}\int_0^{\pi/2}\cos(\phi)\sin(\phi)\,\mathrm{d}\phi $$ lo que conduce a $$ \frac1{r_{\text{pn}}}-\frac1{r_{\text{eq}}} =\frac{\omega^2}{2 g} $$ y $$ 1-\frac{r_{\text{pn}}}{r_{\text{eq}}} =\frac{\omega^2r_{\text{pn}}}{2g} $$ Sin embargo, la evaluación numérica y Wikipedia parecen indicar que este debe ser el doble de lo que estoy consiguiendo. Es decir, $$ 1-\frac{r_{\text{pn}}}{r_{\text{eq}}} =\frac{\omega^2r^3}{Gm} =\frac{\omega^2r}{g} $$ ¿Qué estoy haciendo mal?