By applying the continuity equation at branch i, we loop over all connecting pipes and get

(1)

where is the flow rate in the pipe connecting junctions i and j.