function S=NCotesJF(f1,a,b,n,e)
f=inline(f1);
syms t
for k=0:n;
    S1=1;
    S2=1;
    S3=1;
    for i=k:-1:1
       S2=S2*i;
    end
    for i=n-k:-1:1
       S3=S3*i;
    end
    for i=0:n
        if i~=k
            S1=(t-i)*S1;
        end
    end
    f2=inline(S1);
    C(k+1)=(-1)^(n-k)/(n*S2*S3)*quad(f2,0,n);
end
n1=1;
h=(b-a);
h1=h/n;
S=0;
for i=0:n
    S=S+C(i+1)*f(a+i*h1);
end
S=(b-a)*S;
n1=2;
h=(b-a)/n1;
h1=h/n;
R=0;
for j=1:2
    for i=0:n
        R=R+C(i+1)*f(a+h*(j-1)+i*h1);
    end
end
R=(b-a)*(R+C(1)*f(a)+C(n+1)*f(b))/2;
S
R
%while abs(R-S)>e
    %S=R;
    %S
    %n1=2*n1;
    %h=(b-a)/n1;
    %h1=h/n;
    %R=0;
    %for j=1:n1
        %for i=0:n
            %R=R+C(i+1)*f(a+h*(j-1)+i*h1);
            %end
            %end
    %R=(b-a)*R;
    %end