norbert
05 de Noviembre del 2004

necesito un programa (código fuente) en C que resuelva la órbita de un planeta de masa "m" alrededor del sol, usar unidades astronómicas (UA)
Utilizando la tercera ley de Newton (Fg =m*a) planteo el sistema de dos ec. dif. ordinarias de segundo órden:
dVx/dt = G*M*x/r^3
dVy/dt = G*M*y/r^3 donde: r = (x^2+y^2)^(1/2)
Vx: velocidad orbital según abcisas
Vy: velocidad orbital según ordenadas
t: tiempo
G*M: igual a 1
r: distancia planeta-sol
Utilizar condiciones iniciales convenientes
Resolver solamente aplicando método de RUNGE-KUTTA de cuarto órden
gracias...