Cinética Química

Cinética química

A=1.0
P=0
B=0
k1=0.1
k2=0.3
k3=0.2
dt=0.01
t= 0

AA=c()
tt=c()
PP=c()
BB=c()

while(t<50){
  
  AA=c(AA,A)
  tt=c(tt,t)
  PP=c(PP,P)
  BB=c(BB,B)
  
  dA=-k1*A*dt + k2*B*dt
  dB=k1*A*dt- k2*B*dt - k3*B*dt
  dP=k3*P*dt
  
  A=A+dA
  B=B+dB
  t=t+dt
  P=P+dP
  }
  plot(tt,AA,type="l", col="red",xlab="t(min)", ylab="[A]", main="Cinética química reacciones consecutivas con equilibrios", ylim=c(0,1.1),lwd=2.5)
lines(tt,PP,type="l",col="blue", lwd=2.5)
lines(tt,BB,type="l",col="green", lwd=2.5)