Nonlinear moment matching-based model order reduction