Initial conditions

using ModelingToolkit
using OrdinaryDiffEq, SteadyStateDiffEq, DiffEqCallbacks
using Plots
using CaMKIIModel
Plots.default(lw=2)

sys = build_neonatal_ecc_sys(simplify=true, reduce_iso=true, reduce_camk=true)
prob = SteadyStateProblem(sys, [])
alg = DynamicSS(Rodas5P())

sol = solve(prob, alg; abstol=1e-10, reltol=1e-10)

for (k, v) in zip(unknowns(sys), sol.u)
    println(k, " => ", v, ",")
end

k_i(t) => 152559.99968427862,
na_i(t) => 12237.062940236663,
vm(t) => -67.4575076695362,
(Cai(t))[43] => 0.12242837623321787,
(Cai(t))[42] => 0.12242837701334203,
(Cai(t))[41] => 0.12242837780008216,
(Cai(t))[40] => 0.12242837859372387,
(Cai(t))[39] => 0.12242837939456436,
(Cai(t))[38] => 0.12242838020291295,
(Cai(t))[37] => 0.12242838101909173,
(Cai(t))[36] => 0.12242838184343618,
(Cai(t))[35] => 0.12242838267629598,
(Cai(t))[34] => 0.1224283835180356,
(Cai(t))[33] => 0.12242838436903525,
(Cai(t))[32] => 0.12242838522969163,
(Cai(t))[31] => 0.12242838610041896,
(Cai(t))[30] => 0.1224283869816498,
(Cai(t))[29] => 0.12242838787383617,
(Cai(t))[28] => 0.12242838877745067,
(Cai(t))[27] => 0.1224283896929876,
(Cai(t))[26] => 0.12242839062096424,
(Cai(t))[25] => 0.12242839156192224,
(Cai(t))[24] => 0.12242839251642906,
(Cai(t))[23] => 0.12242839348507949,
(Cai(t))[22] => 0.12242839446849736,
(Cai(t))[21] => 0.12242839546733726,
(Cai(t))[20] => 0.12242839648228655,
(Cai(t))[19] => 0.12242839751406741,
(Cai(t))[18] => 0.12242839856343904,
(Cai(t))[17] => 0.12242839963120009,
(Cai(t))[16] => 0.1224284007181912,
(Cai(t))[15] => 0.12242840182529789,
(Cai(t))[14] => 0.1224284029534536,
(Cai(t))[13] => 0.12242840410364281,
(Cai(t))[12] => 0.12242840527690474,
(Cai(t))[11] => 0.12242840647433724,
(Cai(t))[10] => 0.12242840769710077,
(Cai(t))[9] => 0.12242840894642315,
(Cai(t))[8] => 0.12242841022360443,
(Cai(t))[7] => 0.1224284115300222,
(Cai(t))[6] => 0.12242841286713757,
(Cai(t))[5] => 0.1224284142365017,
(Cai(t))[4] => 0.12242841563976256,
(Cai(t))[3] => 0.12242841707867275,
(Cai(t))[2] => 0.12242841855509803,
(Cai(t))[44] => 0.12242837545943497,
(Cai(t))[1] => 0.12242842007102644,
CaMKOX(t) => 0.0,
CaMKAOX(t) => 0.0,
CaMKA2(t) => 0.0006225205847209173,
CaMKA(t) => 0.0024900011407603785,
CaMKPOX(t) => 0.0,
CaMKP(t) => 0.001244981353589065,
CaMKBOX(t) => 0.0,
CaMKB(t) => 0.009189828208149755,
i_g(t) => 0.5604326397678525,
i_b(t) => 0.003900084792963966,
i_fca(t) => 1.0034738183232468,
i_f(t) => 0.9997191241989334,
i_d(t) => 0.00039848918796449393,
i_Naj(t) => 0.45560478714105657,
i_Nah(t) => 0.3177101545122867,
i_Nam(t) => 0.025916959753948296,
i_y(t) => 0.14577202425452118,
i_IK(t) => 0.0005627934020656607,
i_OK(t) => 0.0017682977530854307,
i_CK2(t) => 0.00195470675242204,
i_CK1(t) => 0.0030255555190542103,
i_nKs(t) => 0.002312515090567276,
i_sslow(t) => 0.9995208752463377,
i_s(t) => 0.999520862137374,
i_r(t) => 0.007693207433471413,
CaNSR(t) => 757.9220292965225,
CaJSR(t) => 757.7585079294614,
PO1RyR(t) => 0.0009533883894018224,

---

This notebook was generated using Literate.jl.