Hi, in this modified code, I added a prev_I variable to store the previous value of I. After computing I, we compute the ratio as I / prev_I and print it in the third column!
clear
syms f(x)
prev_I=1;
f(x) = exp(-((x).^2));
a =0; % lower limit
b = 1; %upper limit
fprintf('n\t\tIntegral\tRatio\n');
for j= 1:6
n = 2.^j;
h = (b-a)/n;
kamir = vpa(int(f,x,a,b));
s = 0.5*(f(a)+f(b));
for i=1:n-1
s = s+f(a+i*h);
end
I =vpa(h*s);
if i==1
ratio =1;
else
ratio = I / prev_I;
end
fprintf('%d\t\t%.4f\t\t%.4f\n', n, I, ratio)
prev_I = I;
end
