Theta_acceleration_data/code/Fig_2/all_fit.m

function f=all_fit(aux)
autoc_step=.002;
autoc_bins=1:300; 
bins=(autoc_bins)*autoc_step*2; bins=bins(1:150);
    aux=mean(reshape(aux,[],150));
    aux=aux/max(aux);    
    c=do_fit(bins,aux);
    f=c.f;