clear all
clc
syms x
f=x^2
b=20
a=5
s(1)=0
n=20
delta_x=b-a/n
i=1
for x=a:b
t=subs(f,x)
s(i+1)=s(i)+t*delta_x
end