%y

%y

%y

%y

%y

t

y

%y

t

y

ymax2 = 5.349622232407849
ymin2 = -1.4203355405805136
yscale2: ymax2 - ymin2
threshold: 0.01
threshold_y2: (threshold * yscale2)
dy: 1
xmax2 = 7.012841091492777
ymax1 = 3.571184894694351
ymin1 = -1.5629868121604424
yscale1: ymax1 - ymin1
threshold_y1: (threshold * yscale1)
c2point1 = DynamicCoordinates(c2point1_sub, 0, Min[ymax2, Max[ymin2, If[abs((round(y(c2point1_sub) / dy) * dy) - y(c2point1_sub)) < threshold_y2, (round(y(c2point1_sub) / dy) * dy), y(c2point1_sub)]]])
c1point3_sub = (1.6853932584269664, 1.9048313709038034)
c1point2_sub = (5.215088282504013, 0.733879477748281)
xmax1 = 7.012841091492777
c2point3 = DynamicCoordinates(c2point3_sub, Min[(xmax2 * 0.8), Max[0, x(c2point3_sub)]], Min[ymax2, Max[ymin2, If[(x(c2point3_sub) > Min[(xmax2 * 0.8), Max[0, x(c2point2_sub)]] ∨ abs(y(c2point3_sub) - y(c2point1)) < threshold_y2) ∧ abs((round(y(c2point3_sub) / dy) * dy) - y(c2point3_sub)) < threshold_y2, (round(y(c2point3_sub) / dy) * dy), y(c2point3_sub)]]])
c1point1_sub = (1617.1987947306184, 1.281632864384151)
c1point1 = DynamicCoordinates(c1point1_sub, 0, Min[ymax1, Max[ymin1, If[abs((round(y(c1point1_sub) / dy) * dy) - y(c1point1_sub)) < threshold_y1, (round(y(c1point1_sub) / dy) * dy), y(c1point1_sub)]]])
c1point3 = DynamicCoordinates(c1point3_sub, Min[(xmax1 * 0.8), Max[0, x(c1point3_sub)]], Min[ymax1, Max[ymin1, If[(x(c1point3_sub) > Min[(xmax1 * 0.8), Max[0, x(c1point2_sub)]] ∨ abs(y(c1point3_sub) - y(c1point1)) < threshold_y1) ∧ abs((round(y(c1point3_sub) / dy) * dy) - y(c1point3_sub)) < threshold_y1, (round(y(c1point3_sub) / dy) * dy), y(c1point3_sub)]]])
c1point2 = DynamicCoordinates(c1point2_sub, Min[(xmax1 * 0.8), Max[0, x(c1point2_sub)]], Min[ymax1, Max[ymin1, If[(x(c1point2_sub) > Min[(xmax1 * 0.8), Max[0, x(c1point3_sub)]] ∨ abs(y(c1point2_sub) - y(c1point1)) < threshold_y1) ∧ abs((round(y(c1point2_sub) / dy) * dy) - y(c1point2_sub)) < threshold_y1, (round(y(c1point2_sub) / dy) * dy), y(c1point2_sub)]]])
xmin1 = -1.0128410914927768
xscale: xmax1 - xmin1
threshold_x: (threshold * xscale)
xlabel1displayed: "" + (LaTeX[xlabel1]) + ""
ylabel1displayed: "" + (LaTeX[ylabel1]) + ""
calpha: -10 / xmax1
t12: x(c1point2)
t13: x(c1point3)
y11: y(c1point1)
y12: y(c1point2)
y13: y(c1point3)
b1alpha: ((((y12 - y11) * t13) * ℯ^((calpha * t13))) - (((y13 - y11) * t12) * ℯ^((calpha * t12)))) / ((((ℯ^((calpha * t12)) - 1) * t13) * ℯ^((calpha * t13))) - (((ℯ^((calpha * t13)) - 1) * t12) * ℯ^((calpha * t12))))
a1alpha: (y12 - y11 - ((ℯ^((calpha * t12)) - 1) * b1alpha)) / ((t12 * ℯ^((calpha * t12))))
f1alpha(x) = y11 - b1alpha + (((a1alpha * x) + b1alpha) * ℯ^((calpha * x)))
r1(x) = (y12 - y11) / (ℯ^((t12 * x)) - 1) - (y13 - y11) / (ℯ^((t13 * x)) - 1)
c1exp: x(Element[RemoveUndefined[{Root[r1, 1.0E-4], Root[r1, -1.0E-4]}], 1])
b1exp: (y12 - y11) / (ℯ^((c1exp * t12)) - 1)
f1exp(x) = y11 - b1exp + (b1exp * ℯ^((c1exp * x)))
c1monotonic = Defined(b1exp)
f1(x) = If[x ≥ 0, If[c1monotonic, f1exp, f1alpha]]
y1right = If(t12 > t13, y12, y13)
yinf1 = If(c1exp > 1 / xmax1, If[f1(xmax1) > ymax1, Infinity, If[f1(xmax1) < ymin1, -Infinity]], If[abs(y1right - f1(xmax1)) < threshold_y1, y1right])
c2point2 = DynamicCoordinates(c2point2_sub, Min[(xmax2 * 0.8), Max[0, x(c2point2_sub)]], Min[ymax2, Max[ymin2, If[(x(c2point2_sub) > Min[(xmax2 * 0.8), Max[0, x(c2point3_sub)]] ∨ abs(y(c2point2_sub) - y(c2point1)) < threshold_y2) ∧ abs((round(y(c2point2_sub) / dy) * dy) - y(c2point2_sub)) < threshold_y2, (round(y(c2point2_sub) / dy) * dy), y(c2point2_sub)]]])
t22: x(c2point2)
t23: x(c2point3)
y21: y(c2point1)
y22: y(c2point2)
y23: y(c2point3)
b2alpha: ((((y22 - y21) * t23) * ℯ^((calpha * t23))) - (((y23 - y21) * t22) * ℯ^((calpha * t22)))) / ((((ℯ^((calpha * t22)) - 1) * t23) * ℯ^((calpha * t23))) - (((ℯ^((calpha * t23)) - 1) * t22) * ℯ^((calpha * t22))))
a2alpha: (y22 - y21 - ((ℯ^((calpha * t22)) - 1) * b2alpha)) / ((t22 * ℯ^((calpha * t22))))
r2(x) = (y22 - y21) / (ℯ^((t22 * x)) - 1) - (y23 - y21) / (ℯ^((t23 * x)) - 1)
c2exp: x(Element[RemoveUndefined[{Root[r2, 1.0E-4], Root[r2, -1.0E-4]}], 1])
b2exp: (y22 - y21) / (ℯ^((c2exp * t22)) - 1)
f2alpha(x) = y21 - b2alpha + (((a2alpha * x) + b2alpha) * ℯ^((calpha * x)))
f2exp(x) = y21 - b2exp + (b2exp * ℯ^((c2exp * x)))
f2(x) = If[x ≥ 0, If[Defined[b2exp], f2exp, f2alpha]]
c2monotonic = Defined(b2exp)
y2right = If(t22 > t23, y22, y23)
LL1 = Corner(1, 1)
UR1 = Corner(1, 3)
LL2 = Corner(2, 1)
UR2 = Corner(2, 3)
xmin2 = -1.0128410914927768
y01: y11
y02: y21
xlabel2displayed: "" + (LaTeX[xlabel2]) + ""
ylabel2displayed: "" + (LaTeX[ylabel2]) + ""
y1left = If(t12 > t13, y13, y12)
y2left = If(t22 > t23, y23, y22)
initialdir1 = If(y1left > y11, 1, If[y1left < y11, -1, If[y1right > y11, -1, If[y1right < y11, 1, 0]]])
initialdir2 = If(y2left > y21, 1, If[y2left < y21, -1, If[y2right > y21, -1, If[y2right < y21, 1, 0]]])
yinf2 = If(c2exp > 1 / xmax2, If[f2(xmax2) > ymax2, Infinity, If[f2(xmax2) < ymin2, -Infinity]], If[abs(y2right - f2(xmax2)) < threshold_y1, y2right])