eik = (-1, 2)
Z_1 = (0, 0)
Z_2 = (3, 0)
a = Segment(eik, Z_1)
b = OrthogonalLine(Z_1, a)
c = Circle(Z_1, eik)
D = Intersect(c, b, 2)
d = Segment(D, Z_2)
e = Circle(Z_2, D)
f = OrthogonalLine(Z_2, d)
E = Intersect(e, f, 2)
g = Segment(eik, E)
schat = Midpoint(eik, E)
h = Segment(Z_1, D)
i = Segment(Z_2, E)