clear;
format long;
a1 = 0;
b1 = 1;
x = [5,9,17,33];
y(1) = oula(a1, b1, x(1));
y(2) = oula(a1, b1, x(2));
y(3) = oula(a1, b1, x(3));
y(4) = oula(a1, b1, x(4));
for s = 1:4
z(s) = 0.9460831;
end
y(2,:) = z;
plot(x,y);

function result = oula(a, b, n)

h = (b-a)/n;
x = 0.0;

for m = 1:n
x0 = a + (m-1)*h;
x1 = a + m*h;
if(x0 == a)
y0 = 1;
else
y0 = (sin(x0))/x0;
end
y1 = (sin(x1))/x1;
x = x + y0 + y1;
end

result = x * (h/2);
end