Ejemplo 2: Calculo de las concentraciones en equilibrio (R-script)

#A 1280 °C la constante de equilibrio (Kc) para la reacción:
#Br2 (g)    <-->      2Br (g)
#es = 1.1 x 10-3.  
Si la concentración inicial es [Br2] = 0.063 M 
# y [Br] = 0.012 M,  calcula las concentraciones de las especies en equilibrio:
# Las concentraciones de equilibrio estan dadas por: 
#[Br2]eq= 0.063-x
#[2Br]eq= 0.012+2x
#Entonces Kc=(0.012+2x)^2)/(0.063-x)=1.1e-3
# Resuelve  x:
# Solución para la ecuación cuadrática:
#x= ((-b)+-sqrt((b^2)-(4*a*c)))/(2*a)
# Si:  4x2 + 0.0491x + 0.0000747 = 0 
#  Entonces 
a=4
b=0.0491
c=0.0000747
# Resolvindo la expresión: sqrt((b^2)-(4*a*c))=y  , entonces: 
y=(sqrt((b^2)-(4*a*c)))
#Tenemos que
x1=((-b)+(y))/(2*a)
x2=((-b)-(y))/(2*a)
print(x1)
print(x2)
# Para calcular las concentraciones en equilibrio [Br]=0.012+2x
Brx1=0.012+(2*x1)
Brx2=0.012+(2*x2)
print(Brx1)
print(Brx2)
# Debido a que Brx2 es negativa y este es un resultado imposible, se considera como Br1 la concentracion correcta y se considera entonces el valor de "x1" para calcular la concentracion de [Br2]

Br2=0.063+x1

print (Br2)

================================

Fin del script 

Encuentra los valores para las concentraciones de Br y Br2 en equilibrio: