velocity.m

来自「he basic idea of the method of bisection」· M 代码 · 共 10 行

 L=200/(3.28084);     
 D=6/(39.370);
 Rho=1000;              
 Mu=1e-3; 
 E=(4*(10^-4))/(3.28084);
 g=9.81;  
 C= E/D;
 Pdrop = 10600;
 
 y = incsearch2( @(v) (((-1.8 .* log10((6.9./(Rho.*v.*D./Mu)) + (C/3.7).^1.11)).^(-2)) .* L .* (v.^2).*(Rho .* g) ./ (D .* 2 .* g)) - (Pdrop ),0,23)