streamplot[{x(x^2-9)(x^2-16)^2,y(1-x^2)}]