clear all exact_vlue=2/3; for n=[17 33 65 129 257 361 1025 4001 40001 400001] % n 必須是4的倍數加1 h=1/(n-1); x=linspace(0,1,n); y_h=(x(2:2:n-1)).^(1/2); y_2h=(x(3:4:n-2)).^(1/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 = %9d , h= % 1.7f , p1 = %1.7f\n',n,h,p1) end