求解逆矩阵源代码:＃include 〈stdio.h>＃include 〈stdlib.h〉#include #define TINY 1。0e—20void inverse（double＊*，int);void ludcmp（double**， int, int＊， double*)；void lubksb(double＊*, int， int*, double＊）；double *＊matrix(int，int,int,int);double *vector(int,int);void free_matrix（double＊*,int，int，int，int);void free_vector（double＊,int，int）；void inverse（double **mat, int dim){ int i，j，*indx； double ＊＊y,d,＊col; y = matrix(0,dim—1，0,dim—1)； indx = （int *）malloc((unsigned）(dim＊sizeof(int))）; col = vector（0,dim—1)； ludcmp（mat，dim，indx,＆d)； for (j=0;j〈dim；j++） ｛ for （i=0;i= big） ｛ big = dum； imax = i； } ｝ if (j ！= imax） ｛ for (k=0；k〈n;k++） { dum = a［imax][k］； a[imax][k] = a[j］［k]； a[j][k］ = dum; } ＊d = -(*d）； vv［imax] = vv[j］； } indx［j] = imax; if （a[j][j］ == 0。0） a[j]［j] = TINY; if （j != n-1) { dum = 1。0 / a［j][j］; for （i=j+1；i〈n;i++） a[i］[j] ＊= dum; ｝ ｝ free（vv）;}void lubksb（double ＊＊a, int n, int *indx， double ＊b）｛ int i，ip,j，ii=-1； double sum； for （i=0;i=0；i—-) { sum = b［i]； for （j=i+1;j
