pts = {Re[z], Im[z]} /. 
        NSolve[(z^9) + 2  (z^5) + 4  (z^3) + 9 == 0]

lst = 
      ListPlot[pts, AspectRatio -> 1, 
        PlotStyle -> PointSize[ .03]]

circ1 =  
        ParametricPlot[{Cos[t], Sin[t]}, {t, 0, 
            2  π}];
circ2 =  
        ParametricPlot[2 {Cos[t], Sin[t]}, {t, 0, 2  π}];

Show[lst, circ1, circ2]