测绘程序设计c语言矩阵求逆调用,矩阵求逆用C语言实现

矩阵求逆的C/C++源程序

刚做的作业：

#include 

#include 

#include 

#include 

#define N 4

double A[N][N]={{1,3,2,13},{7,2,1,-2},{9,15,3,-2},{-2,-2,11,5}};

inline void swap(double &a,double &b){double c=a;a=b;b=c;};

int DinV(double A[N][N],int n)

{

int i,j,k;

double d;

int JS[N],IS[N];

for (k=0;k

{

d=0;

for (i=k;i

for (j=k;j

{

if (fabs(A[i][j])>d)

{

d=fabs(A[i][j]);

IS[k]=i;

JS[k]=j;

}

}

if (d+1.0==1.0) return 0;

if (IS[k]!=k)

for (j=0;j

swap(A[k][j],A[IS[k]][j]);

if (JS[k]!=k)

for (i=0;i

swap(A[i][k],A[i][JS[k]]);

A[k][k]=1/A[k][k];

for (j=0;j

if (j!=k) A[k][j]=A[k][j]*A[k][k];

for (i=0;i

if (i!=k)

for (j=0;j

if (j!=k) A[i][j]=A[i][j]-A[i][k]*A[k][j];

for (i=0;i

if (i!=k) A[i][k]=-A[i][k]*A[k][k];

}

for (k=n-1;k>=0;k--)

{

for (j=0;j

if (JS[k]!=k) swap(A[k][j],A[JS[k]][j]);

for (i=0;i

if (IS[k]!=k) swap(A[i][k],A[i][IS[k]]);

}

for (i=0;i

{

for (j=0;j

printf("  %1.4f",A[i][j]);

puts("");

}

}

void main()

{

DinV(A,4);

getch();

}