% Harjutustund 12, ylesanne 1
% integraali leidmine ristkylikvalemi abil
f=@(x) sin(x)./x;
fplot(f,[0,pi])
a=pi/4;
b=pi/2;
h=pi/40;
% leiame osaloikude arvu
n=(b-a)/h % osaloike 10
% a) leiame integraali vasakpoolse ristkylikvalemi abil
% x=[x_0,x_1,...,x_11], x_0->x(1)
x(1)=pi/4
x(2)=pi/4+h
for i=1:11
 x(i)=a+(i-1)*h;
end
x
y=f(x)
disp('Solmed ja vastavad funktsiooni v22rtused:')
disp([x;y])

% nyyd integraal vasakpoolse ristkylikvalemi abil
pindala=0;
for i=1:10
  pindala=pindala+h*y(i);
end
integraal1=pindala

% b) leiame integraali parempoolse ristkylikvalemi abil
pindala=0;
for i=2:11
  pindala=pindala+h*y(i);
end
integraal2=pindala