clear clc delete demo25117c.log diary demo25117c.log format long n = 4083 x = 2; target = log(x) trapezoidSum = (x^2-1)/(2*n*x) + sum( (x-1)./(n+[1:n-1]*(x-1))) error = abs( trapezoidSum - target ) if error < 10^(-8) disp('Yes, error is less than 10^-8'); end diary off