I have system and parametres:

```
function HR!(du, u, p, t)
function sigma(x)
return 1.0 / ( 1.0 + exp( -10.0 * ( x - ( - 0.25 ) ) ) )
end
a, b, c, d, s, xr, r, I, vs, k1, k2, el_link = p
x1, y1, z1, x2, y2, z2 = u
du[1] = y1 + b * x1 ^ 2 - a * x1 ^3 - z1 + I - k1 * ( x1 - vs ) * sigma(x2) + el_link * ( x2 - x1 )
du[2] = c - d * x1 ^2 - y1
du[3] = r * ( s * ( x1 - xr ) - z1 )
du[4] = y2 + b * x2 ^ 2 - a * x2 ^3 - z2 + I - k2 * ( x2 - vs ) * sigma(x1) + el_link * ( x1 - x2 )
du[5] = c - d * x2 ^2 - y2
du[6] = r * ( s * ( x2 - xr ) - z2 )
return SVector(du[1], du[2], du[3],
du[4], du[5], du[6])
end
tspan = (0., 50000.)
a = 1.
b = 3.
c = 1.
d = 5.
xr = -1.6
r = 0.01 # 0.01
s = 5.
I = 4.
xv = 2.
k1= -0.17
k2 = k1
```

I try calculate lyapunov spectrum with this code

```
p = [a, b, c, d,
s, xr, r, I, xv, k1, k2, 0]
initialcondition = [-0.469412067339509, -0.5839635169884322, 4.025111076102515, -0.469412067339509, -0.5839635169884322, 4.025111076102515]
ds_HR = ContinuousDynamicalSystem(HR!, initialcondition, p )
spectrum = lyapunovspectrum(ds_HR, tspan[2]; diffeq = (alg = AutoVern9(Rodas5()), maxiters = 10000000))
```

and get next error: