// ***** ALGORITHME 31 ***** /// <summary> /// Calcule l’anomalie vraie et le rayon vecteur des comètes périodiques rapportés à l’équinoxe J2000 à partir de l’anomalie moyenne /// des comètes, ainsi que de l’excentricité et du périhélie des orbites des comètes. /// </summary> private void CalculerAnomalieVraieEtRayonVecteur() { // Déclaration des variables de la méthode double E; // variable de calcul // Calcul de l'anomalie excentrique en radians E = Maths.CalculerAnomalieExcentriqueAvecKelper(anomalieMoyenne.Decimale, parametresOrbitaux.Exentricite); // Calcul de l'anomalie vraie anomalieVraie = new Angle(2 * Maths.RadToDeg(Math.Atan(Math.Sqrt((1 + parametresOrbitaux.Exentricite) / (1 - parametresOrbitaux.Exentricite)) * Math.Tan(E / 2))), TypeAngle.ANGLE_DEGRES_360); // Calcul du rayon vecteur rSoleil = (parametresOrbitaux.Perihelie / (1 - parametresOrbitaux.Exentricite)) * ((1 - parametresOrbitaux.Exentricite * parametresOrbitaux.Exentricite) / (1 + parametresOrbitaux.Exentricite * Math.Cos(Maths.DegToRad(anomalieVraie.Decimale)))); }