使用Python（Scipy）查找常数可变的方程的分数分解

Question

我试图找到the following equation: (x^2+2x+a)/((x-1)^3(x^2+1))的分数分解，然后求解“ a”，其中其2到4的积分等于1。

我已经能够简单地使用Scipy.signal.residue进行分数运算，但是我不能在它们中使用诸如“ a”之类的变量（不能是多变量吗？）。如果有人知道我在做什么错，我很乐意接受帮助。

这是我到目前为止的代码：

import scipy.signal
import numpy as np
from sympy.abc import a
List = []
b = [1, 2, 6]  #numerator, q = random number
a = [1, -3, 4, -4, 3, -1]  #denominator coefficients
[r, p, k] = scipy.signal.residue(b, a)
Listr = r.tolist()
Listp = p.tolist()
Listk = k.tolist()

for i in range(len(Listr)):
    List.append(str(Listr[i])+"/"+str(Listp[i]))
List.append(str(Listk[0]))
res = "+".join(List)
print(res)

它用“ a”提供错误，但我不知道有任何解决方法。 感谢您的帮助，JW。