function eGlobal = berechneEGlobal(x)

wMess = xlsread('messwerte.xlsx','I2:I22')';
AMess = xlsread('messwerte.xlsx','J2:J22')';
tSim = 500;
ts = 0.001;
f = wMess/(2*pi);
A = 1;

% Testwerte
K    = x(1,1);        % K
zeta = x(1,2);      % zeta
w0   = x(1,3);        % w0

abweichung = zeros(1,length(wMess));
amplitude1 = zeros(1,length(wMess));
AdB1     = zeros(1,length(wMess));
phase1 = zeros(1,length(wMess));

for i = 1:length(wMess)
    omega     = wMess(i);
    options = simset('SrcWorkspace','current');
    sim('testEinlesenAusSimulink',[],options)
    simulationszeit= tout';
    eingangssignal = A*sin(wMess(i).*simulationszeit);
    ausgangssignal = yout';
    [amplitude1(i),phase1(i),~] = berechneAmpPhs(simulationszeit,f(i),eingangssignal,ausgangssignal);
    AdB1(i) = 20*log10(amplitude1(i)/A);
    abweichung(1,i) = AdB1(i)-AMess(i);
    disp(['Die Abweichung in dB bei i = ',num2str(i),' ist ',num2str(abweichung(1,i))])
end
eSquared1 = abweichung.^2;
eSquaredSum1 = sum(eSquared1);
eGlobal = eSquaredSum1;
disp(['eGlobal bei i = ',num2str(i),' ist ',num2str(eGlobal)])