在Python中使用NumPy对Legendre级数进行积分的完整攻略：

1. 引入必要的库

import numpy as np
from scipy.integrate import fixed_quad
from numpy.polynomial.legendre import leggauss

2. 定义需要计算的Legendre级数

def legendre_series(x, n):
series_sum = 0
for i in range(n + 1):
series_sum += ((-1)i * np.math.factorial(2 * n – 2 * i)) / \
(2 n * np.math.factorial(i) * np.math.factorial(n – i) * np.math.factorial(n – 2 * i)) * \
x ** (n – 2 * i)
return series_sum

3. 定义需要积分的函数

def legendre_integral_func(x, n):
return legendre_series(x, n) ** 2

4. 设置积分下限

a = -1

5. 使用fixed_quad函数对函数进行积分

result, _ = fixed_quad(legendre_integral_func, a, 1, args=(2, ), n=5)

6. 输出积分结果

print(‘The integral of Legendre series squared from {} to 1 is {:.8f}’.format(a, result))

在以上代码中，我们首先引入了必要的库，分别是NumPy库、Scipy库的fixed_quad函数，以及Numpy库的Legendre-Gauss积分点的计算函数。然后，我们定义了需要计算的Legendre级数和需要积分的函数。在积分之前，我们需要先设置积分的下限，代码中我们选取了下限为-1。最后使用fixed_quad函数进行积分，并输出积分结果。

下面我们再用另一个例子说明如何设置别的积分下限，比如0：

定义需要积分的函数

def f(x):
return x**2

设置积分下限

a = 0

使用fixed_quad函数对函数进行积分

result, _ = fixed_quad(f, a, 1, n=5)

输出积分结果

print(‘The integral of x^2 from {} to 1 is {:.8f}’.format(a, result))

在以上代码中，我们定义了一个需要积分的函数，即f(x) = x^2。然后我们设置了积分的下限为0，并使用fixed_quad函数进行积分。最后输出积分结果。