Problema 4-Semestral

Problema 4

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)
  dS=-K2*S*dt
  dX=(K2*C*dt)+(K2*S*dt)
  dX=-K3*X*dt
  dS=K3*S*dt
  dD=K3*D*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="red",ylim=c(0,0.37))
lines(tt,BB,type="l",col="blue")
lines(tt,CC,type="l",col="green")
lines(tt,SS,type="l",col="orange")
lines(tt,XX,type="l",col="purple")
lines(tt,DD,type="l",col="black")