test_b = np.array([8, 18])
print(nnls(test_A, test_b))
How can I obtain the same 'special' solutions to underdetermined linear systems that Matlab's `A \ b` (mldivide) operator returns using numpy/scipy?