b)

O = (21.0854, -1.6314)
C = (29.8942, -1.6556)
F_2 = Vector(O, C)
a = OrthogonalLine(O, F_2)
D = Point(a)
F_1 = Vector(O, D)
β = Angle(C, O, D)