PROBLEMA DE LOS DOS CUERPOS

[cypascal]

Les dejo en el siguiente enlace el programa:
http://cypascal.blogspot.com.es/2012/11/orbita-de-la-tierra-entorno-al-sol.html

y aqui mismo:

Código [Seleccionar] Expandir

#include <stdio.h>
#include <math.h>


const int L=31536000;/*SEGS. EN UN AÑO*/
const double dt=1;/*DT ES UN SEG.*/
const double m1=5.97e24, m2=1.99e30, G= 6.674e-11;

int i;
double F,coseno,seno,d;

struct estado{
double X;
double Y;
double Vx;
double Vy;
double Ax;
double Ay;
}Solido1[L];

//Situacion inicial;


int main()
{
Solido1[0].X=1.5e11;
Solido1[0].Y=0;
Solido1[0].Vx=0;
Solido1[0].Vy=2.98e4;

d=sqrt(pow(Solido1[0].X,2)+pow(Solido1[0].X,2));
F=(G*m1*m2)/(pow(d,2));

Solido1[0].Ax=-F/m1;
Solido1[0].Ay=0;

printf("\n\n     ORBITA DE LA TIERRA EN TORNO AL SOL\n\n\n");

for (i=1;i<=L;i++)
{
Solido1[i].X=Solido1[i-1].X+Solido1[i-1].Vx*dt;
Solido1[i].Y=Solido1[i-1].Y+Solido1[i-1].Vy*dt;

d=sqrt(pow(Solido1[i].X,2)+pow(Solido1[i].Y,2));
coseno=Solido1[i].X/d;
seno=Solido1[i].Y/d;

Solido1[i].Vx=Solido1[i-1].Vx+Solido1[i-1].Ax*dt;
Solido1[i].Vy=Solido1[i-1].Vy+Solido1[i-1].Ay*dt;

F=(G*m1*m2)/(pow(d,2));
Solido1[i].Ax=-F*coseno/m1;
Solido1[i].Ay=-F*seno/m1;

if ((i%86400)==1){/*PARA NO IMPRIMIR DEMASIDADOS DATOS*/
printf(" DIA %4d)  X=%10.4f MKm",i/86400,Solido1[i].X/1000000000);
printf("   Y=%10.4f MKm",Solido1[i].Y/1000000000);
printf("   d=%10.4f MKm\n",d/1000000000);
}
}

printf("\n\n MKm significa MegaKilometro o Millones de kilometros.");
return 0;
}



Un Saludo!

[saraou]

Muchisimas gracias! ahora ya puedo avanzar con el trabajo.
Un saludo!