function xi = xi2D(ff,p)
% Eqn (13) in the paper for nonlinear fitting. You don't need to touch this function.
A = p(1); B = p(2);

xi = A./sqrt(ff).*exp(B./sqrt(ff));
