I want to produce a phase portrait of the system:
x\'=x(y-1) and y\'=4-y^2-x^2
My code below is what I have attempted so far. My issue is the ics and duration
