如何在python中从头开始编写Simplex方法(以tableau/矩阵形式)？

import numpy as np

def simplex_method(tableau):
    m, n = tableau.shape
    while np.min(tableau[-1, :-1]) < 0:
        col_index = np.argmin(tableau[-1, :-1])
        if np.all(tableau[:, col_index] <= 0):
            raise Exception("Linear programming problem is unbounded.")
        row_index = np.argmin(tableau[:-1, -1] / tableau[:-1, col_index])
        tableau[row_index, :] /= tableau[row_index, col_index]
        for i in range(m):
            if i == row_index:
                continue
            tableau[i, :] -= tableau[i, col_index] * tableau[row_index, :]
    return tableau[-1, -1]

# 创建初始表格
tableau = np.array([[1, -3, -1, 0, 0, 0],
                    [1, 1, 0, 1, 0, 4],
                    [2, -1, 0, 0, 1, 3]])

# 运行Simplex方法
optimal_solution = simplex_method(tableau)
print("Optimal solution: ", optimal_solution)