diff --git a/models/doping_3way_oscillator.jl b/models/doping_3way_oscillator.jl
new file mode 100644
index 0000000000000000000000000000000000000000..5e4be8eb6491b8eea92503409622f2efb7b7ab7e
--- /dev/null
+++ b/models/doping_3way_oscillator.jl
@@ -0,0 +1,16 @@
+
+doping_3way_oscillator = @network_model begin
+    R1: (A + B => 2B, rA*A*B)
+    R2: (B + C => 2C, rB*B*C)
+    R3: (C + A => 2A, rC*C*A)
+    R4: (DA + C => DA + A, rC*C*DA)
+    R5: (DB + A => DB + B, rA*A*DB)
+    R6: (DC + B => DC + C, rB*B*DC)
+end "doping 3way oscillator"
+
+set_x0!(doping_3way_oscillator, [:A,:B,:C,:DA,:DB,:DC], [333,333,333,10,10,10])
+set_param!(doping_3way_oscillator, [:rA,:rB,:rC], [1.0, 1.0, 1.0])
+set_time_bound!(doping_3way_oscillator, 0.2)
+
+export doping_3way_oscillator
+
diff --git a/models/repressilator.jl b/models/repressilator.jl
new file mode 100644
index 0000000000000000000000000000000000000000..dbd5a36aaee3f368ea1a1f85a085e7a15f95db4e
--- /dev/null
+++ b/models/repressilator.jl
@@ -0,0 +1,18 @@
+
+repressilator = @network_model begin
+    R1: (G1 => G1 + m1, α/(1+p3^n) + α0)
+    R2: (G2 => G2 + m2, α/(1+p1^n) + α0)
+    R3: (G3 => G3 + m3, α/(1+p2^n) + α0)
+    R4: (m1 => m1 + p1, β * m1)
+    R5: (m2 => m2 + p2, β * m2)
+    R6: (m3 => m3 + p3, β * m3)
+    R7: (m1 => 0, m1)
+    R8: (m2 => 0, m2)
+    R9: (m3 => 0, m3)
+end "Repressilator pkg"
+
+set_observed_var!(repressilator, [:m1, :m2, :m3, :p1, :p2, :p3])
+set_time_bound!(repressilator, 15.0)
+
+export repressilator
+