clear ll
clc
a=1;
b=4;
n=20;
x=1;
A(1)=0
delt_x=b-a/n
for i=1:n
x(i)=(a+2*i-1)*delt_x/2
f=x(i)^2+3
A(i+1)=A(i)+f*delt_x
end
plot(x,i)