Statistics of Computation of damped nonlinear normal modes with internal resonances: a boundary value approach