clear all

exact_value=cos(0)-cos(1);
for n=[32 64 128]
    h=1/n;
    xg_h=0:h:1;
    y_h=sin(xg_h);
    tra_h = (h/2)*(y_h(1) + 2*sum(y_h(2:1:end-1)) + y_h(end));
    
    hh=2*h;
    xg_hh=0:hh:1;
    y_2h=sin(xg_hh);
    tra_2h = (h)*(y_2h(1) + 2*sum(y_2h(2:1:end-1)) + y_2h(end));
    
    p1=log2((exact_value-tra_2h)/(exact_value-tra_h));
    fprintf('n = %3d , h= % 1.7f , p1 = %1.7f\n',n,h,p1)
end