X=[0,2.5,3.1]
Y=[0,2.4,1]
x_pol=np.linspace(0,3.1)
pol_aprox=np.poly1d([-1,3.45,0])
pol_sci=interpolate.interp1d(X,Y,kind='quadratic')
plt.plot(X,Y,'bo')
plt.plot(x_pol,pol_aprox(x_pol),label='Polinomio aproximado en clase: -x^2+3.5x')
plt.plot(x_pol,pol_sci(x_pol),label='Polinomio con scipy')
plt.legend()
plt.show()