A = (3, 3)
B = (12, 1)
C = (2, -5)
t1, c, a, b = Polygon(A, B, C)
f = Line(B, C)
P = Point(f)
g = Line(A, B)
h = Line(A, C)
i = OrthogonalLine(P, f)
D = Intersect(g, i)
E = Intersect(i, h)
d = Circle(A, B, C)
e = Circle(A, D, E)
l = Tangent(A, d)
m = Tangent(A, e)
F = Intersect(d, e, 2)
j, n = Tangent(F, d)
k, p = Tangent(F, e)
γ = Angle(j, k)
ε = Angle(m, l)
α = Angle(h, g)
δ = Angle(A, B, C)
ζ = Angle(B, C, A)
O1 = Intersect(m, k)
O2 = Intersect(l, j)
AO1 = Vector(A, O1)
AO2 = Vector(A, O2)
ps: (x(AO1) * x(AO2)) + (y(AO1) * y(AO2))