def sh1(n, q):
    return sj(n, q) + 1j * sy(n, q)