The PHENIX experiment at RHIC has a special ability to detect muons at the region from heavy quarks which are studied extensively to understand their production and modification in high energy hadron-hadron collisions. Quarkonia states like and of different binding energies are expected to have their yields modified differently by the medium or comoving particles. The is significantly more suppressed than in heavy-ion collisions by color screening effect. Quarkonia production in small collision systems like + collisions is also modified due to the initial state and final state effects, and a relative modification between and can provide important information on final state effects. We have been studied the multiplicity dependence of the production of and in + collisions at PHENIX to analyze the behavior of the potentially smallest QGP (collision system). In this poster, we will present the status of this study along with recent measurements of and in + collisions at GeV with other experiments and theoretical models.