def _eval_expand_func(self, **hints):
        n, m, theta, phi = self.args
        rv = (sqrt((2*n + 1)/(4*pi) * factorial(n + m)/factorial(n - m)) *
                exp(I*m*phi) * assoc_legendre(n, m, cos(theta)))
        # We can do this because of the range of theta
        return rv.subs(sqrt(-cos(theta)**2 + 1), sin(theta))
dataset/ETHPy150Open sympy/sympy/sympy/functions/special/spherical_harmonics.py/Ynm._eval_expand_func/swapped operands of `/`