 This is discussed further on the Hardy Cross method site.The modern method is simply to create a set of conditions from the above Kirchhoff laws (junctions and head-loss criteria).Then, use a Root-finding algorithm to find "Q" values that satisfy all the equations.