A=[ -1,1,1 ; 2,-2,2 ; 3,3,-3; 4,4,4 ]

disp("QRdec:"); [Q,R]=QRdec(A)

QR=Q*R, QTQ=Q'*Q

b=[1;2;3;4]
x=QRback(Q,R,b)
QRx=(Q*R)*x
