function f = misfit(p)
global t w      
    a = p(1); b = p(2); c=p(3);
    f = (a*sech(b*(t-c)).^2 - w)'; 
end