function [res]=kurt(x) 
res=(sum(x.^4)/length(x))-3*(sum(x.^2)/length(x))^2; 
res=abs(res);