X(14961), labeled as O1
Sketch by Peter Moses, October 25, 2017

A = (7.492727272727278, -4.016363636363641)
B = (6.529090909090913, -7.870909090909099)
C = (11.874545454545464, -7.925454545454554)
a = Segment(B, C)
b = Segment(C, A)
c = Segment(A, B)
L1 = TriangleCurve(A, B, C, (B * (a^(4) - ((2 * a^(2)) * b^(2)) + b^(4) + ((2 * a^(2)) * c^(2)) + ((2 * b^(2)) * c^(2)) - (3 * c^(4)))) + (A * (a^(4) + ((2 * a^(2)) * b^(2)) - (3 * b^(4)) - ((2 * a^(2)) * c^(2)) + ((2 * b^(2)) * c^(2)) + c^(4))) + (((-3 * a^(4)) + ((2 * a^(2)) * b^(2)) + b^(4) + ((2 * a^(2)) * c^(2)) - ((2 * b^(2)) * c^(2)) + c^(4)) * C) = 0)
L2 = TriangleCurve(A, B, C, (A * (a^(4) - ((2 * a^(2)) * b^(2)) + b^(4) + ((2 * a^(2)) * c^(2)) + ((2 * b^(2)) * c^(2)) - (3 * c^(4)))) + (B * ((-3 * a^(4)) + ((2 * a^(2)) * b^(2)) + b^(4) + ((2 * a^(2)) * c^(2)) - ((2 * b^(2)) * c^(2)) + c^(4))) + ((a^(4) + ((2 * a^(2)) * b^(2)) - (3 * b^(4)) - ((2 * a^(2)) * c^(2)) + ((2 * b^(2)) * c^(2)) + c^(4)) * C) = 0)
Brocard = TriangleCurve(A, B, C, (((((-a^(2)) * B) * (a - c)) * c^(2)) * (a + c)) + ((((A * b^(2)) * (b - c)) * c^(2)) * (b + c)) + ((((a^(2) * (a - b)) * b^(2)) * (a + b)) * C) = 0)
O = TriangleCenter(A, B, C, 3)
O1 = Barycenter({A, B, C}, {((a^(2) * (a^(2) - b^(2) - c^(2))) * ((a^(4) * b^(2)) - b^(6) + (a^(4) * c^(2)) - (((2 * a^(2)) * b^(2)) * c^(2)) + (b^(4) * c^(2)) + (b^(2) * c^(4)) - c^(6))), ((b^(2) * ((-a^(2)) + b^(2) - c^(2))) * ((-a^(6)) + (a^(2) * b^(4)) + (a^(4) * c^(2)) - (((2 * a^(2)) * b^(2)) * c^(2)) + (b^(4) * c^(2)) + (a^(2) * c^(4)) - c^(6))), ((c^(2) * ((-a^(2)) - b^(2) + c^(2))) * ((-a^(6)) + (a^(4) * b^(2)) + (a^(2) * b^(4)) - b^(6) - (((2 * a^(2)) * b^(2)) * c^(2)) + (a^(2) * c^(4)) + (b^(2) * c^(4))))})
L0 = TriangleCurve(A, B, C, (A * ((-3 * a^(4)) + ((2 * a^(2)) * b^(2)) + b^(4) + ((2 * a^(2)) * c^(2)) - ((2 * b^(2)) * c^(2)) + c^(4))) + (B * (a^(4) + ((2 * a^(2)) * b^(2)) - (3 * b^(4)) - ((2 * a^(2)) * c^(2)) + ((2 * b^(2)) * c^(2)) + c^(4))) + ((a^(4) - ((2 * a^(2)) * b^(2)) + b^(4) + ((2 * a^(2)) * c^(2)) + ((2 * b^(2)) * c^(2)) - (3 * c^(4))) * C) = 0)
U = Barycenter({A, B, C}, {((a^(2) * (a^(2) - b^(2) - c^(2))) * (((3 * a^(4)) * b^(2)) - ((2 * a^(2)) * b^(4)) - b^(6) - (a^(4) * c^(2)) - (((2 * a^(2)) * b^(2)) * c^(2)) + ((3 * b^(4)) * c^(2)) + ((2 * a^(2)) * c^(4)) - (b^(2) * c^(4)) - c^(6))), (((-b^(2)) * ((-a^(2)) + b^(2) - c^(2))) * (a^(6) - ((2 * a^(4)) * b^(2)) + (a^(2) * b^(4)) + (a^(4) * c^(2)) + (((2 * a^(2)) * b^(2)) * c^(2)) - ((3 * b^(4)) * c^(2)) - ((3 * a^(2)) * c^(4)) + ((2 * b^(2)) * c^(4)) + c^(6))), ((c^(2) * ((-a^(2)) - b^(2) + c^(2))) * ((-a^(6)) + ((3 * a^(4)) * b^(2)) - (a^(2) * b^(4)) - b^(6) - ((2 * a^(4)) * c^(2)) - (((2 * a^(2)) * b^(2)) * c^(2)) + ((2 * b^(4)) * c^(2)) + ((3 * a^(2)) * c^(4)) - (b^(2) * c^(4))))})
P = Barycenter({A, B, C}, {((a^(2) * (a^(2) - b^(2) - c^(2))) * ((a^(4) * b^(2)) - ((2 * a^(2)) * b^(4)) + b^(6) - ((3 * a^(4)) * c^(2)) + (((2 * a^(2)) * b^(2)) * c^(2)) + (b^(4) * c^(2)) + ((2 * a^(2)) * c^(4)) - ((3 * b^(2)) * c^(4)) + c^(6))), (((-b^(2)) * ((-a^(2)) + b^(2) - c^(2))) * ((-a^(6)) - ((2 * a^(4)) * b^(2)) + ((3 * a^(2)) * b^(4)) + ((3 * a^(4)) * c^(2)) - (((2 * a^(2)) * b^(2)) * c^(2)) - (b^(4) * c^(2)) - (a^(2) * c^(4)) + ((2 * b^(2)) * c^(4)) - c^(6))), ((c^(2) * ((-a^(2)) - b^(2) + c^(2))) * (a^(6) + (a^(4) * b^(2)) - ((3 * a^(2)) * b^(4)) + b^(6) - ((2 * a^(4)) * c^(2)) + (((2 * a^(2)) * b^(2)) * c^(2)) + ((2 * b^(4)) * c^(2)) + (a^(2) * c^(4)) - ((3 * b^(2)) * c^(4))))})
E = Midpoint(U, P)
d = Circle(E, U)