MATLAB代码以计算雷诺数（RE）和相对粗糙度系数（Epsilon）的给定值的管道中的摩擦因子。

Syntax:
F = Colebrook（RE，Epsilon）

Example 1: Single Re, single epsilon
re = 1e5;
epsilon = 1e-4;
F = Colebrook（RE，Epsilon）

示例2：多个RE，单个Epsilon
RE = 5000：1000：100000;
epsilon = 1e-4;
F = Colebrook（RE，Epsilon）;
情节（RE，F）

Example 3: Single Re, multiple epsilon
re = 1e5;
epsilon = linspace（1e-4,1e-1,100）;
F = Colebrook（RE，Epsilon）;
情节（Epsilon，F）

示例4：多个RE，多个Epsilon
Re = logspace(4,8,100);
epsilon = linspace（1e-4,1e-1,100）;
[re，epsilon] = meshgrid（re，epsilon）;
F = Colebrook（RE，Epsilon）;
冲浪（RE，Epsilon，F）

参考：
[1] Colebrook，C。F.和White，C。M.（1937）。液体实验
粗糙的管道中的摩擦。皇家学会论文集
伦敦。A系列 - 数学和物理科学，161（906），
367-381.
[2] Colebrook，C。（1939）。管道中的湍流，特别是
参考光滑和粗糙之间的过渡区域
管道法。土木工程师学会杂志，11（4），
133-156。