# Find the Centroid of a triangle by finding the intersection
# of two of the three "point to midpoint-of-side" lines.
#
# Problem posed by Lee Bradley in the Facebook "not just tiny-c
# Programming Group".
#
# This solution coded by Jim McClanahan W4JBM (Jan 2023)
#
# Figure below is for reference:
#
# Points do not need to be in any particular order as
# long as they form a triangle.
#
#             ^ (Point 2)
#            / \
#    Side 1 /   \ Side 0
#          /     \
# Point 0 +-------+ Point 1
#           Side 2
#
# Side 0 is across the Angle formed at Point 0, etc.

# cosmetic mods, user may now enter coordinates of 3 points - 7/9/23 - lrb

from sympy import symbols, Eq, solve

def line_between_points(P1, P2):
    # print(P1)
    # print(P2)
    # Calculate Slope
    y_delta = P2[1] - P1[1]
    x_delta = P2[0] - P1[0]
    slope = y_delta / x_delta

    # Y = Slope*X + Intercept
    # Intercept = Y - Slope*X
    intercept = P1[1] - slope*P1[0]

    print("Line through Point 1 and Point 2 is:",
           "y = ", slope, "* x + ", intercept)
    return(slope,intercept)

# Definition of the three original points

print("\nEnter x coordinate of Point 0")
x0 = int(input())
print("\nEnter y coordinate of Point 0")
y0 = int(input())

print("\nEnter x coordinate of Point 1")
x1 = int(input())
print("\nEnter y coordinate of Point 1")
y1 = int(input())

print("\nEnter x coordinate of Point 2")
x2 = int(input())
print("\nEnter y coordinate of Point 2")
y2 = int(input())

# x_given = [0, 4, 10]
# y_given = [0, 0, 6]

x_given = [x0, x1, x2]
y_given = [y0, y1, y2]

x_mid_side0 = (x_given[1] + x_given[2]) / 2.0
y_mid_side0 = (y_given[1] + y_given[2]) / 2.0

x_mid_side1 = (x_given[0] + x_given[2]) / 2.0
y_mid_side1 = (y_given[0] + y_given[2]) / 2.0

print("\nEvaluating Point 0 to Mid-Point of Side 0...")
slope0, intercept0 = line_between_points([x_given[0], y_given[0]], [x_mid_side0, y_mid_side0])
print("\nEvaluating Point 1 to Mid-Point of Side 1...")
slope1, intercept1 = line_between_points([x_given[1], y_given[1]], [x_mid_side1, y_mid_side1])
print()

x, y = symbols('x,y')

# Y = Slope*X + Intercept
# Y - Slope*X = Intercept
equation0 = Eq((y-slope0*x), intercept0)
# print(equation0)
equation1 = Eq((y-slope1*x), intercept1)
# print(equation1)

print("Solution for unknowns...")
print(solve((equation0, equation1), (x, y)))