ASTRONOMIE

Tout est dans le ZIP garanti sans virus. ce n'est pas un exe mais juste un .bas
Pour vous donner une idée voici le code d'une des fonctions

'cette fonction renvoie la longitude de la planète en radian.
Function Calcul_Longitude(Planete As Long, JJ As Double) As Double
    'Pour calculer la longitude (l) par rapport à l'écliptique, il faut connaître
    'i, l'inclinaison
    'u, la latitude de la planète
    'Omega, la longitude du noeud ascendant
    'tan(l - Omega)=cos(i)*tan(u)   formule 21.4 page 78
    
    'Pour cela on utilise les fonctions correspondantes qui dépendent de la planète et du
    'temps écoulé de puis le 1er janvier 1900 à minuit (T) ainsi que de T^2 et T^3
    
    Dim T As Double
    T = (JJ - 2415020#) / 36525
    
    Dim T_2 As Double
    T_2 = T * T
    
    Dim T_3 As Double
    T_3 = T * T_2
    
    If Planete = Lune Then
        Calcul_Longitude = Calcul_Longitude_Lune(T, T_2, T_3)
        Exit Function
    End If
    
    'L'inclinaison
    'i=a0+a1*T+a2*T_2+a3*T_3 tableau 20a, page 73
    Dim i As Double
    i = Inclinaison(Planete, T, T_2, T_3)
    
    'u, L'argument (?) de la latitude de la planète dépend de
    'LM, la longitude moyenne
    'AV, l'anomalie vraie
    'AM, l'anomalie moyenne
    'Omega, la longitude du noeud ascendant
    'u=LM+AV-AM-Omega   formule 21.3 page 78
    
    'LM, la longitude moyenne
    'LM=a0+a1*T+a2*T_2+a3*T_3 tableau 20a, page 73 (dans le tableau il s'agit de L)
    Dim LM As Double
    LM = LongitudeMoyenne(Planete, T, T_2, T_3)
    
    'AV, l'anomalie vraie dépend de l'excentricité e et de l'anomalie excentrique AE
    'tan(AV/2)=Racine((1+e)/(1-e))*tan(AE/2)    formule 21.1 page 78
    
    'e, l'excentricité
    'e=a0+a1*T+a2*T_2+a3*T_3 tableau 20a, page 73
    Dim e As Double
    e = Excentricite(Planete, T, T_2, T_3)
    
    'AE, l'anomalie excentrique dépend de
    'e, l'excentricité
    'AM, l'anomalie moyenne
    'Elle se trouve en résolvant l'équation de Kepler
    '2ème méthode modifiée pour obtenir des radians, page 69
    'AE1=AE0+(AM+e*sin(AE0)-AE0)/(1-e*cos(AE0))
    
    'AM, l'anomalie moyenne dépend de
    'AM=a0+a1*T+a2*T_2+a3*T_3 page 77 (55 pour la terre)
    
    'Omega, la longitude du noeud ascendant
    'Omega=a0+a1*T+a2*T_2+a3*T_3 tableau 20a, page 73 (dans le tableau il s'agit de grand oméga)
    Dim Omega As Double
    Omega = LongitudeNoeudAscendant(Planete, T, T_2, T_3)
    
    Dim AM As Double
    AM = AnomalieMoyenne(Planete, T, T_2, T_3)
    
    Dim AE As Double
    AE = AnomalieExcentrique(e, AM)
    
    Dim AV As Double
    AV = AnomalieVraie(e, AE)
    
    Dim u As Double
    u = latitude(LM, AV, AM, Omega)
    
    Dim L As Double
    L = LongitudeEcliptique(i, u, Omega)
    
    Calcul_Longitude = L

End Function