Moore-penrose-inverse with Sage?
How to compute the moore-penrose inverse of a matrix?
What I have found is using Scipy (or Numpy):
M = matrix(2,2,(1,2,2,4))
from scipy import linalg
matrix(linalg.pinv(M))
Any other way in Sage?