semestal problema 3

A=0.2
B=0.15
S=0.02
C=0.35
X=0.37
k1=0.1
k2=0.2
k3=0.05
dt=0.1
t=0

AA=c()
tt=c()
SS=c()
BB=c()

while(t<50){
  AA=c(AA,A)
  SS=c(SS,S)
  BB=c(BB,B)
  tt=c(tt,t)
  
  dA=-k1*A*dt+k2*B*dt
  dB=k1*A*dt-k2*B*dt-k3*C*dt*S*dt
  dS=k3*S*dt
  
    
  A=A+dA
  B=B+dB
  S=S+dS
  t=t+dt
}

plot(tt,AA,type="l",col="blue",ylim=c(0,1.1),xlab="t(min)",ylab="[A]",main="Cinética de segundo orden")
lines(tt,BB,type="l",col="green")
lines(tt,SS,type="l",col="red")