.46e-03
Iter   650, f(x) = 5.861851e+02, ||grad||_infty = 9.60e-03
Iter   660, f(x) = 5.861851e+02, ||grad||_infty = 2.05e-02
Iter   664, f(x) = 5.861851e+02, ||grad||_infty = 3.68e-02

exitmsg =

    'CONVERGENCE: REL_REDUCTION_OF_F_&lt;=_FACTR*EPSMCH.'


fit =

   99.0247

</pre><pre class="codeinput"><span class="comment">% Check the answer (1 is perfect)</span>
scr = score(M_plus, M_true)
</pre><pre class="codeoutput">
scr =

    0.9983

</pre><h2 id="12">Nonnegative factorization</h2><p>We can employ l