Antecedentes del problema:
Estoy tratando de hacer una simulación aproximada de la segunda ley de Kepler (áreas iguales en tiempos iguales) y para hacerlo he dividido el área de la elipse en varias partes. Quiero averiguar cuánto necesito aumentar el ángulo para que barra un arco que forme un sector que tenga un área igual al área del corte (que es (Área de elipse/n)
Estoy representando esta elipse con la ecuación polar para una elipse de:
http://mathworld.wolfram.com/Ellipse.html
donde e es la excentricidad de la elipse.
Esta pregunta:
Cómo calcular el área del sector de elipse *desde un foco*
describe resolver el área dado el ángulo, ¿cómo me acerco a esto si estoy tratando de resolver el cambio en theta mientras conozco el área que quiero resolver?
Acabo de escribir una nueva respuesta a la pregunta a la que hace referencia . Esa publicación describe cómo puedes salir de una verdadera anomalía. vía anomalía excéntrica significar anomalía que es proporcional tanto al área como al tiempo.
Para encontrar el cambio en la anomalía verdadera correspondiente a un cambio dado en el área, primero giraría el ángulo original en alguna anomalía media , luego agregaría a eso e intentaría revertir el proceso, es decir, calcular de . Pero la relación entre y es la ecuación de Kepler :
Esta es una ecuación trascendental. Wikipedia escribe :
Esta ecuación no tiene una solución de forma cerrada para dado . Suele resolverse mediante métodos numéricos, por ejemplo, el método de Newton-Raphson .
El lema sobre la ecuación de Kepler tiene una sección completa sobre este problema inverso , seguida de una sobre aproximaciones numéricas. Así que supongo que ese sería el camino a seguir.