Harmonic_Sandpile_Dynamics.gif ( 255 × 255 píxeles, tamaño de archivo: 661 KB, tipo MIME: imagen/gif , en bucle, 61 fotogramas, 30 s)
D(t)=(I-piso(t*dH))°. Por tanto, I corresponde a la identidad del grupo de pilas de arena en el dominio, dH es el laplaciano discreto de H, (.)° es un operador que mapea una configuración dada a su respectiva configuración recurrente en la misma clase de equivalencia (correspondiente al operador de relajación si dH es no positivo), y t es una variable real (mod 1) interpretada como el tiempo.
La animación muestra la dinámica de la pila de arena inducida por la función armónica H=x*y en un dominio cuadrado de 255x255 de Z². Por lo tanto, x e y corresponden a las coordenadas estándar de Z², con el origen asumido en el centro del dominio.