c语言程序设计三元一次方程组求解

如题所述

第1个回答推荐于2017-09-12

不知道你知不知道行列式
算法:
定义运算|a b c|
|d e f|=aei+bfg+cdh-afh-bdi-ceg
|g h i|
三元一次方程组
ax+by+cz=d
ex+fy+gz=h
ix+jy+kz=l
|a b c|
令M=|e f g|
|i j k|

|d b c|
M1=|h f g|
|l j k|

|a d c|
M2=|e h g|
|i l k|

|a b d|
M3=|e f h|
|i j l|
则x=M1/M
y=M2/M
z=M3/M
程序:
#include<stdio.h>
int fx(int a,int b,int c,
int d,int e,int f,
int g,int h,int i);
void main()
{int a,b,c,d,e,f,g,h,i,j,k,l,x,y,z,m;
printf("a,b,c,d\n");
scanf("%d,%d,%d,%d",&a,&b,&c,&d);
printf("e,f,g,h\n");
scanf("%d,%d,%d,%d",&e,&f,&g,&h);
printf("i,j,k,l\n");
scanf("%d,%d,%d,%d",&i,&j,&k,&l);
m=fx(a,b,c,e,f,g,i,j,k);
if(m==0)
{printf("x=0,y=0,z=0\n");
return ;
}
x=fx(d,b,c,h,f,g,l,j,k)/(float)m;
y=fx(a,d,c,e,h,g,i,l,k)/(float)m;
z=fx(a,b,d,e,f,h,i,j,l)/(float)m;
printf("x=%d,y=%d,z=%d\n",x,y,z);
}

int fx(int a,int b,int c,
int d,int e,int f,
int g,int h,int i)
{int result;
result=a*e*i+b*f*g+c*d*h-a*f*h-b*d*i-c*e*g;
return result;
}本回答被提问者采纳

第2个回答 2010-04-26

行列式好点啊，呵呵

相似回答

求解方程组Ax=b; 用C语言编程算法:答：float B[3]={-0.05,-0.04,0.02}; //三元一次方程组的常数项 float X[3]; //三元一次方程组的解:X1=-0.750,X2=0.000,X3=-0.417 float DETA ( ) //行列式计算 { return (A[0][0]*A[1][1]*A[2][2]+A[1][0]*A[2][1]*A[0][2]+A[2][0]*A[0][1]*A[...

解三元一次方程的c语言程序答：float a[3][3],detA,detA1,detA2,detA3;float b[3],x,y,z;int i,j;printf("3元一次方程组的格式如下：\n");printf("a[0][0]x+a[0][1]y+a[0][2]z = b[0]\n");printf("a[1][0]x+a[1][1]y+a[1][2]z = b[1]\n");printf("a[2][0]x+a[2][1]y+a...

用c语言实现高斯消去法,解三元一次方程组。求具体程序!!答：{ cout<<"Gauss 消去法得到的原方程组的解为"<<endl;for(int k=0;k<n;k++)cout<<x[k]<<" ";} int main(){ double a[MAX][MAX+1],x[MAX],sum,max,t;int n,i,j,k,max_i;cout<<"输入原方程组的阶"<<endl; cin>>n;input(a,n);for(k=0;k<n-1;k++)//选主元素 ...

c语言程序设计 三元一次方程组求解

c语言程序设计三元一次方程组求解