Wolfram Language Paclet Repository
Community-contributed installable additions to the Wolfram Language
RobertNachbar`CompartmentalModeling`
StochasticSolve  |  |
|
| | ||||
Details and Options
Examples
(4)
Basic Examples
(1)
Define the model and rate constants, and and compute the model's properties:
In[1]:=
model=ℰ,ℰℐ,ℐℛ;params={β0.000714,ζ0.196,γ0.143,1000,I010};modelData=
[model,t]
βλ[t]
→
ζ
→
γ
→
 | KineticCompartmentalModel |
Out[1]=
Variables{ℰ,ℐ,ℛ,},ComponentTransitionsℰℐ,ℐℛ,ℰ,StoichiometricMatrix{{-1,1,0,0},{0,-1,1,0},{1,0,0,-1}},Rates{ζℰ[t],γℐ[t],β[t]λ[t]},Equations{[t]-ζℰ[t]+β[t]λ[t],[t]ζℰ[t]-γℐ[t],[t]γℐ[t],[t]-β[t]λ[t]},TimeSymbolt
ζ
→
γ
→
βλ[t]
→
′
ℰ
′
ℐ
′
ℛ
′
Define the initial state vector and simulate the system:
In[2]:=
state=Replace[modelData["Variables"],{-I0,ℐI0,_0},{1}];
[modelData["Rates"]/.λ[t]ℐ[t]/.params,modelData["StoichiometricMatrix"],state/.params,modelData["Variables"],{t,0,150}]
 | StochasticSolve |
Out[2]=
ℰInterpolatingFunction
,ℐInterpolatingFunction
,ℛInterpolatingFunction
,InterpolatingFunction
 |  |
|
|  |
|
|  |
|
|  |
|
Options
(2)
Applications
(1)
 |
|
 |
|
""