求n阶满秩方阵的阶梯形

版权声明：由陈庚汝编辑 https://blog.csdn.net/ChenGengru/article/details/84593115

#include<stdio.h>
#include<stdlib.h>
int main()
{
	int n, i, j, k;
	
	printf("Dimension: ");
    scanf("%d", &n);   // 得到矩阵空间的维度 // 
	
	int iArray[4][4];  // 定义n维矩阵空间 // 

    for (i = 1; i <= n; i++)  //赋值矩阵// 
	{
		for (j = 1; j <= n; j++)
		{
			printf("iArray[%d][%d] = ", i, j);
			scanf("%d", &iArray[i-1][j-1]);
		}
		printf("\n");
	}
	
    for (k = 0; k < n; k++)
	{
	    for (i = k; i < n; i++)
	    {
	    	if (iArray[i][k] != 0)
	    	{
		    for (j = 0; j < n; j++)	
		    	{
			        iArray[i][j] = iArray[i][j]*1.0/iArray[i][k];
			    }
			}
		}
		    
		    for (i = k+1; i < n; i++)
		    {
		    	if (iArray[i][k] != 0)
				{
					for (j = k; j < n; j++)
					{
		    	        iArray[i][j] = iArray[i][j] - iArray[k][j];
		        	}
		        }
			}
    }
    
	for (i = 0; i < n; i++)
	{
		for (j = 0; j < n; j++)
		{
			printf("%d\t", iArray[i][j]);
		}
		printf("\n");
	}
	
	return 0;
}