clc, clear all, close all
format rat
%*** Problem aus Aufgabe 12 b) 
A = [1,1;1,2;1,2;1,5];
b = [1;2;1;1];

%*** berechne QR_Zerlegung mit Givens-Rotationen
B = qrGivens(A);
%*** berechen Q'*b
b = matVecQt(B,b);
%*** berechne x ( extrahiere R)
R=triu(B);
%*** Löse oberes LGS exakt
x=R(1:size(B,2),:)\b(1:size(B,2))
%*** berechne residuum
norm(R*x-b)
