import scipy as S
import pylab as P
coefs = [  1. ,  -3.5,  -0.5,  11. , -17. ,   6. ]
pol = S.poly1d(coefs)
x = P.frange(-4,4,npts=400)
y = pol(x)
P.figure()
P.plot(x,y)
P.axhline(0)
P.ylim(-100,100)
print 'Coefficients of p(x):',coefs
print 'Roots of p(x):',pol.r
P.show()