% F = COLEBROOK(R,K)快速、准确、稳健的计算
根据Colebrook方程，达西-维斯巴赫摩擦因数F:
% - - - - - - -
% 1 | k
% --------- = - 2 * Log_10 | ----- + ------------- |
%根号(F) | 3.7 R *根号(F) |
% - - - - - - -
%的输入:
% R:雷诺数(应该是>= 2300)。
% K:等效砂面粗糙度高度除以水力
%直径(默认K=0)。
％
%输出:
% F:摩擦系数。
％
%的格式:
% R, K和F是标量或兼容数组。
％
%的准确性:
%机器精度为所有R > 3和所有K>=0，
%，即物理兴趣的所有值。
％
%示例:F = colebrook([3e3,7e5,1e100]，0.01)

这个m文件比通过Lambert w函数或简单的近似来解Colebrook方程更有效。

PS:这个m-file第一次提交时，作者不小心上传了一个测试版本。这次重新提交应该是正确的。