N=input('N. sottointervalli: ');
for k=1:2
    [t,u1,u2,u4]=prTaylor(k*N);
    y=t+exp(-t);
    e1(k)=max(abs(y-u1));
    e2(k)=max(abs(y-u2));
    e4(k)=max(abs(y-u4));
end
p1=log(e1(1)/e1(2))/log(2)
p2=log(e2(1)/e2(2))/log(2)
p4=log(e4(1)/e4(2))/log(2)