Simulation made faster by analysing the matrix A first off. Now am about to make it faster by changing b each time rather than A!