import scipy
import pylab

a = scipy.poly1d([3,4,5])
print a, ' is the polynomial'
print a*a, 'is its square'
print a.deriv(), ' is its derivative'
print a.integ(), ' is its integral'
print a(0.5), 'is its value at x = 0.5'

x = scipy.linspace(0,5,100)
b = a(x)
pylab.plot(x,b)
pylab.show()
