code: def skew(self):
          "skewness"
          n, a, b = self.n, self.a, self.b
          t1 = (a+b+2*n) * (b - a) / (a+b+2)
          t2 = sqrt((1+a+b) / (n*a*b * (n+a+b)))
          return t1 * tf skew(self):
          "skewness"
          n, a, b = self.n, self.a, self.b
          t1 = (a+b+2*n) * (b - a) / (a+b+2)
          t2 = sqrt((1+a+b) / (n*a*b * (n+a+b)))
          return t1 * t2
docs: skewness