# valeurs à entrer
I2_0 = input('quantité initiale en diiode en mol :')
nI2_0 = float(I2_0)
S2O3_0 = input('quantité initiale en thiosulfate en mol :')
nS2O3_0 = float(S2O3_0)

limitant =''  
x=0  # avancement initial
a=0.001 #pas d'avancement

qI2=[nI2_0]
qS2O3=[nS2O3_0]

# écriture de la boucle while
while qI2[-1] > 0 and qS2O3[-1] > 0:
    x=x+a
    qI2.append(nI2_0-x)
    qS2O3.append(nS2O3_0-2*x)

#résolution du problème et affichage du résultat
if qI2[-1]<=0:
    limitant = 'le diiode'
if qS2O3[-1]<=0:
    limitant = 'le thiosulfate'
if qI2[-1]<=0 and qS2O3[-1]<=0 :
    limitant='à la fois le thiosulfate et le diiode : le mélange est stoechiometrique'

#print(limitant)
print('Le réactif limitant est ',limitant,'\n Avancement maximum : ',round(x,2),'mol' )