function Re_theta = Seyb(lambda,Tu)
    % Seyb 1972

    Re_theta = 1000./(1.2+0.7*Tu)+10*((0.09+lambda)./(0.0106+0.036*Tu));

end

