jzchu.c

适用于初学者学习

 #include "stdio.h"
 #define MAX 10

 int main(void)
 {
  float s[MAX][MAX],a[MAX][MAX],b[MAX][MAX];
  int i,j,k; float t,n;
  printf("A=B*C,知道A和B,求C\n");
  printf("请输入方阵的阶,不大于%d: ",MAX);
  scanf("%f",&n);
  printf("请按行序输入矩阵A:\n");
  for(i=0;i<n;i++)
   for(j=0;j<n;j++)
    scanf("%f",&s[i][j]);
  printf("请按行序输入矩阵B:\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("错误!B矩阵没有逆阵!\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];
   }
  }
  for(i=0;i<n;i++)
   for(j=0;j<n;j++)
   {
    a[i][j]=0;
    for(k=0;k<n;k++) a[i][j]+=s[k][j]*b[i][k];
   }
  printf("结果是:\n");
  for(i=0;i<n;i++)
  {
   for(j=0;j<n;j++) printf("%f\t",a[i][j]);
   printf("\n");
  }
  return 1;
 }