Добрый день. На языке Python нужно разработать программу, которая посчитает SVD разложение (Сингулярное разложение) любой прямоугольной вещественной матрицы, размерность может быть 20х15 и более (до 100х100, но не надо это прямо ограничивать в коде), числа в матрице могут быть большими (6-значными и более). Без использования библиотек, только чистый Python, ни одного импорта!

Матрицу можно задать прямо в коде, но, естественно, алгоритм должен работать для любой прямоугольной и квадратной матрицы (в том числе с нулевым определителем).
Какой алгоритм разложения выбрать? Желательно самый простой, который подходит для условия задачи, это обсуждаемо.
В отклике просьба написать с помощью какого алгоритма вы бы выполнили задачу.