I have something like this:
coefs = [28, -36, 50, -22] print(numpy.roots(coefs))
Of course the result is:
[ 0.35770550+1.117926
I hope this help.
roots = np.roots(coeff); for i in range(len(roots)): if np.isreal(roots[i]): print(np.real(roots[i]))