Matrix simulation
by @oliverg