a=1
b=5
n=50
h=(b-a)/n
x=a:h:b
y=sin(x).^2
s=1
for k=1:n
s=((y(k)+y(k+1))/2)*h

end
syms x
int(sin(x).^2)
w=x/2 - sin(2*x)/4
z=b/2 - sin(2*b)/4
u=a/2 - sin(2*x)/4
v=z-u

o1=6.716720366386277e-001, n1=5
o2=3.813915150480616e-001, n2=10
o3=7.508196606063232e-002, n3=50
o4=3.718934741559301e-002, n4=100
j=[o1 o2 o3 o4]
l=[n1 n2 n3 n4]
plot(l,j)