Computing Closest Stable Non-negative Matrices