  We propose a novel method for a solution of a system of linear equations with
the non-negativity condition. The method is based on the Tikhonov functional
and has better accuracy and stability than other well-known algorithms.
