Update compats

Signed-off-by: ErikQQY <2283984853@qq.com>

Project.toml

@@ -13,12 +13,12 @@ Roots = "f2b01f46-fcfa-551c-844a-d8ac1e96c665"
 SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
 
 [compat]
-ControlSystemIdentification = "2.2.4"
-ControlSystems = "0.11, 0.12"
+ControlSystemIdentification = "2.9"
+ControlSystems = "1.10"
 ForwardDiff = "0.10"
-FractionalDiffEq = "0.2"
-Roots = "1.3.14, 2"
-SpecialFunctions = "2"
+FractionalDiffEq = "0.3.1"
+Roots = "2"
+SpecialFunctions = "2.3"
 julia = "1.6"
 
 [extras]

examples/RFLE.jl

@@ -1,4 +1,4 @@
-using FractionalSystems, Plots
+using FractionalDiffEq, Plots
 
 function RF(du, u, p, t)
  du[1] = u[2]*(u[3]-1+u[1]*u[1])+0.1*u[1]
@@ -10,4 +10,9 @@ LE, tspan=FOLyapunov(RF, 0.999, 0, 0.02, 300, [0.1; 0.1; 0.1], 0.005, 1000)
 
 plot(tspan, LE[1, :])
 plot!(tspan, LE[2, :])
-plot!(tspan, LE[3, :])
+plot!(tspan, LE[3, :])
+
+# Or we can compute the Lyapunov exponents with FractionalDiffEq.jl
+using FractionalDiffEq
+sys=FODESystem(RF, [0.999, 0.999, 0.999], [0.1; 0.1; 0.1], (0, 300))
+LE, tspan = FOLyapunov(sys, 0.02, 0.005, 1000)

src/FOLE.jl

@@ -104,7 +104,7 @@ function FOLyapunov(fun, order, t_start, h_norm, t_end, u0, h, out)# TODO: Gener
  return LE, tspan
 end
 
-FOLyapunov(sys::FODESystem, h_norm, h, out) = FOLyapunov(sys.f, sys.α, tspan[1], h_norm, tspan[2], sys.x0, h, out)
+#FOLyapunov(sys::FODESystem, h_norm, h, out) = FOLyapunov(sys.f, sys.α, sys.tspan[1], h_norm, sys.tspan[2], sys.u0, h, out)
 
 mutable struct M
  an

src/fotf/show.jl

@@ -28,7 +28,7 @@ function polyshow(p, np)
  P = ""
  p = p[ii]
 
- for i = 1:length(p)
+ for i in eachindex(p)
  P = P*"+"*"$(p[i])"*"s^"*"{$(np[i])}"
  end