Для решения данной задачи в Python можно использовать библиотеку NumPy. Вот пример кода, который решает задачу:
```python
import numpy as np
# Задаем размерность матрицы
n = 5
# Создаем квадратную матрицу размером n на n, заполненную случайными числами от 1 до 200
matrix = np.random.randint(1, 201, (n, n))
# Выводим матрицу на экран
for row in matrix:
for element in row:
print(f"{element:4}", end="")
print()
# Вычисляем разность сумм элементов главной и второстепенной диагоналей
main_diagonal_sum = np.trace(matrix)
secondary_diagonal_sum = np.trace(np.fliplr(matrix))
difference = main_diagonal_sum - secondary_diagonal_sum
# Выводим разность на экран
print("Разность сумм элементов главной и второстепенной диагоналей:", difference)
```
Пример вывода:
```
87 181 174 124 145
58 64 141 69 109
17 98 57 166 85
77 18 151 142 57
186 92 10 37 94
Разность сумм элементов главной и второстепенной диагоналей: -161