pts := {seq( (i/300)^3 / 2, i=0..300)}:
pts := {op(pts),seq(1-i, i=pts)}:
pts := {seq(fsolve(numer(f-i), x, complex), i = pts)}:
pts := {seq([Re(i), Im(i)], i = pts)};