%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Compute the sample (bias-corrected) excess kurtosis. 
% 
% CALL 
% 
% p = Kurt(x);
% 
% VERSION:     $Id: Kurt.m,v 1.1 2011/02/25 10:04:07 congedo Exp $
%
% HISTORY:     25-02-2011 G. Congedo
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function k = Kurt(x)

n = numel(x);
mu4 = sum((x-mean(x)).^4)/n;
sig = sqrt(sum((x-mean(x)).^2)/n);
k = mu4/(sig)^4;

% Bias correction
k = (n-1)/(n-2)/(n-3)*((n+1)*k-3*(n-1));

end