@@ -76,6 +76,9 @@ def f(t):
return 1.0025+0.002*np.sin(1.25*3.14*t/30)
#return x ** 2
+def sig(x):
+ return 1/(1 + np.exp(-12*(x-0.431)))
+
strength_min = 0.4
strength_max = 0.6