Plot[x^2+1/x,{x,-3,3},
Prolog->{PointSize[0.3],Red,Point[{2,3}]},Epilog->{PointSize[0.2],Green,Point[{-2,3}]}]