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

function opti7(qn, nFlows, nEdges, QFlows, flowSourceTarget, Mu)
    Aeq = zeros((qn + 2) * nFlows, nEdges * nFlows)
    Beq = zeros((qn + 2) * nFlows)
    
    for j in 0:nFlows-1
        rows_range = (qn + 2)*j + 1:(qn + 2)*(j + 1)
        cols_range = nEdges*j + 1:nEdges*(j + 1)
        Aeq[rows_range, cols_range] = QFlows[:, :, j + 1]
        Beq[rows_range] = flowSourceTarget[:, j + 1]
    end
    
    Lambda0 = zeros(nFlows * nEdges)
    A = zeros(nEdges, nEdges * nFlows)
    
    for z in 1:nEdges, y in 0:nFlows-1
        A[z, z + y * nEdges] = 1
    end
    
    B = copy(Mu)
    
    return Aeq, Beq, Lambda0, A, B
end