Added fixed v's

Hey Riccardo, 

I have added a function assessing whether the algorithms works for a fixed v, varying A and b. 

It can either happen that

* np.linalg stops with `singular matrix` error or
* the algorithm does not completely solve the system. (1/500 or so)

Ben