Este codigo fue realizado en un cuaderno de Jupyter Notebook utlizando el lenguaje de programacion de Python. Por lo tanto, para poder utilizar este codigo debe ser ejecutado usando un jupyter notebook.
Para el uso de este codigo se debe ejecutar las celdas que definen todas las funciones en el algorimo antes de intentar ejecutarlas.
Las funciones de donde se obtuvieron las graficas fueron siguientes:
*graficaOrbital(1/24,0,0.22,0.08,1000000,1000) para poder graficar las trayectorias de un sistema con una energia de 1/24, x=0, y0=0.22, vy0 =0.08 con 10000000 puntos y un tiempo de 1000.
*PoincareSection(1/8,0,0.2,0.2,1000000,10000,taj,3,3) grafica la seccion de poincaré correspondiente al plano x=taj para un sistema con energia de 1/8, x=0, y0=0.2, vy0 =0.2 con 10000000 puntos y un tiempo de 10000.
*Exponentes2([1/24,0,0.128846,0.018450],[0.01,0],100000) para calcular los exponentes de lyapnuv para un sistema con los siguientes condiciones iniciales E=1/24,x=0,y=0.128846 y vy=0.018450, donde la diferencia entre los puntos iniciales estarán dados por delta x=0.01 y delta y =0 par aun total de 100000
Si hay alguna duda sobre el codigo se puede comunicar a los autores: Brayan Garcia Vargas- [email protected] Miguel Angel Gonzales- [email protected]