function fatorDeAtrito
  input Real Re;
  input Real epsilon;
  output Real f;
protected
  parameter Real N = 10000;
  parameter Real erro = 0.01;
  Real f0 = 0.02;
  Real s = 1/sqrt(f0);
  Real s0 = 0;
  Real i = 1;
algorithm
  if Re > 1000 then
    while (i<N and abs(s-s0)>erro) loop
      s0 := s;
      s := -2*log10(epsilon/3.7 + 2.51*s0/Re);
      i := i + 1;  
    end while;
    f := 1/s^2;
  elseif Re > 100 then
    f := 64/Re;
  else
    f := 0.60;
  end if;
end fatorDeAtrito;