用C语言编写程序：N元一次方程组的解

1、./*2013 年 4 月 9 日 12:23:20功能：已知有一三元一次方程组如下,方程系数由键盘输入，求其3 个根。*/# include <stdio.h>int main(void)int i, j, m, n, r, k=0, t; /i,j,m,k,t作为循环计数变量；n 表示有 n 个未知数，即n 元； r 表示有 r 个方程组成方程组。/ n = r = 3;printf(" 请输入未知数个数 n 和方程的个数 :n");scanf("%d %d", &n, &r); /n表示有 n 个未知数，即n 元； r

2、表示有 r 个方程组成方程组。float a1120;/用二维数组存储方程组参数。11 和 100 课改动。这里的程序限 10 元 20 个方程的方程组。float b200;/原来存放之间变量。printf(" 请输入方程组参数n");for(j=0; j<n; j+)for(i=0; i<r+1; i+)scanf("%f", &aji);/*for(j=0; j<n; j+)/此段程序作为检测该程序正确性的试数作用。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(" 方程有无穷解 n");if(r=n)printf(" 方程有唯一解 n");for(t=0; t<n-1; t+)/循环次数：当t=0 时，仅将a10,a20,a21的值转换为0； t=1 时，才将其转换成最简型阶梯矩阵。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+) /因为一行的每个元素均要进行运算等式两边才会相等。aji -= ami*bk;bk=0;/此处是必不可少的。因为内存中会有垃圾值，必须要处理。k+;for(j=0; j<n; j+)/为了将其转换成单位矩阵，见线性代数第98 页。ajn /= ajj;ajj /=ajj;.for(j=0; j<r; j+)printf("x%d = %fn", j, ajn); /输出最终结果，即未知数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其正确结果是： x=1,y=3,z=2在 VC+6.0 中的输出结果是：请输入未知数个数 n 和方程的个数 :.3 3请输入方程组参数22-161-24357128方程有唯一解x0 = 1.000000x1 = 3.000000x2 = 2.0000001.000000 -0.000000 0.000000 1.000