clear all exact_vlue=cos(0)-cos(1); for n=[17 33 65 129 257] % n 必須是4的倍數加1 h=1/(n-1); x=linspace(0,1,n); y_h=sin(x(2:2:n-1)); y_2h=sin(x(3:4:n-2)); midp_h = 2*h*sum(y_h); midp_2h = 4*h*sum(y_2h); p1=log2((exact_vlue-midp_2h)/(exact_vlue-midp_h)); fprintf('n = %3d , h= % 1.7f , p1 = %1.7f\n',n,h,p1) end