三元一次方程怎麼解程序

㈠ 用c語言實現高斯消去法，解三元一次方程組。求具體程序！！

#include<iostream>

#include<cmath>

usingnamespacestd;

#defineMAX50

voidinput(doublea[MAX][MAX+1],intn)

{

cout<<"輸入原方程組的增廣矩陣"<<endl;

for(inti=0;i<n;i++)

for(intj=0;j<n+1;j++)

cin>>a[i][j];

}

voidoutput(doublex[],int稿灶n)

{

cout<<"Gauss消去法得到的原方程組的解為"<<endl;

for(intk=0;k<n;k++)

cout<<x[k]<<"";

}

intmain()

{

doublea[MAX][MAX+1],x[MAX],sum,max,t;

intn,i,j,k,max_i;

cout<<"輸入原方程組的階"<<endl;cin>>n;

input(a,n);

for(k=0;k<n-1;k++)//選主元素

{max=a[k][k];

max_i=k;

for(i=k+1;i<n;i++)

if(fabs(a[i][k])>fabs(max))

{

max=a[i][k];

max_i=i;

}

if(max==0)

break;

if(max_i!=k)//交換兩行

for(j=k;j<n+1;j++)

{

t=a[k][j];

a[k][j]=a[max_i][j];

鍵罩扮a[max_i][j]=t;

}

for(i=k+1;i<n;i++)

{

a[i][k]=a[i][k]/-a[k][k];

for(j=k+1;j<n+1;j++)

a[i][j]=a[i][j]+a[i][k]*a[k][j];

悶培}//消元

}

if(max==0)cout<<"原方程組無解"<<endl;

else

{

for(k=n-1;k>=0;k--)

{

sum=0;

for(j=k+1;j<n;j++)

sum=sum+a[k][j]*x[j];

x[k]=(a[k][n]-sum)/a[k][k];

}//回代

output(x,n);

cout<<endl;

}

return0;

}

㈡ 已知一個三元一次方程，且方程兩邊都是整數，請寫出求根的程序。

解答過程如下：7(x-2)=2x+3；解：7x-7×2=2x+3（先把7×（x-2）用乘法分配率）；7x-14=2x+3（算出7×2）；7x-14+14=2x+3+14（左右兩邊同時加14）；7x=2x+3+14（左邊只剩7x）；-2x=3+14（先把7x和2x分別看做兩個整體）；5X=17；X=17／5

擴展內容：

使方程左右兩邊相等的未知數的值，叫核塵棚做方程的解。求方程的解的過程叫做解方程。必須含有未知數等式的等式才叫方程。等式不一定是方程，方程一定是改則等式。

相關概念

1.含有未知數的等式叫方程，也可以說是含有未知兄飢數的等式是方程。

2.使等式成立的未知數的值，稱為方程的解，或方程的根。

3.解方程就是求出方程中所有未知數的值的過程。

4.方程一定是等式，等式不一定是方程。不含未知數的等式不是方程。

5.驗證：一般解方程之後，需要進行驗證。驗證就是將解得的未知數的值代入原方程，看看方程兩邊是否相等。如果相等，那麼所求得的值就是方程的解。

6.注意事項：寫「解」字，等號對齊，檢驗。

7.方程依靠等式各部分的關系，和加減乘除各部分的關系（加數+加數=和，和-其中一個加數=另一個加數，差+減數=被減數，被減數-減數=差，被減數-差=減數，因數×因數=積，積÷一個因數=另一個因數，被除數÷除數=商，被除數÷商=除數，商×除數=被除數）

㈢ c語言程序設計 三元一次方程組求解

不知道你知不知道行列式
演算法:
定義運算|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;
}

㈣ 解三元一次方程的c語言程序

