function X=mc_svd(A,B)

[U,S,V,r] = svd(A); //la diagonale de S décroît, r="rang numérique de A"
s=diag(1 ./diag(S(1:r,1:r)));
psinv=V(:,1:r)*s*U(:,1:r)'; //construction de pinv(A) vue en TD
X=psinv*B;

endfunction