% FBERN A function f needed in Bernstein's inequality % for Jacobi polynomials % function y=fbern(t,a,b,x) y=b-a-(a+b+1)*t+2*(1-t^2)*sum(1./(t-x));