%% MNČ s QR rozkladem
%

A=randi([-10,10],10,3)
b=[1:10]';

%

x1=A\b
x2=MNC_QR(A,b)
norm(x1-x2)

%
%% Pseudoinverze s pomocí QR rozkladu

[Q1,R1]=qr(A,0);
X1=inv(R1)*Q1'
X2=pinv(A)
norm(X1-X2)

%

function x=MNC_QR(A,b)
% function x=MNC_QR(A,b)
% Funkce hledá řešení ve smyslu NČ pomocí QR rozkladu

[Q1,R1]=qr(A,0);
c=Q1'*b;
% c=c(1:length(R1));
x=inv(R1)*c;

end