Fixed a small bug in MatrixOperations.cpp on line 38. b =1.0 not 0.0