Epsilon and Delta Limits
Author: Scott Farrar

At x=c, can you always find a small enough
delta-neighborhood so that the outputs are 
in an epsilon-neighborhood around L?

A = (2, 2)
g = Function(x^(2) - x(A), 0, x(A))
B: (x(A), y(A) + 1)
h = Function(x + y(B) - 2, x(B), 10)
C: (x(A), 0)
L: (0, y(A))
a = Segment(A, C)
b = Segment(L, A)
δ = 0.42 (slider from 0 to 1, step 0.01)
ε = 1.54 (slider from 0 to 2, step 0.01)
D: (x(C) - δ, 0)
E: (x(C) + δ, 0)
c = Segment(D, E)
d = OrthogonalLine(D, xAxis)
e = OrthogonalLine(E, xAxis)
F = Intersect(g, d)
G = Intersect(h, e)
i = OrthogonalLine(F, d)
j = OrthogonalLine(G, e)
H = Intersect(j, yAxis)
I = Intersect(i, yAxis)
k = Segment(H, G)
l = Segment(I, F)
m = Segment(G, E)
n = Segment(F, D)
J: (0, y(L) + ε)
K: (0, y(L) - ε)
p = Segment(J, K)
q = OrthogonalLine(B, yAxis)
M = Intersect(yAxis, q)
poly1, m_1, h_1, g_1, b_1 = Polygon(M, H, G, B)
poly2, l_1, a_1, f, i_1 = Polygon(L, A, F, I)
poly3, d_1, c_1, a_2, f_1 = Polygon(D, C, A, F)
poly4, b_2, c_2, e_1, g_2 = Polygon(B, C, E, G)