% Emacs mode header: -*-Octave-*- source. %
% G for Greenwood tuning function
% G lat analyt

function ansf = Glat_analyt(x,w)
d = 7; %[mm]
w0 = 20000; %[Hz]
% out = w./sqrt(w_0.*exp(- x./d) - w); %Eq.[11]

ansf = (2*d*atanh((w0./exp(x/d) - w).^(1/2)./(-w).^(1/2)))./(-w).^(1/2);
endfunction
