@@ -1196,7 +1196,7 @@ end
11961196"""
11971197```julia
11981198hindmarshrose_two_coupled(u0=[0.1, 0.2, 0.3, 0.4, 0.5, 0.6];
1199- a = 1.0, b = 3.0, d = 5.0, r = 0.001, s = 4.0, xr = -1.6, I = 4.0,
1199+ a = 1.0, b = 3.0, c=1.0, d = 5.0, r = 0.001, s = 4.0, xr = -1.6, I = 4.0,
12001200k1 = -0.17, k2 = -0.17, k_el = 0.0, xv = 2.0)
12011201```
12021202```math
@@ -1214,7 +1214,7 @@ The default parameter values are taken from article "Dragon-king-like extreme ev
12141214coupled bursting neurons", DOI:https://doi.org/10.1103/PhysRevE.97.062311.
12151215"""
12161216function hindmarshrose_two_coupled (u0= [0.1 , 0.2 , 0.3 , 0.4 , 0.5 , 0.6 ];
1217- a = 1.0 , b = 3.0 , d = 5.0 , r = 0.001 , s = 4.0 , xr = - 1.6 , I = 4.0 ,
1217+ a = 1.0 , b = 3.0 , c = 1.0 , d = 5.0 , r = 0.001 , s = 4.0 , xr = - 1.6 , I = 4.0 ,
12181218 k1 = - 0.17 , k2 = - 0.17 , k_el = 0.0 , xv = 2.0 )
12191219 return CoupledODEs (hindmarshrose_coupled_rule, u0, [a, b, c, d, r, s, xr, I, k1, k2, k_el, xv])
12201220end
@@ -1225,11 +1225,11 @@ function hindmarshrose_coupled_rule(u, p, t)
12251225 a, b, c, d, r, s, xr, I, k1, k2, k_el, xv = p
12261226 x1, y1, z1, x2, y2, z2 = u
12271227
1228- du1 = y1 + b * x1 ^ 2 - a * x1 ^ 3 - z1 + I - k1 * ( x1 - vs ) * sigma (x2) + el_link * ( x2 - x1 )
1228+ du1 = y1 + b * x1 ^ 2 - a * x1 ^ 3 - z1 + I - k1 * ( x1 - xv ) * sigma (x2) + k_el * ( x2 - x1 )
12291229 du2 = c - d * x1 ^ 2 - y1
12301230 du3 = r * ( s * ( x1 - xr ) - z1 )
12311231
1232- du4 = y2 + b * x2 ^ 2 - a * x2 ^ 3 - z2 + I - k2 * ( x2 - vs ) * sigma (x1) + el_link * ( x1 - x2 )
1232+ du4 = y2 + b * x2 ^ 2 - a * x2 ^ 3 - z2 + I - k2 * ( x2 - xv ) * sigma (x1) + k_el * ( x1 - x2 )
12331233 du5 = c - d * x2 ^ 2 - y2
12341234 du6 = r * ( s * ( x2 - xr ) - z2 )
12351235 return SVector (du1, du2, du3, du4, du5, du6)
0 commit comments