%Simpsons och trapetsregeln
b = 3.2;
h = 10;
N = b/h;
a = 0;
trapets = zeros();
x = a:h:b;
%V = pi*int(((y(x)).^2),[a b]);
for i = 1:N
    xk = a + (i - 1)*h;
    trapets(i) = (y(xk) + y(a + i*h))/2;
end

trapetssum = sum(trapets);
figure(1)
plot(0:N-1,trapets)


%simpsons

for i = 1:N
    x(i) = i*h;


end




function f = y(x)
    f = (7 + cos(x.^2 - 1))/(2 - cos(pi.*x));
end