function err = fitfun2(lambda)
%FITFUN2 Used by Myfit.
  
global err gain0 tau0 delay tauFreq threshold power

gain0 = lambda(1);
tau0 = lambda(2);
delay = lambda(3);
tauFreq = lambda(4);
threshold = lambda(5);

if (delay < 0) | (tau0 <= 0) | (tauFreq <= 0)
  err = Inf;
  return
end

sim('LPRealFilt', 190);

err = mean((expout-modelout).^2)/(mean(expout.^2)-mean(expout)^2)*100;
disp([transpose(lambda), err])