function xstar = myAusgleich(A,b)
%
% Uebungsblatt: 7
% Aufgabe: 7.2
%
% Matlab-Version: 7.1
%
% Eingabe-Parameter:
%   A  :   Rechteckige Matrix
%   b  :   Rechte Seite
%
% Ausgabe-Parameter:
%   xstar  :  Naeherungsloesung nach der Methode der kleinsten Quadrate
%
% Unterprogramme:  keine
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Bestimmung der Dimension von A
[m,n]=size(A);
% Anlegen des Loesungsvektors
xstar=zeros(n,1);
% Bestimmung der QR-Zerlegung
[Q,R]=qr(A);
% b1=Q^T*b
b1=Q'*b;
% Rueckwaertseinsetzen
xstar(n)=b1(n)/R(n,n);
% Schleife ueber x_n-1...x_1
for i=n-1:-1:1
    % Berechnung von x_i
    xstar(i)=(b1(i)-R(i,i+1:n)*xstar(i+1:n))/R(i,i);
end