Nonnegative Matrix Factorization based on the Solution of Systems of ODEs