function f=quadfehler(y,y_s)

% f=quadfehler(y,y_s) (T8A2_Ute)
% Bekommt als Eingaben eine Matrix y mit Messwerten zu den Parameterwerten,
% (wobei die Messwerte zum gleichen Parameterwert jeweils in einer Zeile
% stehen) sowie einen Spaltenvektor y_s mit einem Schätzwert für jeden
% Parameterwert;
% gibt zurück, wie groß der mittlere quadratische Abweichung ist.

[N,n]=size(y);

y_s_mat=repmat(y_s,1,n);
r=(y-y_s_mat).^2;
f=sum(r(:))/(N*n);