Hola a todos
En esta ocasión voy a tratar de responder a Richard para explicar qué es lo que ha realizado de forma incorrecta en el desarrollo que ha realizado.
1) Nos planteas la función F(x), la fuerza de gravedad que se ejerce sobre una masa m, por una esfera hueca de espesor infinitesimal, siendo x la distancia entre el punto a medir y la superficie exterior de la esfera, a lo largo de una dirección radial.
En un momento determinado, después de resolver la integral se obtiene dentro de la expresión algebraica general un término \( \frac{1}{x} \), el siguiente:
\( \frac{1}{2R+x}-\frac{1}{x} \)
Para cualquier valor de \( x\neq 0 \), no hay ningún problema pero para x=0, la función F(x) no está definida.
Por lo tanto en el momento que escribes \( \frac{1}{x} \) en ese preciso instante, se debe expresar claramente que esta expresión es para todo \( x\neq 0 \). Y todo lo que obtengas después está implícito que \( x\neq 0 \).
Por lo tanto lo que tú has obtenido en tus ecuaciones es sólo válido para \( x\neq 0 \) y la magia final que has hecho es falsa y errónea para x=0.
2) Una vez que observas que en x=0 (r=R), no tenemos un valor definido de la función gravedad, el siguiente paso más probable por las personas que lo descubre es ?
Yo me plantee saber cuánto vale la gravedad para x=0 o r=R.
Pero eso, nadie lo ha realizado desde hasta 300 años hasta este momento que yo sepa y ahora sí se ha realizado por varias personas: Por Masacroso en los inicios del hilo, también lo ha realizado Feriva cuando ha realizado el esquema geométrico en donde llevó la posición r=R, también lo he indicado yo mismo explicando que en el desarrollo de la página original que puse, partimos de r=R (punto sobre la superficie exterior), en la integral que has puesto el término \( \frac{r^2-R^2}{s^2} \) desaparece y no está. Eso mismo obtuvo Feriva y este término desapareció. En todos los casos siempre ha salido el mismo resultado que tanto he compartido con todos:
para r=R \( g=\frac{GM}{2R^2} \) o lo que es lo mismo \( g=2\pi G \sigma \)
A continuación voy a hacer el cálculo de la gravedad de la esfera hueca de espesor infinitesimal para r=R, basándome en el siguiente enlace ya conocido:
http://hyperphysics.phy-astr.gsu.edu/hbasees/Mechanics/sphshell.html#wtls DIBUJO 1

Como podéis ver en esta página la esfera se ha descompuesto en infinitos anillos masivos de anchura diferencial \( d\theta \) y cada anillo está situado respecto del punto de medida a la distancia r del centro de la esfera con 2 ángulos \( \theta \) y \( \alpha \) así como por la distancia s.
La expresión general de la fuerza de la gravedad que genera la esfera hueca de espesor infinitesimal en la masa m es la siguiente:
\( F=2\pi G\sigma mR^2\int_{\theta=0}^{\pi}\frac{1}{s^2}cos\alpha sin\theta d\theta \)
A partir de ahora vamos a obtener el valor de F para r=R (en la superficie de la esfera exterior) y como vemos existe un anillo que corresponde a \( \theta = 0 \), que s=0 y por ello a partir de ahora supondremos que \( s\neq 0 \), para no tener problemas con el término \( \frac{1}{s^2} \) y la función esté siempre definida.
Por ello
\( F=2\pi G\sigma mR^2\lim_{\epsilon \to 0}\int_{\epsilon}^{\pi}\frac{1}{s^2}cos\alpha sin\theta d\theta \) (a)
Aplicando la ley del coseno para \( \theta \) para cualquier punto r:
\( s^2=R^2+r^2-2Rrcos\theta \) sustituyendo para r=R y diferenciando obtenemos:
\( 2sds=2R^2sin\theta d\theta \) despejando
\( sin\theta d\theta=\frac{sds}{R^2} \) (b)
Aplicando la ley del coseno para \( \alpha \) en general para cualquier punto r:
\( R^2=r^2+s^2-2r scos\alpha \) sustituyendo r=R obtenemos:
\( s^2=2R scos\alpha \)
\( s cos\alpha=\frac{s^2}{2R} \) (c)
Entonces sustituyendo (b) en la ecuación (a) obtenemos :
\( F=2\pi G\sigma mR^2\lim_{\epsilon \to 0}\int_{\epsilon}^{2R}\frac{1}{s^2}cos\alpha \frac{sds}{R^2} \) y utilizando el resultado de (c) obtenemos:
\( F=\frac{2\pi G\sigma m}{2R}\lim _{\epsilon\to 0 }\int_{\epsilon}^{2R}ds \) dónde vemos que el término \( \frac{r^2-R^2}{s^2} \) se ha eliminado por completo para r=R.
Finalmente obtenemos que para r=R que la gravedad es:
\( F=2\pi G\sigma m \) o lo que es lo mismo \( F=\frac{GMm}{2R^2} \)
En resumen, llevamos más de 300 años expresando de forma inexacta el resultado de la integral que estamos estudiando, pues en el resultado de la integral de este enlace y en cualquier otro que busques, libros, etc, se les olvidó decir que todo lo dicho, no es válido para r=R, pero esto en matemáticas es un grave error, por lo que mi conclusión es que la integral se ha realizado incorrectamente al faltar este comentario.
\( g=\frac{GM}{r^2} \) para todo r perteneciente a [R,infinito)
ES FALSO pues como mínimo para r=R \( g=\frac{GM}{2R^2} \)
Como ya he explicado en este hilo, este resultado cambia por completo el comportamiento de la gravedad de las esferas huecas de espesor infinitesimal, para los puntos que estén a una distancia del orden del radio de la esfera respecto a lo que oficialmente conocemos y deja invariante al resto de puntos más alejados.
Y por consiguiente afectará al comportamiento de la gravedad de las esferas huecas de espesor finito.
Un saludo.