C++3点求面，两角度求面算法

已知3点坐标，求平面ax+by+cz+d=0;

得出的abcd均在等式的一侧，如上表达式

void get_panel(TPoint p1,TPoint p2,TPoint p3,double &a,double &b,double &c,double &d){

a = ( (p2.y-p1.y)*(p3.z-p1.z)-(p2.z-p1.z)*(p3.y-p1.y) );

b = ( (p2.z-p1.z)*(p3.x-p1.x)-(p2.x-p1.x)*(p3.z-p1.z) );

c = ( (p2.x-p1.x)*(p3.y-p1.y)-(p2.y-p1.y)*(p3.x-p1.x) );

d = ( 0-(a*p1.x+b*p1.y+c*p1.z) );

}

目前认定是所有坐标系通用//右手坐标系
tmp为面法向量
Dip_Angle 是倾角    即平面与XOZ面的夹角  范围为-PI/2 ~PI/2
Inclination是倾向    即平面与XOY面的夹角  范围为-PI/2 ~PI/2
void InfoTo_XYZ(void)
{  

    a = cos(PI / 2 - Dip_Angle)* cos(PI / 2 - Inclination);
    b = cos(PI / 2 - Dip_Angle)* sin(PI / 2 - Inclination);
    c = sin(PI / 2 - Dip_Angle);
    d  =  (  0 - (a * p1.x + b * p1.y + c * p1.z) );

}