nizhen.c

来自「适用于初学者学习」· C语言 代码 · 共 47 行

 #include "stdio.h"
 #define MAX 10

 int main(void)
 {
  float a[MAX][MAX],b[MAX][MAX];
  int i,j,k; float t,n;
  printf("请输入原方阵的阶,不大于%d: ",MAX);
  scanf("%f",&n);
  printf("请按行序输入矩阵:\n");
  for(i=0;i<n;i++)
   for(j=0;j<n;j++)
   {
    scanf("%f",&a[i][j]);
    if(i!=j) b[i][j]=0;
     else b[i][j]=1;
   }
  for(k=0;k<n;k++)
  {
   if(!a[k][k])
    {printf("错误!该矩阵没有逆阵!\n");return 0;}
   t=a[k][k];
   for(j=k;j<n;j++) a[k][j]/=t;
   for(j=0;j<n;j++) b[k][j]/=t;
   for(i=0;i<k;i++)
   {
    t=a[i][k];
    if(!t) continue;
    for(j=k;j<n;j++) a[i][j]-=t*a[k][j];
    for(j=0;j<n;j++) b[i][j]-=t*b[k][j];
   }
   for(i=k+1;i<n;i++)
   {
    t=a[i][k];
    if(!t) continue;
    for(j=k;j<n;j++) a[i][j]-=t*a[k][j];
    for(j=0;j<n;j++) b[i][j]-=t*b[k][j];
   }
  }
  printf("原方阵的逆阵是:\n");
  for(i=0;i<n;i++)
  {
   for(j=0;j<n;j++) printf("%f\t",b[i][j]);
   printf("\n");
  }
  return 1;
 }