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())
SteadyStateDiffEq.DynamicSS{OrdinaryDiffEqRosenbrock.Rodas5P{0, ADTypes.AutoForwardDiff{nothing, Nothing}, Nothing, typeof(OrdinaryDiffEqCore.DEFAULT_PRECS), Val{:forward}(), true, nothing, typeof(OrdinaryDiffEqCore.trivial_limiter!), typeof(OrdinaryDiffEqCore.trivial_limiter!)}, Float64}(OrdinaryDiffEqRosenbrock.Rodas5P{0, ADTypes.AutoForwardDiff{nothing, Nothing}, Nothing, typeof(OrdinaryDiffEqCore.DEFAULT_PRECS), Val{:forward}(), true, nothing, typeof(OrdinaryDiffEqCore.trivial_limiter!), typeof(OrdinaryDiffEqCore.trivial_limiter!)}(nothing, OrdinaryDiffEqCore.DEFAULT_PRECS, OrdinaryDiffEqCore.trivial_limiter!, OrdinaryDiffEqCore.trivial_limiter!, ADTypes.AutoForwardDiff()), Inf)
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.99955525872,
na_i(t) => 12237.06306889442,
vm(t) => -67.45750801042875,
(Cai(t))[43] => 0.1224283772185541,
(Cai(t))[42] => 0.12242837799871027,
(Cai(t))[41] => 0.12242837878548267,
(Cai(t))[40] => 0.12242837957915695,
(Cai(t))[39] => 0.1224283803800303,
(Cai(t))[38] => 0.12242838118841207,
(Cai(t))[37] => 0.12242838200462434,
(Cai(t))[36] => 0.12242838282900265,
(Cai(t))[35] => 0.12242838366189664,
(Cai(t))[34] => 0.12242838450367081,
(Cai(t))[33] => 0.12242838535470542,
(Cai(t))[32] => 0.12242838621539717,
(Cai(t))[31] => 0.12242838708616027,
(Cai(t))[30] => 0.1224283879674273,
(Cai(t))[29] => 0.12242838885965035,
(Cai(t))[28] => 0.12242838976330195,
(Cai(t))[27] => 0.12242839067887645,
(Cai(t))[26] => 0.1224283916068912,
(Cai(t))[25] => 0.12242839254788783,
(Cai(t))[24] => 0.12242839350243384,
(Cai(t))[23] => 0.12242839447112404,
(Cai(t))[22] => 0.12242839545458227,
(Cai(t))[21] => 0.12242839645346319,
(Cai(t))[20] => 0.12242839746845419,
(Cai(t))[19] => 0.12242839850027747,
(Cai(t))[18] => 0.1224283995496922,
(Cai(t))[17] => 0.1224284006174971,
(Cai(t))[16] => 0.12242840170453288,
(Cai(t))[15] => 0.12242840281168507,
(Cai(t))[14] => 0.12242840393988712,
(Cai(t))[13] => 0.12242840509012358,
(Cai(t))[12] => 0.12242840626343374,
(Cai(t))[11] => 0.12242840746091542,
(Cai(t))[10] => 0.12242840868372919,
(Cai(t))[9] => 0.12242840993310293,
(Cai(t))[8] => 0.12242841121033667,
(Cai(t))[7] => 0.12242841251680814,
(Cai(t))[6] => 0.12242841385397851,
(Cai(t))[5] => 0.12242841522339895,
(Cai(t))[4] => 0.12242841662671748,
(Cai(t))[3] => 0.12242841806568684,
(Cai(t))[2] => 0.12242841954217283,
(Cai(t))[44] => 0.1224283764447395,
(Cai(t))[1] => 0.12242842105816358,
CaMKOX(t) => 0.0,
CaMKAOX(t) => 0.0,
CaMKA2(t) => 0.0004420722457396739,
CaMKA(t) => 0.0017681844326322613,
CaMKPOX(t) => 0.0,
CaMKP(t) => 0.0008840674730976597,
CaMKBOX(t) => 0.0,
CaMKB(t) => 0.00798173236532062,
i_g(t) => 0.5604326537656363,
i_b(t) => 0.0039000845480040187,
i_fca(t) => 1.0034738182714618,
i_f(t) => 0.9997191242166692,
i_d(t) => 0.00039848916910502726,
i_Naj(t) => 0.4556048049152983,
i_Nah(t) => 0.3177101719979885,
i_Nam(t) => 0.025916958404658595,
i_y(t) => 0.14577203102220138,
i_IK(t) => 0.0005627933615767597,
i_OK(t) => 0.0017682976596783667,
i_CK2(t) => 0.001954706703028578,
i_CK1(t) => 0.0030255554426356097,
i_nKs(t) => 0.002312514978303931,
i_sslow(t) => 0.9995208752820927,
i_s(t) => 0.9995208621725907,
i_r(t) => 0.007693207255384156,
CaNSR(t) => 757.9220320937478,
CaJSR(t) => 757.7585107157387,
PO1RyR(t) => 0.0009533884511910048,

Runtime information

using InteractiveUtils
InteractiveUtils.versioninfo()
Julia Version 1.12.3
Commit 966d0af0fdf (2025-12-15 11:20 UTC)
Build Info:
  Official https://julialang.org release
Platform Info:
  OS: Linux (x86_64-linux-gnu)
  CPU: 12 × Intel(R) Core(TM) i7-8700 CPU @ 3.20GHz
  WORD_SIZE: 64
  LLVM: libLLVM-18.1.7 (ORCJIT, skylake)
  GC: Built with stock GC
Threads: 2 default, 1 interactive, 2 GC (on 12 virtual cores)
Environment:
  JULIA_CI = true
  JULIA_CONDAPKG_OFFLINE = true
  LD_LIBRARY_PATH = /home/github/actions-runner-1/_work/_tool/Python/3.14.2/x64/lib
  JULIA_NUM_THREADS = 2
  JULIA_PROJECT = /home/github/actions-runner-1/_work/camkii-cardiomyocyte-model/camkii-cardiomyocyte-model/Project.toml
  JULIA_DEPOT_PATH = /home/github/.julia:/home/github/actions-runner-1/_work/_tool/julia/1.12.3/x64/local/share/julia:/home/github/actions-runner-1/_work/_tool/julia/1.12.3/x64/share/julia
  JULIA_CONDAPKG_BACKEND = Null
  JULIA_LOAD_PATH = @:@v#.#:@stdlib
using Pkg
Pkg.status()
Project CaMKIIModel v0.7.0
Status `~/actions-runner-1/_work/camkii-cardiomyocyte-model/camkii-cardiomyocyte-model/Project.toml`
  [336ed68f] CSV v0.10.15
  [a93c6f00] DataFrames v1.8.1
  [459566f4] DiffEqCallbacks v4.11.0
  [f6369f11] ForwardDiff v1.3.1
  [682c06a0] JSON v1.3.0
  [23fbe1c1] Latexify v0.16.10
  [98b081ad] Literate v2.21.0
  [2fda8390] LsqFit v0.15.1
  [1862ce21] MarkdownTables v1.1.0
⌅ [961ee093] ModelingToolkit v10.31.2
  [77ba4419] NaNMath v1.1.3
  [1dea7af3] OrdinaryDiffEq v6.105.0
  [91a5bcdd] Plots v1.41.3
  [9672c7b4] SteadyStateDiffEq v2.8.0
  [bd369af6] Tables v1.12.1
  [ea8e919c] SHA v0.7.0
Info Packages marked with ⌅ have new versions available but compatibility constraints restrict them from upgrading. To see why use `status --outdated`

This notebook was generated using Literate.jl.

Back to top