##Croissance de GaAs par épitaxie
##1. Evolution temporelle de la concentration en AsH3 , CA(t) 

import numpy as np
import matplotlib.pyplot as plt

##1ère partie : Evolution temporelle de [AsH3]= CA  
###Données
CA0= 1                 #concentration initiale en AsH3    mol/L
kd1 = 0.01             # constante de vitesse à T = 580K
Ea= 184                # energie d'activation en kJ/mol
R=8.314                #constante des gaz parfaits J/K/mol

###Calcul de kd à T =750 K
kd2 = 
print("valeur de kd à 750K =",kd2,"s-1")

### Détermination de CA(t) à 750 K  et tracé
# L'intervalle de temps à considérer est 0 < t < 0.1 s . 
from scipy.integrate import odeint

n=       #nombre de valeurs de la variable t dans l'intervalle retenu
T=np.linspace                   #a compléter 

réolution analytique
def CAan(t):
    

CA1=CAan(T)    
    
##résolution à l'aide de Odeint 
def 
    
    
solCA = odeint ( F, [CA0],T)

plt.figure(0)
plt.plot(T,solCA,'--k',label='CA(t)')
plt.plot(T,CA1,'+b',label='CA(T) analytique')
plt.legend()
plt.xlabel("t(s)")
plt.ylabel ("CA(mol/l)")
plt.show()


##2. Epaisseur et taux de croissance 
# Résoudre le système d'équations différentielles associées aux fonctions CT(t) , CM(t) et e(t)  en utilisant la fonction odeint . 
# Tracer le graphe montrant l'écolution de l'épaisseur au cours du temps . 
# L'intervalle des  temps à considérer ici est T2=[0 , 500s ]



#Données
CT0= 0.01                  #concentration initiale en Ga(CH3)3   mol/L
CM0=0                     #concentration initiale en  Ga(CH3)    mol/L
k1 = 5e-2                  # constante de vitesse à T = 580K
k2 =5e-4                   # energie d'activation en kJ/mol
alpha=10**4
mv= 5.3e6                 # masse volumique  en g/m3
M= 145                    #masse molaire de GaAs en g/mol

n=       #nombre de valeurs de la variable t dans l'intervalle retenu
T2=

# Système d'équations différentielles  : utiliser les notations introduites cu dessous 








sol= odeint (F2,[CT0,CM0,0],T2)
CT=sol[:,0]
CM=sol[:,1]
e=sol[:,2]

plt.figure(1)
plt.plot(T2,10**6*e,'--k',label='épaisseur en micromètres')
plt.legend()
plt.xlabel("t(s)")
plt.ylabel ("e")
plt.show()


## 3.Taux de croissance 

# Tracer le graphe montrant l'évolution du taux de croissance en fonction du temps 
def tau.........




plt.figure(2)
plt.plot(T2,Y,'xg',label='taux de croissance')
plt.legend()
plt.xlabel("t(s)")
plt.ylabel ("de/dt")
plt.show()


## 4. temps pour lequel le taux de croissance est maximal , tmax
 
# Faire afficher sur le script Python la valeur tmax obtenue à partir de la résolution analytique  . 
# Cette valeur peut être déterminée  en utilisant la fonction argmax . 
# Regarder l'influence du nombre n de valeurs choisi pour la variable t .
 





tmax = 
print ("tmax analytique = ",tmax,"s")


# utilisation du code précédent 
print ("tmax numérique=",T2[np.argmax(Y)],"s")