def f1(y, t ):
    x, y = y
    return [x-x*y-5*x**2, y-x*y-5*y**2]
def dependXOfT1(i):
    y0 = [0.001 +i* 0.0001,0.001 +i * 0.0003]
    w = odeint(f1, y0, t)
    plt.plot(t, w[:, 0])
def dependYOfT1(i):
    y0 = [-0.001 +i* 0.0001,0.001 +i * 0.0003]
    w = odeint(f1, y0, t)
    plt.plot(t, w[:, 1])
t = np.linspace(0, 5, 100)
for i in range(0, 1):
    dependXOfT1(i)
plt.show()
t = np.linspace(0, 5, 100)
for i in range(0, 1):
    dependYOfT1(i)
plt.show()
def dependXOfT2(i):
    y0 = [0.001 + i * 0.0005, 0.5 + i * 0.001]
    w = odeint(f1, y0, t)
    plt.plot(t, w[:, 0])
def dependYOfT2(i):
    y0 = [0 + i * 0.0005, 0.6 - i * 0.001]
    w = odeint(f1, y0, t)
    plt.plot(t, w[:, 1])
t = np.linspace(0, 5, 100)
for i in range(0, 1):
    dependXOfT2(i)
plt.show()
t = np.linspace(0, 5, 100)
for i in range(0, 1):
    dependYOfT2(i)
plt.show()
def dependXOfT3(i):
    y0 = [0.6 + i * 0.001, 0 + i * 0.001]
    w = odeint(f1, y0, t)
    plt.plot(t, w[:, 0])
def dependYOfT3(i):
    y0 = [0.5 + i * 0.001, 0.001 + i * 0.001]
    w = odeint(f1, y0, t)
    plt.plot(t, w[:, 1])
t = np.linspace(0, 5, 100)
for i in range(0, 1):
    dependXOfT3(i)
plt.show()
t = np.linspace(0, 5, 100)
for i in range(0, 1):
    dependYOfT3(i)
plt.show()
def dependXOfT4(i):
    y0 = [0.1668, 0.1667+ i * 0.001]
    w = odeint(f1, y0, t)
    plt.plot(t, w[:, 0])
    #plt.plot(t, w[:, 1])
def dependYOfT4(i):
    y0 = [0.1667, 0.16]
    w = odeint(f1, y0, t)
    plt.plot(t, w[:, 1])
t = np.linspace(0, 5, 100)
for i in range(0, 1):
    dependXOfT4(i)
plt.show()
t = np.linspace(0, 5, 100)
for i in range(0, 1):
    dependYOfT4(i)
plt.show()