x = -1:0.05:1;
y = 1./(1+25*x.^2);
p=polyfit(x,y,9);
s=vpa(poly2sym(p))
f = polyval(p,x);
plot( x,f,x,y,'o ')