code: def SympyCreate(n):
      
      a = sympy.zeros(n)
      for i in range(1, n):
          a += sympy.sqrt(i) * basis_state(i, n) * basis_state(i-1, n).H
      return a
docs: Creation operator for a Hilbert space of dimension `n`, as an instance
      of `sympy.Matrix`