% Harjutustund 12, ylesanne 2
f=@(x) x.^2.*cos(x.^4);
a=0;
b=10;
n=2000;
h=(b-a)/n;
% a) trapetsvalemiga
for i=1:(n+1)
 x(i)=a+(i-1)*h;
end
x
y=f(x)
disp('Solmed ja vastavad funktsiooni v22rtused:')
disp([x;y])

% integraal trapetsvalemiga
integraal=h/2*(y(1)+2*y(2)+2*y(3)+2*y(4)+y(5))
% Matlabi k2sk trapetsvalemi jaoks - funktsioon diskreetsel kujul
inttrapets=trapz(x,y)

% integraal Simpsoni valemiga
integraal=h/3*(y(1)+4*y(2)+2*y(3)+4*y(4)+y(5))
% Matlabi k2suga - funktsioon on valemi kujul, y=f(x)
integral(f,a,b)
[inttrapets integral(f,a,b)]