The depths for a steady, gradually varied flow in a rectangular open channel with infinite width obtained analytically by the Bresse formula and numerically using trapezoidal rule for integration have been compared. It has been shown that numerical integration may generate significant errors of flow values for looped river networks. The method for choosing the optimum integration step that minimizes numerical integration errors has been presented.