我不时使用sympy，但不是很擅长。目前我坚持定义一个索引变量列表，即n1到nmax并对其进行求和。然后我希望能够进行导数:到目前为止，我尝试了以下方法:numSpecies=10n=IndexedBase('n')i=symbols("i",cls=Idx)nges=summation(n[i],[i,1,numSpecies])但是，如果我尝试对一个变量求导数，则会失败:diff(nges,n[5])我也尽量避免使用IndexedBase.numSpecies=10n=symbols('n0:%d'%numSpecies)k=symbols('k',integer=True)nto