from itertools import combinations as C
from math import sqrt,fabs
from decimal import Decimal
def cross_v(p1,p2):
	return (p1[1]*p2[2]-p1[2]*p2[1],p1[2]*p2[0]-p1[0]*p2[2],p1[0]*p2[1]-p1[1]*p2[0])
def inner_p(p1,p2):
	return sum(p1[i]*p2[i] for i in range(3))
def vec_s(p1,p2):
	return (p1[0]-p2[0],p1[1]-p2[1],p1[2]-p2[2])
def dps(q1,q2,q3):
	n_v = cross_v(vec_s(q2,q1),vec_s(q3,q1))
	return fabs(inner_p(vec_s(P,q1),n_v)) / sqrt(sum(c**2 for c in n_v))
N = input()
f = lambda x: float(x)
P = map(f,raw_input().split())
Q = [map(f,raw_input().split()) for i in range(N)]
print Decimal(sum([dps(*q) for q in C(Q,3)]))