def func(y, t ):
    x, y = y
    return [x-x*y-2*x**2, y-x*y-2*y**2]
def portrait1(i):  # Будуємо фазовий портер для точки x = 0, y = 0
    y0 = [0.00001 -i * 0.00001,-0.00013 -i * 0.00003]
    w = odeint(func, y0, t)
    plt.plot(w[:, 0], w[:, 1])

def portrait2(i):   # Будуємо фазовий портер для точки x = 0, y = 1/a , при a = 2
    y0 = [0 + i * 0.0005, 0.5 + i * 0.001]
    w = odeint(f1, y0, t)
    plt.plot(w[:, 0], w[:, 1])

def portrait3(i):   # Будуємо фазовий портер для точки x = 1/a, y = 0 , при a = 2
    y0 = [0.5 + i * 0.001, 0 + i * 0.001]
    w = odeint(f1, y0, t)
    plt.plot(w[:, 0], w[:, 1])

def portrait4(i):   # Будуємо фазовий портер для точки x = (a -1)/(a^2-1), y = -a(a -1)/(a^2-1)+1 , при a = 2
    y0 = [0.16667 + i * 0.0001, 0.168 + i * 0.001]
    w = odeint(f1, y0, t)
    plt.plot(w[:, 0], w[:, 1])

t = np.linspace(-1, 1, 100)
for i in range(-5, 5):
    portrait1(i)
    portrait2(i)
    portrait3(i)
    portrait4(i)
plt.show()