import numpy as np
import matplotlib.pylab as plt

def wartosc(f,x=0):
    return eval(f,{'x':x})

wielomian = 'x**2+x'
x=np.array([1,2,3])
y=np.array([1,3,8])
z=y+0
for i in range(0,len(x)):
    z[i]=wartosc(wielomian,x[i])

plt.plot(x,y, "ro", ms=10)
plt.plot(x,z,lw=2)
plt.xlim(0,4)
plt.show()