using Model
using ModelingToolkit
using DifferentialEquations
using OrdinaryDiffEqSDIRK
using SteadyStateDiffEq
using Plots
Plots.default(lw=2)@time "Build system" sys = Model.DEFAULT_SYS
@time "Build problem" sprob = SteadyStateProblem(sys, [])
@time "Solve problem" sol = solve(sprob, DynamicSS(KenCarp4()); abstol=1e-8, reltol=1e-8)
for (k, v) in zip(unknowns(sys), sol.u)
println(k, " => ", v, ",")
endBuild system: 0.000023 seconds (12 allocations: 576 bytes)
Build problem: 6.387761 seconds (11.14 M allocations: 584.921 MiB, 1.77% gc time, 96.90% compilation time)
Solve problem: 7.763619 seconds (15.51 M allocations: 737.491 MiB, 2.83% gc time, 99.65% compilation time)
CaNSR(t) => 772.8474709441156,
CaJSR(t) => 772.6302279824113,
PO1RyR(t) => 0.0012456800353167691,
i_y(t) => 0.16332566868178358,
i_IK(t) => 0.000470985253856389,
i_OK(t) => 0.0015515735873586887,
i_CK2(t) => 0.0018361034782616277,
i_CK1(t) => 0.0028420487273897394,
i_nKs(t) => 0.002050251201390719,
i_sslow(t) => 0.999601646893252,
i_s(t) => 0.999600582204528,
i_r(t) => 0.007264580082392783,
i_Naj(t) => 0.5000636823477824,
i_Nah(t) => 0.36248896553993765,
i_Nam(t) => 0.022772305348641907,
i_g(t) => 0.5947420712105905,
i_b(t) => 0.003338233737987029,
i_fca(t) => 1.003321545026503,
i_f(t) => 0.9997597760801543,
i_d(t) => 0.0003544252676668834,
CaMKOX(t) => 0.0,
CaMKAOX(t) => 0.0,
CaMKA2(t) => 0.0005109732587419876,
CaMKA(t) => 0.0020091542224411198,
CaMKPOX(t) => 0.0,
CaMKP(t) => 5.41030013342864e-5,
CaMKBOX(t) => 0.0,
CaMKB(t) => 0.008411975499139548,
(Cai(t))[43] => 0.1253096100065126,
(Cai(t))[42] => 0.12530969621116692,
(Cai(t))[41] => 0.1253097829578949,
(Cai(t))[40] => 0.12530987027446136,
(Cai(t))[39] => 0.12530995818975307,
(Cai(t))[38] => 0.1253100467338366,
(Cai(t))[37] => 0.1253101359380188,
(Cai(t))[36] => 0.12531022583491133,
(Cai(t))[35] => 0.1253103164584997,
(Cai(t))[34] => 0.1253104078442159,
(Cai(t))[33] => 0.1253105000290165,
(Cai(t))[32] => 0.1253105930514653,
(Cai(t))[31] => 0.1253106869518218,
(Cai(t))[30] => 0.12531078177213517,
(Cai(t))[29] => 0.12531087755634499,
(Cai(t))[28] => 0.12531097435038854,
(Cai(t))[27] => 0.12531107220231583,
(Cai(t))[26] => 0.12531117116241197,
(Cai(t))[25] => 0.1253112712833292,
(Cai(t))[24] => 0.12531137262022726,
(Cai(t))[23] => 0.12531147523092484,
(Cai(t))[22] => 0.12531157917606156,
(Cai(t))[21] => 0.1253116845192719,
(Cai(t))[20] => 0.1253117913273728,
(Cai(t))[19] => 0.1253118996705647,
(Cai(t))[18] => 0.12531200962264855,
(Cai(t))[17] => 0.1253121212612596,
(Cai(t))[16] => 0.12531223466811978,
(Cai(t))[15] => 0.12531234992930987,
(Cai(t))[14] => 0.12531246713556424,
(Cai(t))[13] => 0.12531258638258952,
(Cai(t))[12] => 0.1253127077714098,
(Cai(t))[11] => 0.12531283140874122,
(Cai(t))[10] => 0.1253129574073985,
(Cai(t))[9] => 0.12531308588673676,
(Cai(t))[8] => 0.12531321697313258,
(Cai(t))[7] => 0.12531335080050793,
(Cai(t))[6] => 0.1253134875109016,
(Cai(t))[5] => 0.12531362725509387,
(Cai(t))[4] => 0.1253137701932893,
(Cai(t))[3] => 0.1253139164958646,
(Cai(t))[2] => 0.125314066344189,
(Cai(t))[44] => 0.1253095243172358,
(Cai(t))[1] => 0.1253142199315254,
k_i(t) => 152220.8892992791,
na_i(t) => 12574.549962143334,
vm(t) => -68.30153984045008,
This notebook was generated using Literate.jl.