global err  gain0 tau0 delay tauFreq threshold power 

lam = transpose([gain0 tau0 delay tauFreq threshold]);

options = optimset('Display','off','TolX',0.1);
lambda = fminsearch('fitfun2',lam,options);

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

sim('LPRealFilt', 190);