Пусть

- целочисленная симметрическая матрица порядка 150.
Требуется проверить ее на положительную определенность. Это можно сделать разными способами:
Код:
Timing@PositiveDefiniteMatrixQ@A
Timing@Reduce[LeadingPrincipalMinors@A > 0]
Timing[Min@N@Eigenvalues@A > 0]
Здесь LeadingPrincipalMinors - вывод всех угловых миноров

.
С помощью встроенной команды PositiveDefiniteMatrixQ ответ был дан за 0 секунд. Во втором случае примерно 3 секунды, в третьем случае примерно за 10 секунд.
Ясно, что всегда стоит использовать встроенные функции, но есть ли возможность посмотреть, какой алгоритм они используют? И вообще как они запрограммированы?