Anharmonic treatment of vibrational resonance polyads—the diborane: a critical case for numerical methods