def legendre(n, x): e = Integer(1) / (Integer(2)**n * fact(Integer(n))) * diff( (x**2 - 1)**n, x, n) return e.expand()
def legendre(n, x): e = Integer(1)/(Integer(2)**n * fact(Integer(n))) * diff((x**2-1)**n, x, n) return e.expand()