已知如下三元一次方程组的系数值，求该方程组的解。
        a1x+b1y+c1z=d1
        a2x+b2y+c2z=d2
        a3x+b3y+c3z=d3
程序代码：
#include<stdio.h>
#include<math.h>
void main()
{
   double a1,b1,c1,d1,a2,b2,c2,d2,a3,b3,c3,d3,x,y,z;/*定义十五个双精度变量*/
   printf (" input a1,b1,c1,d1,a2,b2,c2,d2,a3,b3,c3 and d3:");/*给出输入提示*/
   scanf ("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f",&a1,&b1,&c1,&d1,&a2,&b2,&c2,&d2,&a3,&b3,&c3,&d3);
   x=((b2*d1-b1*d2)*(b3*c2-b2*c3)-(b3*d2-b2*d3)*(b2*c1-b1*c2))/((a1*b2-a2*b1)*(b3*c2-b2*c3)-(a2*b3-a3*b2)*(b2*c1-b1*c2));
   y=((a2*d1-a1*d2)*(a3*c2-a2*c3)-(a3*d2-a2*d3)*(a2*c1-a1*c2))/((a2*b1-a1*b2)*(a3*c2-a2*c3)-(a3*b2-a2*b3)*(a2*c1-a1*c2));
   z=((a2*d1-a1*d2)*(a3*b2-a2*b3)-(a3*d2-a2*d3)*(a2*b1-a1*b2))/((a2*c1-a1*c2)*(a3*b2-a2*b3)-(a3*c2-a2*c3)*(a2*b1-a1*b2));
   printf("x=%f\n,y=%f\n,z=%f\n",x,y,z);

}
程序本身没有错误，可能算法不正确。。。。。。。。请高手指教一下，谢谢

考虑不怎么完整，如果除数为零呢？和他们有无穷解咋办呢？
还是用矩阵解决这样的方程比较好！！

但是C++上我们还没学到矩阵，同样谢谢你了！