Removed the square constraint from lu decomposition

@bsoetaer 
The built-in matrix class still does lu decomposition, even on rectangular matrices.  See https://github.com/ruby/ruby/blob/trunk/lib/matrix/lup_decomposition.rb