def askew(a,dimension=None):
    denom = N.power(amoment(a,2,dimension),1.5)
    zero = N.equal(denom,0)
    if type(denom) == N.ndarray and asum(zero) <> 0:
        print "Number of zeros in askew: ",asum(zero)
    denom = denom + zero  # prevent divide-by-zero
    return N.where(zero, 0, amoment(a,3,dimension)/denom)
Returns the skewness of a distribution (normal ==> 0.0; >0 means extra
dataset/ETHPy150Open radlab/sparrow/src/main/python/third_party/stats.py/askew