function [F]=T8_A1markus(y,y_schaetz)

F=0;
[N,n]=size(y);
for i=1:N
    for j=1:n
        F=((y(i,j)-y_schaetz(j)).^2)+F;
    end
end
F=F/(n*N);
