Open ChrisRackauckas opened 1 year ago
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())