C++求矩阵的逆程序代码(4)

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

b[ip] = b[i];

if (ii>=0)

for (j=ii;j<i;j++) sum -= a[i][j] * b[j]; else if (sum) ii = i;

b[i] = sum;

}

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

{

sum = b[i];

for (j=i+1;j<n;j++) sum -= a[i][j] * b[j]; b[i] = sum / a[i][i];

}

}