空间四面体体积公式

int T,n,m;
int va[N];

db volume(db x1,db y1,db z1,db x2,db y2,db z2,db x3,db y3,db z3,db x4,db y4,db z4);

signed main()
{
	IOS;
	double x1,y1,z1,x2,y2,z2,x3,y3,z3,x4,y4,z4;
    cin>>x1>>y1>>z1;
    cin>>x2>>y2>>z2;
    cin>>x3>>y3>>z3;
    cin>>x4>>y4>>z4;
    double s = volume(x1,y1,z1,x2,y2,z2,x3,y3,z3,x4,y4,z4);
    printf("%.10lf\n",s);
}

db volume(db x1,db y1,db z1,db x2,db y2,db z2,db x3,db y3,db z3,db x4,db y4,db z4)
{
    db A1, A2, A3, A4, value;
    A1 = x2 * (y3 * z4 - y4 * z3) - x3 * (y2 * z4 - y4 * z2) + x4 * (y2 * z3 - y3 * z2);
    A2 = x1 * (y3 * z4 - y4 * z3) - x3 * (y1 * z4 - y4 * z1) + x4 * (y1 * z3 - y3 * z1);
    A3 = x1 * (y2 * z4 - y4 * z2) - x2 * (y1 * z4 - y4 * z1) + x4 * (y1 * z2 - y2 * z1);
    A4 = x1 * (y2 * z3 - y3 * z2) - x2 * (y1 * z3 - y3 * z1) + x3 * (y1 * z2 - y2 * z1);
    value = (A1 - A2 + A3 - A4) / 6;
    return fabs(value);
}
上一篇:最小二乘法以及拓展的最小二乘法辨识电机参数


下一篇:Java中的Math类