On using triangular matrix factorization in the method of penalty functions for solution of problems in linear dynamic programming

A finite algorithm is proposed for solution of problems in linear dynamic programming that relies on the method of a quadratic penalty function with triangular factorization of unit matrices.