We give a complete classification of the phase portraits in Poincare´ disk for cubic polynomial systemsx˙=1−x−axy2,y˙=−by+axy2, R2 according with values its two parameters and b. These differential systems correspond to Gray-Scott model. Moreover we provide bifurcation diagram parameter plane (a,b) these systems.