From 3719719da76bfc112afe860edd355ed77e765fb2 Mon Sep 17 00:00:00 2001
From: Nathanael Bosch <nathanael.bosch@uni-tuebingen.de>
Date: Sat, 23 Dec 2023 12:21:01 +0100
Subject: [PATCH] Better tests, and some bugfixes

---
 src/chaotic_attractors.jl       | 12 ++++++------
 test/test_chaotic_attractors.jl |  4 ++--
 2 files changed, 8 insertions(+), 8 deletions(-)

diff --git a/src/chaotic_attractors.jl b/src/chaotic_attractors.jl
index 84c069a..06fde79 100644
--- a/src/chaotic_attractors.jl
+++ b/src/chaotic_attractors.jl
@@ -3137,10 +3137,10 @@ function CellCycle()
         V3 = m1 * Vm3
         U3 = m2 * Um3
         du[1] = vi1 * Kim1 / (Kim1 + m2) - vd1 * x1 * c1 / (Kd1 + c1) - kd1 * c1
-        du[2] = vi2 * Kim2 / (Kim2 + m1) - vd2 * x2 * c2 / (Kd2 + c2) - kd2 * c2
-        du[3] = V1 * (1 - m1) / (K1 + (1 - m1)) - V2 * m1 / (K2 + m1)
-        du[4] = U1 * (1 - m2) / (H1 + (1 - m2)) - U2 * m2 / (H2 + m2)
-        du[5] = V3 * (1 - x1) / (K3 + (1 - x1)) - V4 * x1 / (K4 + x1)
+        du[4] = vi2 * Kim2 / (Kim2 + m1) - vd2 * x2 * c2 / (Kd2 + c2) - kd2 * c2
+        du[2] = V1 * (1 - m1) / (K1 + (1 - m1)) - V2 * m1 / (K2 + m1)
+        du[5] = U1 * (1 - m2) / (H1 + (1 - m2)) - U2 * m2 / (H2 + m2)
+        du[3] = V3 * (1 - x1) / (K3 + (1 - x1)) - V4 * x1 / (K4 + x1)
         du[6] = U3 * (1 - x2) / (H3 + (1 - x2)) - U4 * x2 / (H4 + x2)
     end
     u0 = Float64.(ATTRACTOR_DATA["CellCycle"]["initial_conditions"])
@@ -3428,8 +3428,8 @@ function ForcedVanDerPol()
     function rhs(du, u, p, t)
         x, y, z = u
         @unpack a, mu, w = p
-        du[1] = mu * (1 - x^2) * y - x + a * sin(z)
-        du[2] = y
+        du[2] = mu * (1 - x^2) * y - x + a * sin(z)
+        du[1] = y
         du[3] = w
     end
     u0 = Float64.(ATTRACTOR_DATA["ForcedVanDerPol"]["initial_conditions"])
diff --git a/test/test_chaotic_attractors.jl b/test/test_chaotic_attractors.jl
index 0a0df17..41677e0 100644
--- a/test/test_chaotic_attractors.jl
+++ b/test/test_chaotic_attractors.jl
@@ -130,8 +130,8 @@ using SafeTestsets
 
         prob = @test_nowarn System()
         @test prob isa ODEProblem
-        # @test_nowarn solve(prob, Tsit5())
-        @test solve(prob, Tsit5()) isa ODESolution
+        @test_nowarn solve(prob, Tsit5())
+        @test_nowarn solve(prob, Tsit5(), tspan=(0.0, 100.0), abstol=1e-6, reltol=1e-4)
     end
 
 end