Gauss-Jordan法实矩阵求逆
下面是实现Gauss-Jordan法实矩阵求逆。
#include <stdlib.h>
#include <math.h>
#include <stdio.h>
int brinv(double a[], int n)
{ int *is,*js,i,j,k,l,u,v;
double d,p;
is=malloc(n*sizeof(int));
js=malloc(n*sizeof(int));
for (k=0; k<=n-1; k++)
{ d=0.0;
for (i=k; i<=n-1; i++)
for (j=k; j<=n-1; j++)
{ l=i*n+j; p=fabs(a[l]);
if (p>d) { d=p; is[k]=i; js[k]=j;}
}
if (d+1.0==1.0)
{ free(is); free(js); printf("err**not inv\n");
return(0);
}
if (is[k]!=k)
for (j=0; j<=n-1; j++)
{ u=k*n+j; v=is[k]*n+j;
p=a[u]; a[u]=a[v]; a[v]=p;
}
if (js[k]!=k)
for (i=0; i<=n-1; i++)
{ u=i*n+k; v=i*n+js[k];
p=a[u]; a[u]=a[v]; a[v]=p;
}
l=k*n+k;
a[l]=1.0/a[l];
for (j=0; j<=n-1; j++)
if (j!=k)
{ u=k*n+j; a[u]=a[u]*a[l];}
for (i=0; i<=n-1; i++)
你可能喜欢
- 线性方程组的高斯消去法
- 儿童教育方法
- 求逆矩阵
- 高中数学高考知识点总结
- 函授 数值 课件 chap4：解线性方程组的直接法1(高斯消去法)19页
- 用C语言实现解线性方程组的高斯消去法2页
- 线性方程组的高斯消去法34页
- 计算方法-第5章-2、解线性方程组的直接方法(5.5高斯消去法的变形)16页
- 解线性方程组的列主元素高斯消去法和LU分解法实验报告14页
- 高斯消去法和列主元高斯消去法解线性方程组的程序(C语言)5页
- 四岁儿童教育的方法5页
- 儿童教育方法集锦22页
- 小学儿童品德教育的主要方法6页
- 浅谈班主任在留守儿童管理中的教育方法2页
- 情绪行为障碍儿童的教育方法2页
- 少年儿童美术教育方法散论12页
- 求逆矩阵方法2页
- 求逆矩阵的方法3页
- Fortran求逆矩阵子程序3页
- 求逆矩阵3页
- 求逆矩阵3页
- 2014考研数学备考重点解析――求逆矩阵的方法4页
