#include<bits/stdc++.h>
#define ll long long
using namespace std;
struct ve
{
	ll x,y,z;
	ll operator*(ve b)const
	{
		return x*b.x+y*b.y+z*b.z;
	}
	ve operator|(ve b)const
	{
		return (ve){y*b.z-z*b.y,
		z*b.x-x*b.z,
		x*b.y-y*b.x
		};
	}
	ve operator-(ve b)const
	{
		return (ve){x-b.x,y-b.y,z-b.z};
	}
	double len()
	{
		return sqrt(x*x+y*y+z*z);
	}
}a[100];
int main()
{
	int i,j,k,l,n,m;
	ll x,y,z;
	scanf("%d",&n);
	for(i=0;i<n;i++)
	{
		scanf("%d%d%d",&x,&y,&z);
		a[i]=(ve){x,y,z};
	}
	if(n<4)
	{
		printf("0.00\n");
		return 0;
	}
	double ans=1e10;
	for(i=0;i<n;i++)
	{
		for(j=i+1;j<n;j++)
		{
			for(k=j+1;k<n;k++)
			{
				ve norm=(a[i]-a[j])|(a[i]-a[k]);
				if(norm.len()<1e-8)
				{
					continue;
				}
				int po=0,ne=0;
				double ma=0;
				for(l=0;l<n;l++)
				{
					x=norm*(a[i]-a[l]);
					if(x>=0)po++;
					if(x<=0)ne++;
					ma=max(ma,abs(x)*1.0/norm.len());
				}
				if(po==n||ne==n)
				{
					ans=min(ans,ma);
				}
			}
		}
	}
	for(i=0;i<n;i++)
	{
		for(j=i+1;j<n;j++)
		{
			for(k=j+1;k<n;k++)
			{
				for(l=k+1;l<n;l++)
				{
					ve norm = (a[i]-a[j])|(a[k]-a[l]);
					if(norm.len()<1e-8)
					{
						continue;
					}
					x=0;
					for(int p=0;p<n;p++)
					{
						if((norm*(a[i]-a[p])>=0&&norm*(a[k]-a[p])<=0)||(norm*(a[i]-a[p])<=0&&norm*(a[k]-a[p])>=0))
						{
							x++;
						}
					}
					if(x==n)
						ans=min(ans,abs((a[i]-a[k])*norm)/norm.len());
				}
			}
		}
	}
	printf("%.10lf",ans);
	return 0;
}