for (k=0 ; k<9 ; k++)
{
mT = mT + M[k];
}for (int k=0 ; k<9 ; k++)
{
mT = mT + M[k];
}for (int k=0 ; k<9 ; k++)
{
mT = mT + M[k];
cout<< setprecision (12)<< mT <<endl;
}double M[8];
M[0]= 1.0;
M[1]= 0.0000001747;
M[2]=0.0000025759;
M[3]= 0.0000031606;
M[4]= 0.00000033960;
M[5]=0.00095455;
M[6]= 0.00030077;
M[7]= 0.000045931;
M[8]= 0.0000541958;
en faite, je vois tout de suite qu'elle est incorrecte car j'ai convertis en masse solaire, et mon résultat final est inférieur à 1 ..
int N = 9; // Soleil + 8 planetes
double a = 0;
double Masse[N]; // Masse de N (les masses sont exprimees en masse solaire)
Masse[0] = 1.0;
Masse[1] = 0.00000016602;
Masse[2] = 0.0000024478;
Masse[3] = 0.0000030035; //on optimise la precision en prenant 5 chiffres significatifs
Masse[4] = 0.00000032272;
Masse[5] = 0.00095460;
Masse[6] = 0.00028582;
Masse[7] = 0.000043658;
Masse[8] = 0.000051501;
int i; // Variable pour balayer toutes les planetes avec une boucle for
for (int i = 0 ; i<N ; i++)
[INDENT]{
[INDENT]a= a + Masse [i];
cout<< setprecision (12)<<" "<< a <<endl;[/INDENT]
}[/INDENT]
// On va maintenant définir ainsi toutes les coordonnées et velocites de chaque planetes sous forme de tableau
double xG,yG,zG,vxG,vyG,vzG,mT; // Coordonnees du centres de masse et masse totale
xG = 0.0;
yG = 0.0;
zG = 0.0;
vxG = 0.0; // On doit prendre les "coordonnees" du centre de masse initiale
vyG = 0.0; // où on ajoutera ensuite ceux des autres planetes
vzG = 0.0;
mT = 0.0;
// Coordonnees, impulsions et masses totales
int k; // Variable pour balayer toutes les planetes avec une boucle for
for (int k = 0 ; k<N ; k++)
[INDENT] {
[INDENT] mT= mT + Masse[k];
cout<< setprecision (12)<<" "<< mT <<endl;[/INDENT]
}[/INDENT]
cout<< setprecision (12)<<" "<< mT <<endl;
le plus simplifié
bah je met quasi tout le code alors, car c'est un gros tableau apparemment qui obstrue ma valeur
cout<< setprecision(12) << xG << " " << yG << " "<< zG << " " << vxG << " " << vyG << " "<< vzG << " " << mT <<endl;
{
[INDENT] ofstream f("1.txt");
f <<a[0][0]<<" "<< a[1][0]<<" " << a[2][0]<< " " << a[3][0] << " "<< a[4][0]<< " "<< a[5][0];
ofstream g("2.txt");
g <<a[0][1]<<" "<< a[1][1]<<" " << a[2][1]<< " " << a[3][1] << " "<< a[4][1]<< " "<< a[5][1];
......
[/INDENT]
}
f <<a[0][0]<<" "<< a[1][0]<<" " << std::flush;f <<a[0][0]<<" "<< a[1][0]<<" " << std::endl;ofstream f("1.txt");
f << 1 <<" "<< 2<<" " << 3<< " " << 4 << " "<< 5 << " "<< 6 << std::endl;
void ecriture(double t, double a[6][N])
{
[INDENT]ofstream f("1.txt");
f << 1 <<" "<< a[1][0]<<" " << a[2][0]<< " " << a[3][0] << " "<< a[4][0]<< " "<< a[5][0] << std::endl;[/INDENT]
}
Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité
Tu pars déja ?
Identification
Pas encore inscrit ?
Ou identifiez-vous :