Drawing phase diagrams in Matlab with x and y isoclines and solution curves

Question

Good day, I plotted a phase diagram in Matlab for dx/dt = -x and dy/dt = (y-1)^2 * y * (x-1). Is there a better way of doing this and maybe adding some initial conditions and