Use pseudo inverse to match matlab/book behaviour

[mahdiolfat]

No description provided.

[mahdiolfat]

The mrdivide(A, B) matlab function (or simply A / B) serves multiple purposes, one of which is a least squares fit "if A is a rectangular m-by-n matrix with m ~= n, and B is a mtrix with n columns...", see https://www.mathworks.com/help/matlab/ref/mrdivide.html

The equivalent of this in numpy is achieved with numpy.linalg.lstsq(A.T, B.T).T