Оптимизация СМО

function opti4(qn, nFlows, fsource, ftarget, flows)
    fST = zeros(eltype(flows), qn, nFlows)
    for l in 1:nFlows
        fST[fsource[l], l] = -flows[l]
        fST[ftarget[l], l] = flows[l]
    end
    flowSourceTarget = zeros(eltype(flows), qn + 2, nFlows)
    flowSourceTarget[1:qn, :] = fST
    return flowSourceTarget
end