5 lines
167 B
Text
5 lines
167 B
Text
|
pseudoinv(M)=my(sz=matsize(M),T=conj(M))~;if(sz[1]<sz[2],T/(M*T),(T*M)^-1*T)
|
||
|
|
addhelp(pseudoinv, "pseudoinv(M): Moore pseudoinverse of the matrix M.");
|
||
|
|
|
||
|
|
y*pseudoinv(X)
|