-
Notifications
You must be signed in to change notification settings - Fork 0
/
master_file.R
36 lines (28 loc) · 995 Bytes
/
master_file.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
library(deSolve)
state.variables <- c(gp = 1E-9, p0 = 1E-9, xs = 0, m0 = 0,
r0 = 1E-9, mr = 0, z0 = 0,
AA = 1E-9, NT = 1E-9)
parm.values <- c(k1 = 6E9, k2 = 600, k3 = 1E10,
k4 = 6E9, k5 = 135, k6 = 1E10,
k7 = 18)
con.profile <- function(t, state.variables, parm.values){
with(as.list(c(state.variables, parm.values)),{
dgp <- -k1*gp*p0 + k2*xs + k3*NT*xs
dp0 <- -k1*gp*p0 + k2*xs + k3*NT*xs
dxs <- k1*gp*p0 - k2*xs - k3*NT*xs
dm0 <- k3*NT*xs - k4*m0*r0 + k5*mr + k6*AA*mr - k7*m0
dr0 <- -k4*m0*r0 + k5*mr + k6*AA*mr
dmr <- k4*m0*r0 - k5*mr - k6*AA*mr
dz0 <- k6*AA*mr
dAA <- -k6*AA*mr
dNT <- -k3*NT*xs
return(list(c(dgp, dp0, dxs, dm0, dr0, dmr, dz0,
dAA, dNT)))
})
}
out <- ode(y = state.variables,
parms = parm.values,
method = "bdf",
times = seq(0, 30, 0.001),
func = con.profile)
plot(out)