Friction Factor Colebrook Equation (MATLAB) using fzero syntax

#matlab #colebrook #frictionfactor #fluidmechanics
clc; clear all; close all;
%initial guess
f0=0.1;
[f, yval]=fzero(@cb,f0)
%Colebrook Equation
function y=cb(f)
NRe=6.5e4; %Reynold Number
eD = 0.00013; % e/D
y=1/sqrt(f)+0.86*log(eD/3.7+2.51/NRe/sqrt(f));
end