Точка Ферма-Торричелли

from math import sqrt, atan, cos
 
def angle(xa, ya, xb, yb, xc, yc):
    
    ac = (xc-xa, yc-ya)
    bc = (xc-xb, yc-yb)
    dist = lambda x: sqrt(x[0]**2 + x[1]**2)
    cc = (ac[0]*bc[0] + ac[1]*bc[1]) / (dist(ac)*dist(bc))
 
    return cc < -0.5, xc, yc
    
 
def hab(xa, ya, xb, yb, xc, yc):
    
    xab, yab = (xa+xb)/2, (ya+yb)/2
    d = sqrt(3*((xa-xb)**2 + (ya-yb)**2))/2
    dk = lambda f: 1 - 2*(f<0)
    if xa==xb:
        return [xa + d*dk(xa), yab]
    if ya == yb:
        return [xab, ya + d*dk(ya)]
 
    kab = (ya-yb) / (xa-xb)
    kh = -1 / kab
    bh = yab - kh * xab
 
    dx = d * cos(atan(kh))
 
    x1 = xab + dx
    y1 = kh * x1 + bh
    x2 = xab - dx
    y2 = kh * x2 + bh
 
    if (x1-xc)**2 + (y1-yc)**2 < (x2-xc)**2 + (y2-yc)**2:
        x1, y1 = x2, y2
    return [x1, y1]
 
 
def point_ferma(p1, p2):
    xa, ya, xb, yb = p1
    xda, yda, xdb, ydb = p2
    
    aA = yda - ya
    aB = xa - xda
    aC = ya*xda - yda*xa
    
    bA = ydb - yb
    bB = xb - xdb
    bC = yb*xdb - ydb*xb
    
    dd = aA*bB - aB*bA
    dx = bC*aB - aC*bB 
    dy = aC*bA - bC*aA
    
    return dx/dd, dy/dd
 
 
p = list(map(float, input().split()))
for i in range(3):
    p = p[2:] + p[:2]
    cond, *point = angle(*p)
    if cond:
        break
if not cond:
    dp = []
    for i in range(3):
        p = p[2:] + p[:2]
        dp.extend(hab(*p))
 
    point = point_ferma(p[:4], dp[:4])
print(*map(lambda x: f'{x:.6f}', point))