Вариант алгоритма Левенберга–Марквардта решения нелинейной задачи наименьших квадратов

Предлагается вариант реализации итерационного алгоритма Левенберга–Марквардта решения нелинейной задачи наименьших квадратов, использующий $QR$-разложение матрицы якобиана и не требующий факторизации матрицы при определении параметра Левенберга–Марквардта.