Problema 4 Cinética

Introducción

La cinética química es el estudio de las velocidades de las reacciones químicas y de los mecanismos mediante los que tienen lugar. La cinética química introduce la variable tiempo en el estudio de las reacciones químicas y estudia el camino que siguen los reactivos para convertirse en productos.

Desarrollo

A=0.2
B=0.15
C=0
S=0.02
X=0
D=0
K1=0.1
K2=0.2
K3=0.05
dt=0.01
t=0

AA=c()
BB=c()
CC=c()
SS=c()
XX=c()
DD=c()
tt=c()

while (t<50) {
  AA=c(AA,A)
  BB=c(BB,B)
  CC=c(CC,C)
  SS=c(SS,S)
  XX=c(XX,X)
  DD=c(DD,D)
  tt=c(tt,t)
  
  dA=-K1*A*dt
  dB=-K1*B*dt
  dC=(K1*A*dt)*(K1*B*dt)
  dC=(K1*A*dt)+(K1*B*dt)-(K2*C*dt)
  dX=(K2*C*dt)*(K2*S*dt)
  dS=-K2*S*dt
  dD=-K3*D*dt
  dD=(K2*C*dt)+(K2*S*dt)-(K2*S*dt)
  dX=-K3*X*dt
  
  A=A+dA
  B=B+dB
  C=C+dC
  S=S+dS
  X=X+dX
  D=D+dD
  t=t+dt
  
}

plot(tt,AA,type="l",col="black",ylim=c(0,0.20),xlim=c(0,50),lwd=2,main="Cinética de Reacciones")
lines(tt,BB,type="l",lwd=2,col="deeppink")
lines(tt,CC,type="l",lwd=2,col="green4")
lines(tt,SS,type="l",lwd=2,col="blue4")
lines(tt,XX,type="l",lwd=2,col="purple3")
lines(tt,DD,type="l",lwd=2,col="red4")
legend(x = "topright",
       legend = c("[A]", "[B]","[C]","[S]","[X]","[D]"),
       lty = 1, col = c("black", "deeppink","green4","blue4","purple3","red4"), lwd = 2)