clc; clearvars; close all;

krok=0.01;

kvyskyrotoru=8^0.21;

rozsahrotor=(2.5:1:20.5)';
rozsah10m=rozsahrotor/kvyskyrotoru;
strednihodnoty=rozsahrotor(1:(end-1))+0.5;

prychlost=zeros(length(strednihodnoty),1);
for i=1:(length(rozsah10m)-1)
    prychlost(i)=pvetru(rozsah10m(i),rozsah10m(i+1),krok);
end

tabulka=[strednihodnoty prychlost]

a=sum(prychlost)

function Integral=pvetru(dolnimez,hornimez,krok)

a=6.51;
b=1.91;

x=dolnimez:krok:hornimez;

y=b.*x.^(b-1).*exp(-(x/a).^b)./a^b;

 plot(x,y);hold on;
 ylim([0 0.2])
grid on;

Integral=0;

for i=1:(length(x)-1)
    yi=(y(i)+y(i+1))/2;
    dS=krok*yi;
    Integral=Integral+dS;
end

end