Numerical Treatment of The Liouville-Von Neumann Equation For Quantum Spin Dynamics