Numpy and Scipy: Using Sparse Matrices in our Monte Carlo Simulation (part 3)