Page 38 - Tlahuizcalli CB-28_Neat

y3(i) = xi(11); 5. Conclusiones
vy3(i) = xi(12);
Hemos utilizado el lenguaje de programación
end SCILAB para abordar el Problema de los 3-Cuerpos
aproximándose a colisión total. En lugar de utilizar
endfunction
las librerías de SCILAB para la solución de la
ecuación diferencial asociada al problema,
Pasemos ahora a la ejecución de las rutinas.
Como se nota en el código, estamos usando = decidimos construir una rutina propia basada en el
1, 1 = 2 = 3 = 1 y el número de iteraciones del método de Runge-Kutta de orden 4. Lo anterior fue
método es = 20001. La rutina principal, runge2_3, hecho con la finalidad de mostrar como SCILAB
se ejecuta con ta = 0.0, tb = 10.0 y los valores para puede ser utilizado para practicar los aspectos de
ra y va se seleccionan para estar cercanos a una programación relacionados a la solución de
trayectoria de colisión total. Estos valores son ra = problemas de la física que tienen que ver con la
solución numérica de ecuaciones diferenciales.
[x1; y1; x2; y2; x3; y3] = [cos(0); sin(0); cos(2*%pi/3);
sin(2*%pi/3); cos(4*%pi/3); sin(4*%pi/3)] y va = [vx1; Las librerías gráficas de SCILAB brindaron una gran
vy1; vx2; vy2; vx3; vy3] = 0.75*[-cos(0); -sin(0); - ayuda en la búsqueda de las condiciones iniciales
cos(2*%pi/3); -sin(2*%pi/3) + eps; -cos(4*%pi/3); - adecuadas para el problema (ver la dirección de
internet http://www.scilab.org).
sin(4*%pi/3) – 2*eps], donde eps=0.1 (con eps = 0.0
tendríamos una condición cerca de colisión total).
Las desviaciones de la condición de rumbo de
colisión total relacionadas a eps, fueron
seleccionadas a prueba y error tras varios intentos
en los que se utilizó recurrentemente la librería para
la generación de gráficos de SCILAB. Los resultados
finales de trayectoria se muestran en la Figura 3.
Dr. Julio Saucedo Zul Julio Saucedo Zul
Como se observa en la figura 3, en nuestro evento Profesor investigador de tiempo completo. Universidad
cercano a colisión total, dos de las masas forman Autónoma de Coahuila. Facultad de Ciencias Físico
al final una binaria (estado ligado de dos masas) y Matemáticas. Línea de investigación: Sistemas Dinámicos y
Mecánica Celeste
la masa restante sale disparada en sentido
opuesto a la trayectoria de esta binaria.

Dr. Carlos Eduardo Rodríguez García
Profesor investigador de tiempo completo. Universidad
Autónoma de Coahuila. Facultad de Ciencias Físico
Matemáticas. Línea de investigación: Física Aplicada

M.C. José Luis Fraga Almanza
Figura 3. Gráfica de las trayectorias calculadas para el Profesor investigador de tiempo completo. Universidad
evento de casi colisión total de 3 masas. Autónoma de Coahuila. Facultad de Ciencias Físico
Matemáticas. Línea de investigación: Física Aplicada

37
Tlahuizcalli ISSN: 2448-7260 Año 10 Núm. 28 enero-abril 2024