exercise10a.m

% Loesung von Aufgabe 10, Matrixversion
% funktioniert auch fuer beliebige Zahl von Messwerten

% gegebene Werte
alpha = [-64.5 62.0 147.6 -122.7]'*pi/180;  % Spaltenvektor
punkte = [2.7 2.5; 36.1 5.5; 28.4 31.9; 4.0 24.6];

% Berechnung der Steigungen
m = tan(0.5*pi + alpha);

% Geraden sind gegeben durch
% -m*x + y = b = -m*x0 + y0
% also
n = size(alpha, 1);
A = [-m, ones(n,1)];
r = -m.*punkte(:,1) + punkte(:,2);

% Damit lautet die Loesung
x = A \ r