%LUv1
%[L,U]=lu(A)
%lu1=inv(U)*inv(L)*b

%LUv2
%[L,U]=lu(A)
%ylu2=inv(L)*b
%xlu2=inv(U)*ylu2

%QR
%[q,r] = qr(A)
%xgr=inv(r)*q'*B

%QR2
%[Q, R]=qr(A)
%xqr2=inv(R)*Q'*B

%SVD
%[U, S, V]=svd(A);
%xsvd=inv(V)*inv(S)*U'*B