Matlab - Colebrook-White phenomenological equation

f = 0.0390

%% Datos

% Constantes

epsilon = 1.5e-4;       % [m]

Dh      = 0.1;          % [m]

Re      = 5000;         % [-]

% Ecuaciones como simbólicos

syms f;

eq1 = 1 / sqrt(f);

eq2 = -2 * log10(epsilon/(3.71*Dh) + 2.51/(Re*sqrt(f)));

%% Iteración

i       = 1;            % Inicialización del contador para el bucle

x       = 0.1;          % Valor de x inicial para empezar el bucle

error	= 0.2;          % Error absoluto inicial para entrar en el bucle

while error > 0.1

    % Calcular f

    f = 1 / x^2;

    % Evaluar parte izquierda

    eq1_sol = double (subs(eq1, f));

    % Evaluar parte derecha

    eq2_sol = double (subs(eq2, f));

    % Calcular error

    error = abs(eq1_sol - eq2_sol);

    % Aumentar x

    x = x + 0.1;

    % Aumentar contador

    i = i + 1;

end

%% Resultado

fprintf('f = %.3f (%d iteraciones)\n\n',f,i);

syms f positive

e=0.00015;  dh=0.10; re=5000;

g= sqrt(1/f)+2*log10(e/(3.71*dh)+2.51/(re*sqrt(f)));

fplot(g,[0.01 0.1])

raiz=fzero(@(f)sqrt(1/f)+2*log10(e/(3.71*dh)+2.51/(re*sqrt(f)))  ,[0.000001  100])

>> raizces_de_una_funcion_14_01_21_foro

raiz =

    0.0390