#include <stdio.h>
void main()
{
    float a[3][3],b[3][3]={{1,0,0},{0,1,0},{0,0,1}};
    int i,j,k;
    scanf("%f%f%f%f%f%f%f%f%f",&a[0][0],&a[0][1],&a[0][2],&a[1][0],&a[1][1],&a[1][2],&a[2][0],&a[2][1],&a[2][2]);
    for(k=0;k<=2;k++)
    {
        for(j=0;j<=2;j++)
        {
            a[k][j]=a[k][j]/a[k][k];
            b[k][j]=b[k][j]/a[k][k];
        }
        for(i=0;i<=2;i++)
        {
            if(i!=k)
            {
                for(j=0;j<=2;j++)
                {
                    a[i][j]=a[i][j]-a[k][j]*a[i][k];
                    b[i][j]=b[i][j]-b[k][j]*a[i][k];
                }
            }
        }
    }
    printf("%f,%f,%f\n%f,%f,%f\n%f,%f,%f\n",a[0][0],a[0][1],a[0][2],a[1][0],a[1][1],a[1][2],a[2][0],a[2][1],a[2][2]);
    printf("%f,%f,%f\n%f,%f,%f\n%f,%f,%f",b[0][0],b[0][1],b[0][2],b[1][0],b[1][1],b[1][2],b[2][0],b[2][1],b[2][2]);
}
这是我编的一段求逆矩阵的代码，算法应该没问题，因为a输出时为单位矩阵，但是b却原样输出！这是什么问题？？为什么b没变化？！

自己顶一下！

#include <stdio.h>
void main()
{
    float a[3][3],b[3][3]={{1,0,0},{0,1,0},{0,0,1}};
    int i,j,k;
    scanf("%f%f%f%f%f%f%f%f%f",&a[0][0],&a[0][1],&a[0][2],&a[1][0],&a[1][1],&a[1][2],&a[2][0],&a[2][1],&a[2][2]);
    for(k=0;k<=2;k++)
    {
        for(j=0;j<=2;j++)
        {
            a[k][j]=a[k][j]/a[k][k];
            b[k][j]=b[k][j]/a[k][k];    // 上一句使得a[k][k]=1，所以这一句没意义
        }
        for(i=0;i<=2;i++)
        {
            if(i!=k)
            {
                for(j=0;j<=2;j++)
                {
                    a[i][j]=a[i][j]-a[k][j]*a[i][k];
                    b[i][j]=b[i][j]-b[k][j]*a[i][k];
                    //因为 当i=j时b[i][j]=1,当i!=j时b[i][j]=0, 所以b[k][j]*a[i][k]=0
                    //所以这一句b的值也没变，
                }
            }
        }
    }
    printf("%f,%f,%f\n%f,%f,%f\n%f,%f,%f\n",a[0][0],a[0][1],a[0][2],a[1][0],a[1][1],a[1][2],a[2][0],a[2][1],a[2][2]);
    printf("%f,%f,%f\n%f,%f,%f\n%f,%f,%f",b[0][0],b[0][1],b[0][2],b[1][0],b[1][1],b[1][2],b[2][0],b[2][1],b[2][2]);
}
－－－－－－－－－－－－－－－－－－－－
由此我认为你的算法问题大了，自己再好好找找

谢谢，我现在搞明白了！