%% sciety stozek
% szukamy h(VS)i R(VS)
% R(VS)
RVS=((VS.*3)/(h*pi)+1).^(1/3);
% h(VS)
hVS=(3.*VS)./(pi.*(Rs.^2+(Rs.*r)+r.^2));