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*5*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="purple",ylim=c(0,0.37))
lines(tt,BB,type="l",col="yellow")
lines(tt,CC,type="l",col="orange")
lines(tt,SS,type="l",col="blue")
lines(tt,XX,type="l",col="green")
lines(tt,DD,type="l",col="red")