A = (5.185544821032594, 4.0613763066202)
B: (0, 0)
c = Segment(A, B)
C: (197 / 13, 0)
b = Segment(C, A)
a = Segment(B, C)
cubicEqn = TriangleCurve(A, B, C, 0 = (((a^(2) * A^(3)) * b^(4)) * c^(4)) + (((a^(4) * b^(2)) * B^(3)) * c^(4)) - ((((a^(2) * A) * B^(2)) * c^(4)) * (a^(4) - ((2 * a^(2)) * b^(2)) + (2 * b^(4)) - ((2 * a^(2)) * c^(2)) + c^(4))) - ((((A^(2) * b^(2)) * B) * c^(4)) * ((2 * a^(4)) - ((2 * a^(2)) * b^(2)) + b^(4) - ((2 * b^(2)) * c^(2)) + c^(4))) - ((((A^(2) * b^(4)) * c^(2)) * ((2 * a^(4)) + b^(4) - ((2 * a^(2)) * c^(2)) - ((2 * b^(2)) * c^(2)) + c^(4))) * C) - ((((a^(4) * B^(2)) * c^(2)) * (a^(4) + (2 * b^(4)) - ((2 * a^(2)) * c^(2)) - ((2 * b^(2)) * c^(2)) + c^(4))) * C) - ((((((a^(2) * A) * b^(2)) * B) * c^(2)) * ((5 * a^(4)) - ((4 * a^(2)) * b^(2)) + (5 * b^(4)) - ((4 * a^(2)) * c^(2)) - ((4 * b^(2)) * c^(2)) + (5 * c^(4)))) * C) - ((((a^(2) * A) * b^(4)) * (a^(4) - ((2 * a^(2)) * b^(2)) + b^(4) - ((2 * a^(2)) * c^(2)) + (2 * c^(4)))) * C^(2)) - ((((a^(4) * b^(2)) * B) * (a^(4) - ((2 * a^(2)) * b^(2)) + b^(4) - ((2 * b^(2)) * c^(2)) + (2 * c^(4)))) * C^(2)) + (((a^(4) * b^(4)) * c^(2)) * C^(3)))
MoveMe = Point(cubicEqn)
AB = Line(A, B)
BC = Line(B, C)
CA = Line(C, A)
barysP = Element(Transpose[ReducedRowEchelonForm[{{x(A), x(B), x(C), x(MoveMe)}, {y(A), y(B), y(C), y(MoveMe)}, {1, 1, 1, 1}}]], 4)
CIP = Barycenter({A, B, C}, {((((barysP(1) * (barysP(1) + barysP(2) + barysP(3))) * a^(2)) * b^(2)) * c^(2)) + ((a^(2) * ((-a^(2)) + b^(2) + c^(2))) * (((barysP(2) * barysP(3)) * a^(2)) + ((barysP(1) * barysP(3)) * b^(2)) + ((barysP(1) * barysP(2)) * c^(2)))), ((((barysP(2) * (barysP(1) + barysP(2) + barysP(3))) * a^(2)) * b^(2)) * c^(2)) + ((b^(2) * (a^(2) - b^(2) + c^(2))) * (((barysP(2) * barysP(3)) * a^(2)) + ((barysP(1) * barysP(3)) * b^(2)) + ((barysP(1) * barysP(2)) * c^(2)))), ((((barysP(3) * (barysP(1) + barysP(2) + barysP(3))) * a^(2)) * b^(2)) * c^(2)) + ((c^(2) * (a^(2) + b^(2) - c^(2))) * (((barysP(2) * barysP(3)) * a^(2)) + ((barysP(1) * barysP(3)) * b^(2)) + ((barysP(1) * barysP(2)) * c^(2))))})
x3 = TriangleCenter(A, B, C, 3)
x6 = TriangleCenter(A, B, C, 6)
x1316 = TriangleCenter(A, B, C, 1316)
BrocardCircle = Circle(x3, x6, x1316)