out pp;
on rounded;
% --- wyd wydajnosc pompy [m**3/s] ----
wyd:=1/3600;
% --- s1 przekroj rury doprowadz. [m**2] ----
d:=2.54/(2*100);
s1:=pi*d**2/4;
% --- v1 prdkosc w rurze doprowadz. [m/s] ----
v1:=wyd/s1;
% --- s1 przekroj rury za zwezka [m**2] ----
s2:=0.1*s1;

% --- roznica cisnien za i przed zwezka dp = p2 - p1 [Pa] ----
ro:=1000;
dp:=(1 - (s1/s2)**2 )*ro*v1**2/2;
% --- dpatm = p2 - p1 [atm] ----

dpatm:=dp/(1.033*(10**5));

shut pp;
bye;
end;