Un método para generar modelos de galaxias elípticas consiste en partir de un sistema de partículas con una configuración adecuada y seguir su evolución con un código de N-cuerpos hasta que alcanza el equilibrio. En trabajos anteriores hemos utilizado el código de L.A. Aguilar para generar modelos tanto no cuspidales como cuspidales, ofreciendo los segundos una mejor representación de las galaxias elípticas. Sin embargo, el método propuesto por Hernquist & Ostriker es más adecuado para los casos con cúspides centrales, puesto que utiliza una expansión del potencial en funciones radiales que ya tiene caracter cuspidal. Aquí consideramos dos modelos, uno no cuspidal y otro cuspidal, obtenidos con el método de Aguilar y cuyos potenciales ajustamos ahora con la expansión de Hernquist & Ostriker. Obtuvimos las ecuaciones variacionales correspondientes y las empleamos para determinar los exponentes de Lyapunov de ´orbitas ya estudiadas con el método de Aguilar y separar las caóticas. Las órbitas regulares las clasificamos mediante el análisis de sus frecuencias orbitales. Presentamos una comparación de los resultados obtenidos con los dos métodos para ambos tipos de órbitas.

One way to create models of elliptical galaxies is to start from a system of particles with a suitable distribution and to follow its evolution with an N-body code until it reaches equilibrium. In previous works we used the code of L.A. Aguilar to generate noncuspy and cuspy models, the latter giving a better representation of elliptical galaxies. However, the method proposed by Hernquist and Ostriker is more suitable for the cases with central cusps, because it uses an expansion of the potential in radial functions that already reflects the cuspy character. Here we consider two models, one non-cuspy and another cuspy, obtained with the method of Aguilar and whose potentials we now fit with the Hernquist and Ostriker expansion. We obtained the corresponding variational equations and we used them to determine the Lyapunov exponents of orbits already investigated with the method of Aguilar and to separate the chaotic ones. We classified the regular orbits through the analysis of their orbital frequencies. Here we present a comparison of the results obtained with the two methods for both types of orbits.

En este artículo se presentan estrategias didácticas en tono lúdico para explorar la metodología de Monte Carlo y se analizan las puestas en práctica realizadas en diversos escenarios educativos. Las propuestas elaboradas permiten poner en acción conceptos de probabilidad y estadística de un modo divertido y atrayente mediante la realización de “experimentos-juegos” cooperativos con objetos reales y virtuales. Cada estrategia está dirigida a profesores de escuelas secundarias y de institutos de formación docente, y también al ciclo básico de la Universidad, para que sea modificada, recreada y profundizada de acuerdo a las áreas específicas de aplicación.

This paper presents didactic strategies in a ludic tone to explore the methodology of Monte Carlo and analizes the implementations that have been performed in several educational scenarios. These proposals allow us to put into action concepts of probability and statistics in an amusing and attractive way by conducting cooperative “experiments-games” that are based on the use of real and virtual objects. Each strategy is aimed to teachers of high school and teacher training institutes, and also for early stages at University level, in order to be modified, recreated, and enhanced according to the specific areas of application.