矩阵求逆的C实现

#include 
#include 
#include
int lud（double* m int n int* index） {
    int i=0 j=0 k=0 i_max = 0;
    double d_max=0 d_temp=0 d_dum=0;
    double pivot[10];
    for（i=0;i<10;i++）{
        pivot[i]=0;
    }

    for（i=0; i        for（d_max = 0.0 j=0; j            d_temp = fabs（m[i*n+j]）;
            if（ d_temp > d_max ）
		d_max = d_temp;
        }
        if（ d_max == 0.0 ）{
            return 0;
        }
        pivot[i] = 1.0 / d_max;
    }
   
    for（j=0; j        for（i=0; i            for（d_temp = m[i*n+j] k=0; k		d_temp -= m[i*n+k] * m[k*n+j];
            }
            m[i*n+j] = d_temp;
        }
        for（d_max=0.0 i=j; i            for（d_temp = m[i*n+j]