import matplotlib.pyplot as plt
import numpy as np
import numpy.random as rd
from scipy.integrate import odeint
from math import *

M_air = 29e-3
R=8.314
g=9.81
T0=
T1=
Tm=
P_init=   #à completer avec la pression au niveau de la mer
h=           #à completer avec la hauteur de la troposphère
z=np.linspace(0,h,10000)
def F_pression(P,z):
    return (  )

T=T0
P_0=odeint(F_pression,P_init,z)
T=T1
P_1=odeint(F_pression,P_init,z)
T=Tm
P_m=odeint(F_pression,P_init,z)

plt.plot(z,P_0,label="T_0 ")

plt.plot(z,P_1,label="T_1")
plt.plot(z,P_m,label="T_m")
plt.legend()
plt.show()

def F_pression2(P,z):
    return ( )

P_2=odeint(F_pression2,P_init,z)

plt.plot(z,P_2,label="température variable")
plt.legend()