Ha 2013

A = (-1.5, -0.7)
B = (0.92, 2.12)
C = (2.58, -0.64)
D = (5.58, 1.02)
f = Polynomial({A, B, C, D})
P = Point(f)
a = Tangent(P, f)
W = TurningPoint(f)
P' = Mirror(P, W)
a' = Mirror(a, W)
E, G = Intersect(f, a')
H = Intersect(f, a, 2)
b = OrthogonalLine(P, xAxis)
c = OrthogonalLine(P', xAxis)
d = OrthogonalLine(H, xAxis)
e = OrthogonalLine(E, xAxis)
g = Line(W, a')
h = OrthogonalLine(W, xAxis)
I = Intersect(a, e)
J = Intersect(a', d)
Vieleck1, e_1, i, h_1, j = Polygon(E, I, H, J)