Estudio y resolución numérica de las ecuaciones de evolución estelar con el método de Henyey

Las estrellas son cuerpos esféricos masivos cuya estabilidad surge del equilibrio entre la fuerza gravitacional debida a su masa, atractiva hacia el centro, y la presión. La gran mayoría de las estrellas mantienen este equilibrio durante gran parte de sus vidas, convirtiendo hidrógeno en helio y, posteriormente, elementos más pesados a través de procesos de fusión nuclear. Este camino evolutivo define el ciclo de vida de la estrella y moldea su estructura hasta el punto de colapso, en el momento que se rompe el equilibrio entre presión y fuerza gravitatoria. En este estudio, caracterizaremos la evolución estructural de estrellas en diferentes etapas mediante simulaciones numéricas basadas en el algoritmo de Henyey (1959). Para la elaboración de este algoritmo se emplearán los lenguajes de programación Python y Julia, así como técnicas de resolución de ecuaciones diferenciales de manera iterativa.