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

function opti9(nFlows, nEdges, Lambda, qn, source, target, Mu)
    oFlows = reshape(Lambda, nEdges, nFlows)'
    oFlows[oFlows .< 1e-9] .= 0.0
    OFQM = zeros(qn, qn, nFlows)
    for k in 1:nFlows
        for x in 1:nEdges
            OFQM[source[x], target[x], k] = oFlows[k, x]
        end
    end
    
    qL = vec(sum(oFlows, dims=1)) ./ Mu
    qLM = zeros(qn, qn)
    for v in 1:nEdges
        qLM[source[v], target[v]] = qL[v]
    end
    
    return OFQM, qLM
end