function [qv, v]=vypocet_objemoveho_prutoku(p1,p2)
p_0=99991.5;%tlak barometrický
p_1=p_0+p1;%tlak pøed clonou
p_2=p_1-p2;%tlak za clonou
p=p_1-p_2;%rozdíl tlaku pøed a za clonou
D=0.13;%vnitøní prùmìr potrubí
d=0.095;%vnitøní prùmìr clony
B=d/D;%pomìr prùmìrù
K=1.4;%izoentropický exponent vzduchu
ny=15.32*(10^-6);%kinematická viskozita pro 20°C
ro=1.159;%hustota vzduchu v laboratoøi
L_1=1;%pomìr vzdálenost pøedního odbìru tlaku od pøední strany clony
L_2=0.47;%pomìr vzdálenost zadního odbìru tlaku od zadní strany clony
S=(pi*(D^2))/4;%prùtoèná plocha
M_2=(2*L_2)/(1-B);%koeficient M2

v=50;%volím
v_1=40;%volím
while abs(v_1-v)>0.001
v_1=v;
Re=(v*D)/ny;%reynoldsovo èíslo
E=1-(0.351+0.256*(B^4)+0.93*(B^8))*(1-(p_2/p_1)^(1/K));%souèinitel expanze E
A=((19000*B)/Re)^0.8;%koeficient A
C=0.5961+0.0261*(B^2)-0.216*(B^8)+0.000521*(((10^6)*B)/Re)^0.7+(0.0188+0.0063*A)*(B^3.5)*((10^6)/Re)^0.3+(0.043+0.08*exp(-10*L_1)-0.123*exp(-7*L_1))*(1-0.11*A)*((B^4)/(1-B^4))-0.031*(M_2-0.8*(M_2^1.1))*B^1.3;%výpoèet souèinitele prùtoku C
qm=(C/sqrt(1-B^4))*E*(pi/4)*(D^2)*sqrt(2*p*ro);%výpoèet hmotnostniho prùtoku
qv=qm/ro;%výpoèet objemového prùtoku
v=qv/S; %výpoèet rychlosti
end





