Calcul de la position d'un satellite

Bnojour,

je voudrai bien réaliser une application qui permet de visualiser en 3D la terre et aussi les orbites des satellites autour de cette dernière.

Mon problème est le suivant : supposant que la terre est déjà modélisée dans un rendu 3D et qui tourne autour d'elle, ce que je cherche c'est j'ai les fonctions qui permettent de calculer la position du satellite à chaque instant t depuis un temps de début t0, et la je ne sais pas comment faire pour intégrer ce calcul vu que j'ai la date système .

par exemple pour le calcule de la position sur l'axe x je dois utiliser la fonction suivante :

long double Satellite::calculX(const long double tps) const
{
long double res;
long double e=calculAnomalieExcentrique(tps);
res=getDemiGrandAxe()*(cos(e)-getExcentricite())*
(
(cos(getArgument())*cos(getAscension()))
-
(sin(getArgument())*cos(getInclinaison())*sin(getAscension()))
);
res+=getDemiGrandAxe()*(sqrt(1.0-(getExcentricite()*getExcentricite())))
*sin(e)
*
(
(-sin(getArgument())*cos(getAscension()))
-
(cos(getArgument())*cos(getInclinaison())*sin(getAscension()))
);
return res;
}

long double Satellite::calculAnomalieExcentrique(const long double tps) const
{
long double m = calculAnomalieMoyenne(tps);
long double e=0;

long double tolerance = 0.0000000000001;//probleme ???
long double mMod = modulo(m,2*PI);
long double nbTours=m-mMod;

long double en=mMod;
long double enp=mMod+getExcentricite()*sin(en);
long double ecart=valAbsolue(enp-en);

while(ecart>tolerance)
{
en=enp;
enp=mMod+getExcentricite()*sin(en);
ecart=valAbsolue(enp-en);
}
return enp+nbTours;
}

long double Satellite::calculAnomalieMoyenne(const long double tps) const
{
return (calculMouvementMoyen()*(tps-getDateDebut()))+getAnomalie();
}


les fonctions retournent les bonnes valeurs, mais comment je vais faire alors pour mettre ça en oeuvre en temps réel ?

j'utilise en fait l'Open Scene Graph pour le rendu 3D.

merci de me répondre .