#include<bits/stdc++.h>
#define LL long long
#define LD long double
using namespace std;
LD X[5],Y[5],Z[5];
template<class T>void read(T &x)
{
	x=0;int f=0;char ch=getchar();
	while(ch<'0' || ch>'9')f|=(ch=='-'),ch=getchar();
	while(ch>='0' && ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
	x=f? -x:x;return ;
}
inline LD calc(LD x,LD y,LD z)
{
	//printf(" calc %Lf %Lf %Lf\n",x,y,z);
	return x*y*z/6;
}
inline LD pw3(LD x)
{
	return x*x*x;
}
inline void print()
{
	//printf(" (%Lf,%Lf,%Lf) (%Lf,%Lf,%Lf) (%Lf,%Lf,%Lf) (%Lf,%Lf,%Lf)\n",X[1],Y[1],Z[1],X[2],Y[2],Z[2],X[3],Y[3],Z[3],X[4],Y[4],Z[4]);
}
int main()
{
	//freopen("data.out","r",stdin);
	//freopen("my.out","w",stdout);
	for(;~scanf("%Lf%Lf%Lf",&X[1],&Y[1],&Z[1]);)
	{
		for(int i=2;i<=4;++i)scanf("%Lf%Lf%Lf",&X[i],&Y[i],&Z[i]);
		if(X[2]<X[1])for(int i=1;i<=4;++i)X[i]=-X[i];
		if(Y[2]<Y[1])for(int i=1;i<=4;++i)Y[i]=-Y[i];
		if(Z[2]<Z[1])for(int i=1;i<=4;++i)Z[i]=-Z[i];
		if(X[3]>X[4])swap(X[3],X[4]);
		if(Y[3]>Y[4])swap(Y[3],Y[4]);
		if(Z[3]>Z[4])swap(Z[3],Z[4]);
		print();
		for(int i=4;i;--i)X[i]-=X[1],Y[i]-=Y[1],Z[i]-=Z[1];
		if(X[3]>=X[2] || Y[3]>=Y[2] || Z[3]>=Z[2])
		{
			puts("0");
			continue;
		}
		if(X[3]>X[1])
		{
			Y[2]*=(X[2]-X[3])/X[2];
			Z[2]*=(X[2]-X[3])/X[2];
			X[1]=X[3];
		}
		for(int i=4;i;--i)X[i]-=X[1],Y[i]-=Y[1],Z[i]-=Z[1];
		print();
		if(Y[3]>Y[1])
		{
			X[2]*=(Y[2]-Y[3])/Y[2];
			Z[2]*=(Y[2]-Y[3])/Y[2];
			Y[1]=Y[3];
		}
		for(int i=4;i;--i)X[i]-=X[1],Y[i]-=Y[1],Z[i]-=Z[1];
		print();
		if(Z[3]>Z[1])
		{
			X[2]*=(Z[2]-Z[3])/Z[2];
			Y[2]*=(Z[2]-Z[3])/Z[2];
			Z[1]=Z[3];
		}
		for(int i=4;i;--i)X[i]-=X[1],Y[i]-=Y[1],Z[i]-=Z[1];
		print();
		X[4]=min(X[4],X[2]);
		Y[4]=min(Y[4],Y[2]);
		Z[4]=min(Z[4],Z[2]);
		if(X[1]>=X[4] || Y[1]>=Y[4] || Z[1]>=Z[4])
		{
			puts("0");
			continue;
		}
		for(int i=4;i;--i)X[i]-=X[1],Y[i]-=Y[1],Z[i]-=Z[1];
		print();
		LD w=X[4]/X[2]+Y[4]/Y[2]+Z[4]/Z[2];
		LD tmp=0,res;
		if(w>=1)
		{
			//printf(" %Lf %Lf %Lf\n",pw3((X[2]-X[4])/X[2]),pw3((Y[2]-Y[4])/Y[2]),pw3((Z[2]-Z[4])/Z[2]));
			LD x0=(1-Y[4]/Y[2]-Z[4]/Z[2])*X[2],y0=(1-X[4]/X[2]-Z[4]/Z[2])*Y[2],z0=(1-X[4]/X[2]-Y[4]/Y[2])*Z[2],za,zb,zc;
			LD k=-Y[2]/X[2];
			//printf(" x0:%Lf y0:%Lf z0:%Lf\n",x0,y0,z0);
			res=calc(X[2],Y[2],Z[2]);
			tmp+=pw3((Z[2]-Z[4])/Z[2]);
			//printf(" +%Lf\n",(Z[2]-Z[4])/Z[2]);
			za=(1-X[4]/X[2])*Z[2];
			if(za>=Z[4])tmp-=pw3((za-Z[4])/Z[2]);
			//if(za>=Z[4])printf(" -%Lf\n",(za-Z[4])/Z[2]);
			zb=(1-Y[4]/Y[2])*Z[2];
			if(zb>=Z[4])tmp-=pw3((zb-Z[4])/Z[2]);
			//if(zb>=Z[4])printf(" -%Lf\n",(zb-Z[4])/Z[2]);
			tmp+=pw3((X[2]-X[4])/X[2]);
			//printf(" +%Lf\n",(X[2]-X[4])/X[2]);
			tmp+=pw3((Y[2]-Y[4])/Y[2]);
			//printf(" +%Lf\n",(Y[2]-Y[4])/Y[2]);
			//assert(fabs(((X[2]-(X[4]-Y[4]/k))/X[2])-((Y[2]-(Y[4]-k*X[4]))/Y[2]))<1e-9);
			//printf(" (%Lf)\n",X[4]-Y[4]/k);
			zc=max(X[2]-(X[4]-Y[4]/k),(LD)0)/X[2]*Z[2];
			tmp-=pw3(zc/Z[2]);
			//printf(" -%Lf\n",zc/Z[2]);
			//printf("tmp:%Lf\n",tmp);
			tmp=1-tmp;
			res*=tmp;
			//printf("res:%Lf\n",res);
			if(fabs(w-1)<=1e-9)assert(fabs(res-X[4]*Y[4]*Z[4])<=1e-9);
		}
		else res=X[4]*Y[4]*Z[4];
		printf("%.9Lf\n",res);
		/*if(fabs(w-1)<=1e-9)cerr<<(X[4]-X[1])*(Y[4]-Y[1])*(Z[4]-Z[1])-
			calc(X[4]-(1-Y[4]/Y[2]-Z[4]/Z[2])*X[2],
			Y[4]-(1-X[4]/X[2]-Z[4]/Z[2])*Y[2],
			Z[4]-(1-X[4]/X[2]-Y[4]/Y[2])*Z[2])
		<<" = "<<(X[4]-X[1])*(Y[4]-Y[1])*(Z[4]-Z[1])<<'\n';*/
	}
	return 0;
}