def SympyCreate(n):
    """Creation operator for a Hilbert space of dimension `n`, as an instance
    of `sympy.Matrix`"""
    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
Creation operator for a Hilbert space of dimension `n`, as an instance
    of `sympy.Matrix`