Finishing Scenario 1

[ChrisRackauckas]

• Add stratified petri net model
• Question 3
• Question 4
• Question 5
• Question 6

[ChrisRackauckas]

This is that SBML model read in using SBMLToolkit:

using SBMLToolkit, ModelingToolkit, OrdinaryDiffEq

xmlfile = raw"C:\Users\accou\.julia\dev\EasyModelAnalysis\docs\src\assets\Carcione2020.xml"
SBMLToolkit.checksupport_file(xmlfile)
mdl = readSBML(xmlfile, doc -> begin
    set_level_and_version(3, 2)(doc)
    convert_simplify_math(doc)
end)

rs = ReactionSystem(mdl)  # If you want to create a reaction system
odesys = convert(ODESystem, rs)  # Alternatively: ODESystem(mdl)
sys = structural_simplify(odesys)

tspan = (0., 1.)
prob = ODEProblem(odesys, [], tspan, [])
sol = solve(prob, Rodas5())

Carcione2020.zip