function mass_conservation_error = hydraulic_network_eval_with_minor_losses(H)
    global H1 num_edges edges f A r K_M;
    H = [H1; H];
    Q = zeros([num_edges, 1]);
    for e = 1 : num_edges
        from = edges(e, 1);
        to = edges(e, 2);

        Q(e) = sign(H(from)-H(to)) * ...
               compute_flow_with_minor_losses(H(from)-H(to), f(e), K_M(e));
    end

    mass_conservation_error = A * Q + r; % This will be 0 when we have a
                                         % good solution
end

Not enough input arguments.

Error in hydraulic_network_eval_with_minor_losses (line 3)
    H = [H1; H];