An algorithm is developed to solve the fundamental flow cases of fully-developed turbulent flow in a pipe and in a channel. The algorithm uses second-order finite-difference approximations for nonuniform grid spacing and is developed in such a way as to easily facilitate the implementation of several two-equation, Reynolds-Averaged-Navier-Stokes turbulence models. Results are included for the Wilcox 1998 k-ω model.