Solucionador de Órbita de Kepler

Dada a anomalia média M (graus), a excentricidade e e o semieixo maior opcional a (UA), a ferramenta resolve a equação transcendente

M = E − e · sin(E)

para a anomalia excêntrica E. A partir de E derivam-se:

• Anomalia verdadeira ν: tan(ν/2) = √((1+e)/(1−e)) · tan(E/2)
• Raio orbital: r = a · (1 − e·cos E)
• Velocidade relativa (vis-viva): v/v_circ = √(2 − r/a)

Três métodos:

• Newton-Raphson (recomendado): derivada analítica f′(E) = 1 − e·cos E, convergência quadrática. Para e alto (>0.8), início em E = π.
• Bissecção: robusto, delimitado em [M−e, M+e], convergência linear.
• Ponto fixo E ← M + e·sin E: o mais simples, mas diverge quando e → 1.

O relatório mostra a órbita em SVG (elipse centrada no foco, com periélio, afélio, foco e posição atual com seu vetor raio), as grandezas, a classificação da órbita (circular/elíptica/parabólica/hiperbólica) e o histórico completo de iterações com resíduos.