begin name
karin5
end name

begin reactions
v1 {1.0}$EXTERNAL = {1.0}g
v2 {1.0}$EXTERNAL = {1.0}ins
v3 {1.0}$EXTERNAL = {1.0}beta
v4 {1.0}$EXTERNAL = {1.0}mbeta
v5 {1.0}$EXTERNAL = {1.0}tamox
end reactions

begin rate equations
v1 = r0 - g[t]*(ego + si*ins[t])
v2 = (sigma*beta[t])/(1 + (alpha/g[t])^1.7) - gamma*ins[t] + (sigma*mbeta[t])/(1 + (alpha/(k*g[t]))^1.7)
v3 = (beta[t]*(-lambdamg + lambdapg - tamox[t]))/1440
v4 = ((-lambdamgmut + lambdapgmut)*mbeta[t] + beta[t]*tamox[t])/1440
v5 = -0.0007220283130832763*tamox[t]
end rate equations

begin parameters
alpha = 8.4
k = 6.0
default = 1.0
end parameters

begin initial conditions
beta[0] = betai
g[0] = gi
ins[0] = insi
mbeta[0] = mbetai
tamox[0] = tamoxi
end initial conditions

begin initial values
betai = 400.0
gi = 4.966667
insi = 11.42
mbetai = 0.0
tamoxi = 0.27
end initial values

begin assignment rules
lambdapg := 0.1/(1 + 37.262740006283494*(g[t]^(-1))^1.7)
lambdamg := 0.2*((1 + 11390625/g[t]^6)^(-1) + (1 + g[t]^8/65536)^(-1))
lambdamgmut := 0.2*((1 + 11390625/(k^6*g[t]^6))^(-1) + (1 + (k^8*g[t]^8)/65536)^(-1))
gamma := 3/10
r0 := 1/30
ego := 0.001
si := 0.0005
sigma := 0.030000000000000002
lambdapgmut := 0.1/(1 + 37.262740006283494*(1/(k*g[t]))^1.7)
perMutatedBcells := (100*mbeta[t])/beta[t]
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units