Mathematical models of the mamalian urine concentrating mechanism consist of a large system of coupled, nonlinear and stiff equations. An efficient numerical orthogonal collocation method was employed to solve the steady-state formulation of urine concentrating mechanism. This method was used to solve the stiff and high order equations of electrolyte transport in a central core, single nephron model of the renal outer medulla. The presented results were in good agreement with implicit finite difference method's results, but this new method was faster and more stable. Due to the greater stability and larger convergence domain of collocation method over Newton's method, a parametric study on concentrated urine was investigated. The results showed that this model was sensitive to non-ideal countercurrent exchange between medulla interstitium and vasa recta. Although distal tubule lies in the cortex interstitium, it affects on the inflow to the collecting duct. This study showed that the effect of changing membrane transport properties of distal tubule wall on properties of outflow from the outer medulla collecting duct was considerable.