计算NDimensional数组sympy时出错

如何解决计算NDimensional数组sympy时出错

我在python中遇到此错误

import scipy as sc
import numpy as np
import sympy as sp
from sympy import *
from sympy import Matrix
from sympy.matrices import Matrix,eye,zeros,ones,diag,GramSchmidt,Trace
from sympy.interactive.printing import init_printing
from sympy import derive_by_array
from sympy import sin,cos 
import math
from scipy import integrate
from sympy.abc import i,j,k,l

theta,phi = symbols('theta,phi')
H=MutableDenseNDimArray([1 for i in range(81)],(3,3,3)) 
P=MutableDenseNDimArray([1 for i in range(81)],3))
for i in range(3):
    for j in range(3):
        for k in range(3):
          for l in range(3): 
              integrand = sp.lambdify([theta,phi],H[i,l]*sin(phi))
              P[i,l]=sc.integrate.dblquad (integrand,np.pi,2*np.pi) 

当我尝试运行此功能时，出现错误：引发NotImplementedError 我认为当代码计算P [i，j，k，l] = sc.integrate.d ...

时，此错误返回