Требуется составить блок-схему под данный код программы, а также добавление комментариев для каждой строки кода

import numpy
 
def simplex_method(A, b, c):
    m, n = A.shape
    B = numpy.arange(n - m, n)
    N = numpy.arange(n - m)
 
    while True:
        B_inv = numpy.linalg.inv(A[:, B])
        x_B = B_inv @ b
        x = numpy.zeros(n)
        x[B] = x_B
 
        c_B = c[B]
        c_N = c[N]
        r_N = c_N - A[:, N].T @ B_inv @ A[:, B] @ c_B
 
        if numpy.all(r_N >= 0):
            break
 
        j = N[numpy.argmin(r_N)]
 
        d = B_inv @ A[:, j]
 
        if numpy.all(d <= 0):
            return None
 
        x_Bd = x_B / d
        x_Bd[x_Bd < 0] = numpy.inf
        i = B[numpy.argmin(x_Bd)]
 
        B[B == i] = j
        N[N == j] = i
    return x
 
# Функция кода, находящая максимум функции, а также максимальное значение переменных системы линейных уравнений
def simplex_method_maximum(A, b, c):
    maximum = simplex_method(A, b, c)
    if maximum is None:
        print('Функция не ограничена!')
    else:
        print('Максимальное значение функции:', c @ maximum)
        print('Максимальное значение переменных системы линейных уравнений равны:', maximum)
 
A = numpy.array([[3, 5], [5, 2]])
b = numpy.array([15, 10])
c = numpy.array([5, 3])
 
simplex_method_maximum(A, b, c)