from scipy.integrate import odeint
import scipy.integrate as integr
import matplotlib.pyplot as plt
import numpy as np

def f(y,t):
    dydt = y
    return dydt

T1 = np.arange(0, 2.01, 0.01) # à gauche, la valeur initiale; 2.01 est exclu
Y1 = odeint(f, 1, T1)
plt.plot(T1,Y1)

T2 = np.arange(0, -2.01, -0.01) # à gauche, la valeur initiale
Y2 = odeint(f, 1, T2)
plt.plot(T2,Y2)
plt.show()

##
T1 = np.linspace(0, 2, 100) # à gauche, la valeur initiale; 2 est inclus
Y1 = odeint(f, 1, T1)
plt.figure()
plt.plot(T1,Y1)

T2 = np.linspace(0, -2, 100) # à gauche, la valeur initiale
Y2 = odeint(f, 1, T2)
plt.plot(T2,Y2)
plt.show()

##
def f(X,t):
    dXdt = [X[1], (1/(1-t))*X[0]]
    return dXdt

T1 = np.arange(0, 0.96, 0.01) # à gauche, la valeur initiale
X1 = odeint(f, [0,1], T1)
plt.figure()
plt.plot(T1,X1[:,0])

T2 = np.arange(0, -2.01, -0.01) # à gauche, la valeur initiale
X2 = odeint(f, [0,1], T2)
plt.plot(T2,X2[:,0])
plt.show()