x=[0 1 2 3 4 5 6];
px = [0 0.1 0.2 0.4 0.2 0.1 0];
subplot(3,1,1)
bar(x,px)
mux = sum(x.*px)
sig2x = sum((x-mux).^2.*px)
sum((x-mux).^3.*px)
skewx = sum((x-mux).^3.*px)/sig2x.^1.5

y=[0 1 2 3 4 5 6];
py = [0 0.1 0.3 0.4 0 0.1 0.1];
subplot(3,1,2)
bar(y,py)
muy = sum(y.*py)
sig2y = sum((y-muy).^2.*py)
sum((y-muy).^3.*py)
skewy = sum((y-muy).^3.*py)/sig2y.^1.5


z=[0 1 2 3 4 5 6];
pz = [0.1 0.1 0 0.4 0.3 0.1 0];
subplot(3,1,3)
bar(z,pz)
muz = sum(z.*pz)
sig2z = sum((z-muz).^2.*pz)
sum((z-muz).^3.*pz)
skewz = sum((z-muz).^3.*pz)/sig2z.^1.5


