def ecosystem(state, t):
    """ The predator-prey system described with ODEs. """
    
    # Unpack current values of x and y
    x, y = state

    # Differential equations
    d_x = x*(2 - y - x)
    d_y = -y*(1 - 1.5*x)

    return [d_x, d_y]