凌乱的球谐函数给出了nan+nanj

Question

在我当前的代码中，我需要一个函数，我可以把（l，m，θ，φ）或者用scipy的符号（n，m，θ，φ）来表示，它将输出这些参数的球谐函数的值。scipy的sph_harm函数在尝试更高的l和m之前一直很好地工作。例如，如果我使用l=88，m=85，theta=0，phi=pi/2（即sph_harm（85,88,0，np.pi/2号文件))然后输出的值是nan+nanj。在

有什么办法避免这个问题吗？我需要能够投入任何合理的| 在