C++求矩阵的逆程序代码

1、求解逆矩阵源代码：#include<stdio.h>#include<stdlib.h>#include<math.h>#defineTINY1.0e-20voidinverse(double*,int);voidludcmp(double*,int,int*,double*);voidlubksb(double*,int,int*,double*);double*matrix(int,int,int,int);double*vector(int,int);voidfree_matrix(double*,int,int,int,int);voidfree_v

2、ector(double*,int,int);voidinverse(double*mat,intdim)inti,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<dim;i+)coli=0.0;colj=1.0;lubksb(mat,dim,indx,col);for(i=0;i<dim;

3、i+)yij=coli;for(i=0;i<dim;i+)for(j=0;j<dim;j+)matij=yij;free_matrix(y,0,dim-1,0,dim-1);free_vector(col,0,dim-1);free(indx);voidludcmp(double*a,intn,int*indx,double*d)inti,imax,j,k;doublebig,dum,sum,temp;double*vv;vv=(double*)malloc(unsigned)(n*sizeof(double);if(!vv)fprintf(stderr,"ErrorAl

4、locatingVectorMemoryn");exit(1);*d=1.0;for(i=0;i<n;i+)big=0.0;for(j=0;j<n;j+)if(temp=fabs(aij)>big)big=temp;if(big=0.0)fprintf(stderr,"SingularMatrixinRoutineLUDCMPn");for(j=0;j<n;j+)printf("%f",aij);printf("/n");exit(1);vvi=1.0/big;for(j=0;j<n;j+)fo

5、r(i=0;i<j;i+)sum=aij;for(k=0;k<i;k+)sum-=aik*akj;aij=sum;big=0.0;for(i=j;i<n;i+)sum=aij;for(k=0;k<j;k+)sum-=aik*akj;aij=sum;if(dum=vvi*fabs(sum)>=big)big=dum;imax=i;if(j!=imax)for(k=0;k<n;k+)dum=aimaxk;aimaxk=ajk;ajk=dum;*d=-(*d);vvimax=vvj;indxj=imax;if(ajj=0.0)ajj=TINY;if(j!=n-1)dum=1.0/ajj;for(i=j+1;i<n;i+)aij*=dum;free(vv);voidlubksb(double*a,intn,int*indx,double*b)inti,ip,j,ii=-1;doublesum;for(i=0;i<n;i+)ip=indxi;sum=bip;bip=bi;if(ii>=0)for(j=ii;j<i;j+)sum-=aij*bj;elsei