function [s,mu,sig]=skewness(x)
mu=mean(x);
sig=sqrt(var(x));
s=(mean((x-mu).^3))/sig^3;
