function X=mc_qr(A,B)

[Q,R,r,E]=qr(A); //le module décroît sur diag(R), r ="rang numérique adapté"
C=Q'*B;
X=E*(R(1:r,:)\C(1:r,:));

endfunction