# 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?