Skip to content

Instantly share code, notes, and snippets.

@Keno
Created April 3, 2021 02:16
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save Keno/2b54d5baa3f4cd6cc30310afec06c476 to your computer and use it in GitHub Desktop.
Save Keno/2b54d5baa3f4cd6cc30310afec06c476 to your computer and use it in GitHub Desktop.
julia> foreach(println, equations(sys))
0 ~ var"0"(t) - v1₊n₊v(t)
0 ~ var"0"(t) - r1₊n₊v(t)
0 ~ var"0"(t) - c1₊n₊v(t)
0 ~ var"0"(t)
0 ~ d1₊p₊i(t) + v1₊p₊i(t)
0 ~ net1(t) - v1₊p₊v(t)
0 ~ net1(t) - d1₊p₊v(t)
0 ~ c1₊p₊i(t) + d1₊n₊i(t) + r1₊p₊i(t)
0 ~ net2(t) - d1₊n₊v(t)
0 ~ net2(t) - r1₊p₊v(t)
0 ~ net2(t) - c1₊p₊v(t)
0 ~ c1₊n₊v(t) + c1₊vc(t) - c1₊p₊v(t)
0 ~ c1₊c*c1₊vc(t) - c1₊q(t)
Differential(t)(c1₊q(t)) ~ c1₊n₊i(t)
0 ~ c1₊n₊i(t) + c1₊p₊i(t)
0 ~ r1₊n₊v(t) + r1₊vr(t) - r1₊p₊v(t)
0 ~ r1₊vr(t)*(r1₊r^-1) - r1₊n₊i(t)
0 ~ r1₊n₊i(t) + r1₊p₊i(t)
0 ~ d1₊c₊p₊i(t) + d1₊d1₊p₊i(t) + d1₊g₊p₊i(t) - d1₊p₊i(t)
0 ~ d1₊n1(t) - d1₊p₊v(t)
0 ~ d1₊n1(t) - d1₊g₊p₊v(t)
0 ~ d1₊n1(t) - d1₊c₊p₊v(t)
0 ~ d1₊n1(t) - d1₊d1₊p₊v(t)
0 ~ d1₊c₊n₊i(t) + d1₊d1₊n₊i(t) + d1₊g₊n₊i(t) + d1₊r₊p₊i(t)
0 ~ d1₊n2(t) - d1₊d1₊n₊v(t)
0 ~ d1₊n2(t) - d1₊g₊n₊v(t)
0 ~ d1₊n2(t) - d1₊c₊n₊v(t)
0 ~ d1₊n2(t) - d1₊r₊p₊v(t)
0 ~ d1₊r₊n₊i(t) - d1₊n₊i(t)
0 ~ d1₊n3(t) - d1₊n₊v(t)
0 ~ d1₊n3(t) - d1₊r₊n₊v(t)
0 ~ d1₊d1₊n₊v(t) + d1₊d1₊vr(t) - d1₊d1₊p₊v(t)
0 ~ 4.352e-9(exp(20.28428190257992d1₊d1₊vr(t)) - 1) - d1₊d1₊n₊i(t)
0 ~ d1₊d1₊n₊i(t) + d1₊d1₊p₊i(t)
0 ~ d1₊c₊n₊v(t) + d1₊c₊vc(t) - d1₊c₊p₊v(t)
0 ~ 1.0e-6d1₊c₊vc(t) - d1₊c₊q(t)
Differential(t)(d1₊c₊q(t)) ~ d1₊c₊n₊i(t)
0 ~ d1₊c₊n₊i(t) + d1₊c₊p₊i(t)
0 ~ d1₊g₊n₊v(t) + d1₊g₊vr(t) - d1₊g₊p₊v(t)
0 ~ d1₊g₊g*d1₊g₊vr(t) - d1₊g₊n₊i(t)
0 ~ d1₊g₊n₊i(t) + d1₊g₊p₊i(t)
0 ~ d1₊r₊n₊v(t) + d1₊r₊vr(t) - d1₊r₊p₊v(t)
0 ~ d1₊r₊vr(t)*(d1₊r₊r^-1) - d1₊r₊n₊i(t)
0 ~ d1₊r₊n₊i(t) + d1₊r₊p₊i(t)
0 ~ v1₊p₊v(t) - v1₊vo - v1₊n₊v(t) - (v1₊va*sinpi(0.005555555555555556v1₊ϕ + 2t*v1₊freq))
0 ~ v1₊n₊i(t) + v1₊p₊i(t)
julia> foreach(println, states(sys))
var"0"(t)
net1(t)
net2(t)
c1₊q(t)
c1₊vc(t)
c1₊p₊v(t)
c1₊p₊i(t)
c1₊n₊v(t)
c1₊n₊i(t)
r1₊vr(t)
r1₊p₊v(t)
r1₊p₊i(t)
r1₊n₊v(t)
r1₊n₊i(t)
d1₊n1(t)
d1₊n2(t)
d1₊n3(t)
d1₊p₊v(t)
d1₊p₊i(t)
d1₊n₊v(t)
d1₊n₊i(t)
d1₊d1₊vr(t)
d1₊d1₊p₊v(t)
d1₊d1₊p₊i(t)
d1₊d1₊n₊v(t)
d1₊d1₊n₊i(t)
d1₊c₊q(t)
d1₊c₊vc(t)
d1₊c₊p₊v(t)
d1₊c₊p₊i(t)
d1₊c₊n₊v(t)
d1₊c₊n₊i(t)
d1₊g₊vr(t)
d1₊g₊p₊v(t)
d1₊g₊p₊i(t)
d1₊g₊n₊v(t)
d1₊g₊n₊i(t)
d1₊r₊vr(t)
d1₊r₊p₊v(t)
d1₊r₊p₊i(t)
d1₊r₊n₊v(t)
d1₊r₊n₊i(t)
v1₊p₊v(t)
v1₊p₊i(t)
v1₊n₊v(t)
v1₊n₊i(t)
julia> foreach(println, parameters(sys))
c1₊c
r1₊r
d1₊TT
d1₊RS
d1₊d1₊AREA
d1₊d1₊M
d1₊d1₊TEMP
d1₊d1₊N
d1₊d1₊BV
d1₊d1₊IBV
d1₊d1₊IK
d1₊d1₊IKR
d1₊d1₊IS
d1₊c₊AREA
d1₊c₊M
d1₊c₊PJ
d1₊c₊CJO
d1₊c₊CJP
d1₊c₊MJ
d1₊c₊MJSW
d1₊c₊TT
d1₊c₊PB
d1₊c₊PHP
d1₊c₊N
d1₊c₊BV
d1₊c₊IBV
d1₊c₊IK
d1₊c₊IKR
d1₊c₊IS
d1₊g₊g
d1₊r₊r
v1₊vo
v1₊va
v1₊freq
v1₊ϕ
julia> foreach(println, MTK.defaults(sys))
Pair{Any, Any}(d1₊g₊g, 1.0e-13)
Pair{Any, Any}(d1₊c₊IBV, 0.0001)
Pair{Any, Any}(d1₊d1₊TEMP, 27.0)
Pair{Any, Any}(d1₊d1₊IKR, 0.0)
Pair{Any, Any}(d1₊c₊IKR, 0.0)
Pair{Any, Any}(d1₊c₊MJ, 0.03)
Pair{Any, Any}(d1₊RS, 0.6458)
Pair{Any, Any}(d1₊d1₊BV, 110.0)
Pair{Any, Any}(d1₊d1₊N, 1.9060000000000001)
Pair{Any, Any}(d1₊c₊N, 1.9060000000000001)
Pair{Any, Any}(r1₊r, 10000.0)
Pair{Any, Any}(d1₊d1₊M, 1.0)
Pair{Any, Any}(d1₊c₊PB, 0.869)
Pair{Any, Any}(v1₊freq, 1000.0)
Pair{Any, Any}(d1₊d1₊IS, 4.352e-9)
Pair{Any, Any}(d1₊c₊PJ, 0.0)
Pair{Any, Any}(d1₊d1₊AREA, 1.0)
Pair{Any, Any}(d1₊d1₊IK, 0.0)
Pair{Any, Any}(d1₊c₊CJO, 7.048e-13)
Pair{Any, Any}(d1₊d1₊IBV, 0.0001)
Pair{Any, Any}(d1₊c₊M, 1.0)
Pair{Any, Any}(d1₊c₊AREA, 1.0)
Pair{Any, Any}(c1₊c, 9.999999999999999e-6)
Pair{Any, Any}(d1₊c₊PHP, 1.0)
Pair{Any, Any}(v1₊ϕ, 0.0)
Pair{Any, Any}(d1₊TT, 3.48e-9)
Pair{Any, Any}(v1₊va, 2.5)
Pair{Any, Any}(d1₊c₊CJP, 0.0)
Pair{Any, Any}(d1₊c₊MJSW, 0.33)
Pair{Any, Any}(v1₊vo, 2.5)
Pair{Any, Any}(d1₊c₊IS, 4.352e-9)
Pair{Any, Any}(d1₊c₊IK, 0.0)
Pair{Any, Any}(d1₊r₊r, 0.6458)
Pair{Any, Any}(d1₊c₊BV, 110.0)
Pair{Any, Any}(d1₊c₊TT, 3.48e-9)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment