% Uso chebpol_int y chebpol_clenshaw
 f = @(x)1./(1+25*x.^2);
 [c,x] = chebpol_int(f,11);
 t = -1:.01:1;
 plot(t,f(t),'b',t,chebpol_clenshaw(c,t),'r',x,f(x),'ok')
