R_2-nullcline

R_1-nullcline

R_1-nullcline

R_2-nullcline

 stable 
direction

unstable
direction

I: 0
W: 1.5
M: 100
theta: 60
sigma: 10
S(x) = M / (1 + ℯ^(((-(x - theta))) / sigma))
R2nullcline = CurveCartesian(t, S((W * t) + I), t, 0, 200)
R1nullcline = CurveCartesian(S((W * t) + I), t, t, 0, 200)
A = Intersect(R2nullcline, R1nullcline, 4.034695969755146E-6, 4.034695969755146E-6)
B: (36.23098, 36.23098)
C: (99.98761, 99.98761)
D = (78.14918458974459, 17.83598671989934)
E = (67.22208487151445, 26.03131150857195)
u = Vector(D, E)
F = (62.213830833992326, 64.73145634397042)
G = (71.31974726585078, 76.56914770538641)
v = Vector(F, G)
H = (17.13954449629296, 88.40683906680239)
I_1 = (23.968981820186798, 79.30092263494393)
w = Vector(H, I_1)
J = (139.61412050478918, 125.7410964374221)
K = (125.95524585700151, 116.17988418397071)
a = Vector(J, K)
L = (7.123036421248666, 6.908887001669189)
M_1 = (18.050136139478813, 15.559507611934727)
N: (80, 10)
O = (67.22208487151448, 24.665424043793188)
P = (62.213830833992326, 38.779594513173805)
Q = (62.213830833992326, 49.25139840981103)
R = (65.8561974067357, 59.723202306448265)
T = (75.87270548178, 73.83737277582888)
c = PolyLine(N, O, P, Q, R, T, C)
b = Vector(M_1, L)
U: (50, 10)
V = (38.99374393275326, 16.925395076713496)
W_1 = (25.334869284965578, 17.835986719899342)
Z = (13.497177923549584, 12.372436860784266)
h = PolyLine(U, V, W_1, Z, A)
A_1 = Point(R1nullcline)
B_1 = (82.70214280567384, 64.2761605223775)
i = Vector(A_1, B_1)
C_1 = Point(R2nullcline)
D_1 = (66.76678904992156, 82.94328920768734)
j = Vector(C_1, D_1)
E_1 = (11.675994637177894, 26.48660733016488)
F_1 = (11.675994637177894, 16.470099255120573)
k = Vector(E_1, F_1)
G_1 = (25.790165106558497, 10.09595775281965)
H_1 = (15.773657031514198, 10.551253574412573)
l = Vector(G_1, H_1)
I_2 = Point(R2nullcline)
J_1 = (132.78468318089534, 100.24453042821841)
m = Vector(I_2, J_1)
K_1 = Point(R1nullcline)
L_1 = (100.00338402620491, 110.71633432485564)
n = Vector(K_1, L_1)
f = Vector(B + (8 * (1, -1)), B)
d = Vector(B, B + (8 * (1, 1)))
e = Vector(B, B - (8 * (1, 1)))
g = Vector(B - (8 * (1, -1)), B)