Question

三元一次方程组 公式 c语言

#include <stdio.h>
#include <conio.h>
#include <math.h>
main()
{
double a,b,c,d,A,B,C,D,e,f,g,h,x,y,z;
scanf("%lf%lf%lf%lf%lf%lf%lf%lf%lf%lf%lf%lf",&a,&b,&c,&d,&A,&B,&C,&D,&e,&f,&g,&h);
x=(d-(b*y)-(c*z))/a;
y=(D-((d-(b*y)-(c*z))/a*A)-(C*z))/B;
z=(h-((d-(b*y)-(c*z))/a)*e-(f*(D-((d-(b*y)-(c*z))/a*A)-(C*z))/B))/g;
printf("%lf\n%lf\n%lf",x,y,z);
}

方程组：{a*x+b*y+c*z=d
{A*x+B*y+C*z=D
{e*x+f*y+g*z=h
x=(d-(b*y)-(c*z))/a;
y=(D-((d-(b*y)-(c*z))/a*A)-(C*z))/B;
z=(h-((d-(b*y)-(c*z))/a)*e-(f*(D-((d-(b*y)-(c*z))/a*A)-(C*z))/B))/g;

公式有误！
请高手帮忙

Answer 1

不知你那公式怎么来的，建议采用线性代数的克莱姆方法来解：
为方便和完善期间，增加一个变量H,把x=,y=,z=那三行和后面那一行
替换为：
H=a*(B*g-C*f)-A*(b*g-c*f)+e*(b*C-c*B);
if(H==0) printf("次方程组无解！");
else
{x=(d*(B*g-C*f)-D*(b*g-c*f)+h*(b*C-c*B))/H;
y=(d*(A*g-C*e)-D*(a*g-c*e)+h*(a*C-c*A))/H;
z=(d*(A*f-B*e)-D*(a*f-b*e)+h*(a*B-b*A))/H;
printf("%lf\n%lf\n%lf",x,y,z);}

试验一下，不行再问。其实你的变量设置不大好，为了尽量保持原貌，未做改动。

追问

答案错了呀！

追答

没关系，说明程序已经基本正确，仔细检查代码，尤其是大小写。我已经说过，你的变量设置不好，大小写很容易错。调试就行了。
把你的程序发一份给我，有时间给你看一下。

Answer 2

赋值表达式的右值必须确定才能进行赋值，你把X,Y,Z的方程化简一下，让表达式的右边没有未知数。