import numpy as np
import matplotlib.pyplot as plt
import scipy.integrate as scin

# Parametres de l'atmosphère et condition initiale en USI
M =   # masse molaire moyenne de l'air
g =  # champ de pesanteur
R =  # constante des GP
T0 =  # température moyenne en z = 0
alpha =  # gradient thermique
P0 =  # pression en z = 0

# Nombre d’iterations, altitude maximale et pas d'espace en USI
n =
zmax = 
dz = 

# Listes des altitudes
z = 
# Création d’une liste pour les pressions obtenues via la méthode d'Euler
P1 = [0 for i in range(n)]

# Initialisation
P1[0] = 
# Resolution par la methode d’Euler
for i in range(n-1):
    P1[i+1] = 

# Définition de la fonction F(P,z)
def F(P,z):
    return()

# Resolution par la fonction odeint
P2=

# Representation graphique
plt.plot(z,P1,z,P2)
plt.xlabel("z (m)")
plt.ylabel("P (bar)")
plt.legend(["Euler","odeint"],loc="best")
plt.grid()
plt.show()