Theta_acceleration_data/code/Fig_2/do_fit.m

function [c,g]=do_fit(bins,aux)

sel=bins>0.015 & bins<0.3;


fo = fitoptions('Method','NonlinearLeastSquares',...
               'Lower',[0.2,0,0,0,6,-1,-1],...
               'Upper',[2,20,1,Inf,12,1,1],...
               'StartPoint',[1,3,.2,10,8,0,0]);
% ft = fittype('a-b*x+c*cos(2*pi*f*x)*exp(-d*x)','options',fo);
ft = fittype('a*exp(-b*x)+c*cos(2*pi*f*x)*exp(-d*x)+g+h*x','options',fo);


[c, g] = fit(bins(sel)',aux(sel)',ft);