The above set of non-linear equations can be solved numerically in a number of ways using a variety of standard
numerical analysis software packages. However, as for large networks the number of equations and unknowns becomes
extremely difficult. Consequently, a special purpose equation solver has been developed, called QUEENSOD. This solver
fully optimizes the objective function of Equation 13. The formulation of Equation 13 makes a single approximation
namely the Stirling approximation, which has been shown to produce errors less than 1% for the range of values and
derivatives being typically considered in the problem. A sample application of the QUEENSOD software is presented later
in the paper.