function ku = kurtosis(xn,xnmean,xn1sig)
%
%  Name: kurtosis
%
%  Usage: ku = kurtosis(xn,xnmean,xn1sig)
%
%  Computes kurtosis (4th moment) of distribution
%  (approaches zero for normal distribution)
%

%
%  Version SWD970422
%

if exist('xnmean') ~=1, xnmean = mean(xn); end
if exist('xn1sig') ~=1, xn1sig = std(xn); end

ku = mean(((xn - xnmean)/xn1sig).^4)-3;

return
