### Taking a large number of points and plotting them in a scatter plot

I'm solving newton's cubic for all the black points that occur within the fractal:

I have the code:

LA(x) = x - (x^3.0-1)/(3.0*x^2.0)
S = solve(LA==0,x)
for s in S:
solve(LA == s.rhs(),x)

and i want to take all those points I'm generating and plot them in a scatter plot. I think I have to use an array? Can anyone tell me exactly how I tell sage to put these points in an array then plot them? Thanks everyone for your help