Формирование симметрической положительно определенной матрицы

import random
import time
 
n = 5  # размерность матрицы
 
# Замеряем время начала операции
start_time = time.time()
 
# Генерируем матрицу
matrix = [[random.randint(0, 9) for _ in range(n)] for _ in range(n)]
 
# Делаем матрицу симметричной
for i in range(n):
    for j in range(i, n):
        matrix[i][j] = matrix[j][i]
 
# Проверяем, что матрица положительно определена
is_positive_definite = True
for k in range(1, n+1):
    submatrix = [row[:k] for row in matrix[:k]]
    determinant = 1
    for i in range(k):
        determinant *= submatrix[i][i]
        for j in range(i+1, k):
            ratio = submatrix[j][i] / submatrix[i][i]
            for l in range(i, k):
                submatrix[j][l] -= ratio * submatrix[i][l]
    if determinant <= 0:
        is_positive_definite = False
        break
 
# Замеряем время окончания операции
end_time = time.time()
 
if is_positive_definite:
    print("Сгенерирована положительно определенная матрица:")
    for row in matrix:
        print(row)
else:
    print("Сгенерирована не положительно определенная матрица.")
 
# Выводим время, затраченное на выполнение операции
print("Время выполнения операции:", end_time - start_time, "секунд.")

import random
import time
import matplotlib.pyplot as plt
 
n = 5  # размерность матрицы
 
# Список для хранения временных затрат
times = []
 
for i in range(10):  # выполняем операцию 10 раз
    # Замеряем время начала операции
    start_time = time.time()
 
    # Генерируем матрицу
    matrix = [[random.randint(0, 9) for _ in range(n)] for _ in range(n)]
 
    # Делаем матрицу симметричной
    for i in range(n):
        for j in range(i, n):
            matrix[i][j] = matrix[j][i]
 
    # Проверяем, что матрица положительно определена
    is_positive_definite = True
    for k in range(1, n+1):
        submatrix = [row[:k] for row in matrix[:k]]
        determinant = 1
        for i in range(k):
            determinant *= submatrix[i][i]
            for j in range(i+1, k):
                ratio = submatrix[j][i] / submatrix[i][i]
                for l in range(i, k):
                    submatrix[j][l] -= ratio * submatrix[i][l]
        if determinant <= 0:
            is_positive_definite = False
            break
 
    # Замеряем время окончания операции
    end_time = time.time()
 
    if is_positive_definite:
        print("Сгенерирована положительно определенная матрица:")
        for row in matrix:
            print(row)
    else:
        print("Сгенерирована не положительно определенная матрица.")
 
    # Вычисляем и добавляем временные затраты в список
    times.append(end_time - start_time)
 
# Построение графика временных затрат
plt.plot(range(1, 11), times, 'o-')
plt.xlabel('Номер операции')
plt.ylabel('Время выполнения операции, сек')
plt.title('График временных затрат')
plt.show()