Monte Carlo algorithms, such as the Metropolis-Hastings one, are commonly used to simulate Ising spin systems. This thesis explores an alternative approach by applying the replicator equation to simulate Ising models, replacing minimisation of the system's energy with the maximisation of the player's strategy. The analysis focuses on an Ising model for which spins vectors can take orientations corresponding to the principal axis of the systems. The replicator equation is modified to take into account interactions between neighbouring spins, both nearest and next-nearest neighbours, as well as external magnetic fields. Results show that both the Metropolis-Hastings and the proposed hybrid algorithm reach similar outcomes at a local level.