#On détermine laquelle des deux hypothèses est la bonne
#On détermine alors le réactif limitant et on l'affiche
if xf1<xf2 :
    xf=xf1
    print("A est le réactif limitant")
if xf1>xf2 :
    xf=xf2
    print("B est le réactif limitant")
if xf1==xf2 :
    xf=xf1
    print("A et B ont été introduits dans les proportions stoechiométriques")

#On calcule les quantités de réactifs et de produits dans l'état final
nfA=n0A-a*xf
nfB=n0B-b*xf
nfC=n0C+c*xf
nfD=n0D+d*xf

#On affiche la composition de l'état final
print("Composition de l'état final:")
print("nf(A)="+str(nfA+" mol")
print("nf(B)="+str(nfB)+" mol")
print("nf(C)="+str(nfC)+" mol")
print("nf(D)="+str(nfD)+" mol")