function Iapp=Eser2(f,a,b,N)
x=-cos(pi*[0:N]/N);
x=(a+b)/2+(b-a)/2*x;
fx=feval(f,x);
p=polyfit(x,fx,N);
q=polyint(p);
Iapp=polyval(q,b)-polyval(q,a);
