RecipePythonBeginner Draw a phase portrait of a two-variable ODE system leaving = sum(np.any((x < xlim[0]) | (x > xlim[1]) | (y < ylim[0]) | (y > ylim[1]))