function maxh = FOBJ_ksi(lambda)
%
global Plothandle H freq
%
meq=lambda(1);                  % massa equivalente do sistema de 1 GL
keq=lambda(2);                  % rigidez equivalente do sistema de 1 GL
ksieq=lambda(3);                % fator de perda do sistema de 1 GL
ceq=2*meq*(keq/meq)^0.5*ksieq;  % coeficiente de amortecimento viscoso equivalente
%
%
% Calculo da inertância - Modelo Viscoso ____________________________ 
for ii=1:length(freq)
  omeg=2*pi*freq(ii);
  Hnum(ii)=-omeg^2/(-omeg^2*meq+1i*omeg*ceq+keq);
end
  H12=20*log10(abs(Hnum)); 
% _______________________________________________________________________
%
set(Plothandle,'ydata',H12)
drawnow
%
erro=H-Hnum;
%
emq = conj(erro)*erro.';
%
% Função objetivo ______________________________________________________
maxh = emq;
