% Dieses Programm stellt die Normalengleichungen fuer das
% Zerfallsproblem auf, und loest diese.

clear t y
load daten

g = [exp(-t/10);
     exp(-t/20);
     exp(-t/5)];

A = zeros(3,3);
b = zeros(3,1);
for i=1:3
  for j=1:3
    A(i,j) = sum(g(i,:) .* g(j,:));
  end
  b(i) = sum(y .* g(i,:));
end

% Die Normalenmatrix ist schlecht konditioniert
% (dies ist typisch fuer diese Art von Problemen):
cond(A)

% Berechne den Loesungsvektor.
% Er wird deutlich von der Loesung des ungestoerten Problems
% [7; 2; 30] abweichen.
A\b

