Simulation of many-spin system dynamics via sparse matrix methodology