Update reactionsystem.jl

Author: ChrisRackauckas

test/reactionsystem.jl

@@ -1,4 +1,4 @@
-using ModelingToolkit, LinearAlgebra
+using ModelingToolkit, LinearAlgebra, Test
 
 @parameters t k[1:13]
 @variables A(t) B(t) C(t) D(t)
@@ -65,4 +65,4 @@ sf = SDEFunction{false}(sdesys, states(rs), parameters(rs))
 du2 = sf.f(u,p,t)
 @test norm(du-du2) < 100*eps()
 G2 = sf.g(u,p,t)
-@test norm(G-G2) < 100*eps()
+@test norm(G-G2) < 100*eps()