用C語(yǔ)言編寫程序：N元一次方程組的解

1、./*2013 年 4 月 9 日 12:23:20功能：已知有一三元一次方程組如下,方程系數(shù)由鍵盤輸入，求其3 個(gè)根。*/# include <stdio.h>int main(void)int i, j, m, n, r, k=0, t; /i,j,m,k,t作為循環(huán)計(jì)數(shù)變量；n 表示有 n 個(gè)未知數(shù)，即n 元； r 表示有 r 個(gè)方程組成方程組。/ n = r = 3;printf(" 請(qǐng)輸入未知數(shù)個(gè)數(shù) n 和方程的個(gè)數(shù) :n");scanf("%d %d", &n, &r); /n表示有 n 個(gè)未知數(shù)，即n 元； r

2、表示有 r 個(gè)方程組成方程組。float a1120;/用二維數(shù)組存儲(chǔ)方程組參數(shù)。11 和 100 課改動(dòng)。這里的程序限 10 元 20 個(gè)方程的方程組。float b200;/原來(lái)存放之間變量。printf(" 請(qǐng)輸入方程組參數(shù)n");for(j=0; j<n; j+)for(i=0; i<r+1; i+)scanf("%f", &aji);/*for(j=0; j<n; j+)/此段程序作為檢測(cè)該程序正確性的試數(shù)作用。for(i=0; i<n+1; i+)printf("%f ", aji);pri

3、ntf("n");.float a1020=2,2,-1,6,1,-2,4,3,5,7,1,28;*/if(r<n)printf(" 方程有無(wú)窮解 n");if(r=n)printf(" 方程有唯一解 n");for(t=0; t<n-1; t+)/循環(huán)次數(shù)：當(dāng)t=0 時(shí)，僅將a10,a20,a21的值轉(zhuǎn)換為0； t=1 時(shí)，才將其轉(zhuǎn)換成最簡(jiǎn)型階梯矩陣。for(m=0; m<n; m+)/m表示第 m 列for(j=0; j<n; j+)/j表示第 j 行if(j!=m)bk = ajm/amm;for(i=0

4、; i<n+1; i+) /因?yàn)橐恍械拿總€(gè)元素均要進(jìn)行運(yùn)算等式兩邊才會(huì)相等。aji -= ami*bk;bk=0;/此處是必不可少的。因?yàn)閮?nèi)存中會(huì)有垃圾值，必須要處理。k+;for(j=0; j<n; j+)/為了將其轉(zhuǎn)換成單位矩陣，見(jiàn)線性代數(shù)第98 頁(yè)。ajn /= ajj;ajj /=ajj;.for(j=0; j<r; j+)printf("x%d = %fn", j, ajn); /輸出最終結(jié)果，即未知數(shù)x,y,z.的值。for(j=0; j<n; j+)for(i=0; i<n+1; i+)printf("%f ", aji);printf("n");return 0;/*輸入為： 22-161-24357128即方程組為： 2x+2y- z=6x-2y+4z=35x+7y+ z=28其正確結(jié)果是： x=1,y=3,z=2在 VC+6.0 中的輸出結(jié)果是：請(qǐng)輸入未知數(shù)個(gè)數(shù) n 和方程的個(gè)數(shù) :.3 3請(qǐng)輸入方程組參數(shù)22-161-24357128方程有唯一解x0 = 1.000000x1 = 3.000000x2 = 2.0000001.000000 -0.000000 0.000000 1.000