code: def headloss_manifold(FlowRate, Diam, Length, KMinor, Nu, PipeRough, NumOutlets):
      
      #Checking input validity - inputs not checked here are checked by
      #functions this function calls.
      ut.check_range([NumOutlets, ">0, int", 'Number of outlets'])
      return (headloss(FlowRate, Diam, Length, Nu, PipeRough, KMinor).magnitude
              * ((1/3 )
                 + (1 / (2*NumOutlets))
                 + (1 / (6*NumOutlets**2))
                 )
              )
docs: Return the total head loss through the manifold.