#include <stdio.h>#define EPSINON 0.000001void main(){
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[2][2]z = b[2]\n");printf("請依正數次輸入方程系數:\n"悶猛);
for(i=0; i<3; i++){ for(j=0; j<3; j++) { printf("a[%d][%d]=", i, j); scanf("%f", &a[i][j]); }}printf("依次輸入方式右邊參數:\n");for(i=0; i<3; i++){ printf("b[%d]=",i); scanf("%f",&b[i]);}
detA = a[0][0]*a[1][1]*a[2][2]+a[0][1]*a[1][2]*a[2][0]+ a[0][2]*a[1][0]*a[2][1]-a[0][0]*a[1][2]*a[2][1]- a[0][1]*a[1][0]*a[2][2]-a[0][2]*a[1][1]*a[2][0];
detA1 = b[0]*a[1][1]*a[2][2]+a[0][1]*a[1][2]*b[2]+ a[0][2]*b[1]*a[2][1]-b[0]*a[1][2]*a[2][1]- a[0][1]*b[1]*a[2][2]-a[0][2]*a[1][1]*b[2];
detA2 = a[0][0]*b[1]*a[2][2]+b[0]*a[1][2]*a[2][0]+ a[0][2]*a[1][0]*b[2]-a[0][0]*a[1][2]*b[2]- b[0]*a[1][0]*a[2][2]-a[0][2]*b[1]*a[2][0];
detA3 = a[0][0]*a[1][1]*b[2]+a[0][1]*b[1]*a[2][0]+ b[0]*a[1][0]*a[2][1]-a[0][0]*b[1]*a[2][1]- a[0][1]*a[1][0]*b[2]-b[0]*a[1][1]*a[2][0];
x = detA1/detA;y = detA2/detA;z = detA3/detA;if ((detA >= -EPSINON) && (detA <= EPSINON)){ printf("方程無唯一解!\n"); }else{ printf("方程組的根為：x=%f,y=%f,z=%f\n", x, y, z); }}

網上搜的用法是：3x-y+z=3 方程系數依次是輸入3-11 方式右邊參數是輸入3

㈤ c++ 奼傝В涓夊厓涓嬈℃柟紼嬬粍

瑙ｄ笁鍏冩柟紼

鏂圭▼褰㈠紡
a1*X + b1*Y + c1*Z = d1
a2*X + b2*Y + c2*Z = d2
a3*X + b3*Y + c3*Z = d3

#include <stdio.h>

void equation2(double *x, double *y, double a1, double b1, double c1, double a2, double b2, double c2);
void equation3(double *x, double *y, double *z, double a1, double b1, double c1, double d1, double a2, double b2, double c2, double d2, double a3, double b3, double c3, double d3);

void equation2(double *x, double *y, double a1, double b1, double c1, double a2, double b2, double c2)
{
*x = (c1 * b2 - c2 * b1) / (a1 * b2 - a2 * b1);
*y = (c1 - a1 * *x) / b1;

}

void equation3(double *x, double *y, double *z, double a1, double b1, double c1, double d1, double a2, double b2, double c2, double d2, double a3, double b3, double c3, double d3)
{
equation2(x, y, a1 * c2 - a2 * c1, b1 * c2 - b2 * c1, d1 * c2 - d2 * c1, a1 * c3 - a3 * c1, b1 * c3 - b3 * c1, d1 * c3 - d3 * c1);
*z = (d1 - a1 * *x - b1 * *y) / c1;

}

int main(void)
{
double x, y, z;
double a1, b1, c1, d1;
double a2, b2, c2, d2;
double a3, b3, c3, d3;

a1 = 75982;
b1 = 56291;
c1 = 32609;
d1 = 191000900;

a2 = 43608;
b2 = 35117;
c2 = 34569;
d2 = 152065100;

a3 = 73648;
b3 = 35117;
c3 = 28039;
d3 = 161731100;

equation3(&x, &y, &z, a1, b1, c1, d1, a2, b2, c2, d2, a3, b3, c3, d3);

printf("x=%15.12f\ny=%15.12f\nz=%15.12f\n",x,y,z);
}