function x=mysolve(A,b)
%
% Uebungsblatt: 4
% Aufgabe: 4.1
%
% Matlab-Version: 7.1
%
% Eingabe-Parameter:
%   A  :   Regulaere quadratische Matrix deren LR-Zerlegung berechnet
%          werden soll
%   b  :   Rechte Seite des linearen Gleichungssystems
%
% Ausgabe-Parameter:
%   x  :  Loesung des linearen Gleichungssystems
%
% Unterprogramme:  mydecomp
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Berechnung der LR-Zerlegung von A
A=mydecomp(A);
% Bestimmung der Dimension der rechten Seite des Gleichungssystems
n=length(b);
% Anlegen des Loesungsvektors
x=zeros(n,1);
% Vorwaerts einsetzen
% x_1=b_1
x(1)=b(1);
% Schleife ueber x_2 ... x_n
for i=2:n
   % Berechnung von x_i
   x(i)=b(i)-A(i,1:i-1)*x(1:i-1);
end
% Rueckwaerts einsetzen
% Initialisierung von x_n
x(n)=x(n)/A(n,n);
% Schleife ueber x_n-1...x_1
for i=n-1:-1:1
    % Berechnung von x_i
    x(i)=(x(i)-A(i,i+1:n)*x(i+1:n))/A(i,i);
end