El proceso de como se calcula la simulación de Monte Carlo en autómatas celulares se estudió en base a el código de NxLcau [4]. Este proceso se explica a continuación, así como su algoritmo.
Primero creamos una configuración inicial por medio de una simulación de un proceso Bernoulli de la siguiente manera:
Para cada célula 
para 
entero y M finito tenemos
donde q=1-p, x es un número aleatorio y p es la densidad que queremos en la configuración inicial, por lo que q varía de 0 a 1, y corresponde a uno y sólo un punto del eje de las X.
Por lo tanto dividiremos el eje de las X en N puntos y p estará dado por 
para 
.
Para cada configuración inicial aplicamos la regla del autómata hasta la R-ésima generación que deseamos conocer su densidad. Obtenemos la concentración de la R-ésima generación y la normalizamos dividiendo entre el tamaño de la red M de la forma 
y graficamos este valor en las cordenadas (p, y) actualizamos p y repetimos el proceso hasta terminar.