矩阵反码(c、c++、fortran或tksolver)

Question

我需要一个代码来使用行操作来确定nxn矩阵A的逆矩阵。我在写这段代码时感到很可笑。任何事都有帮助，我是一名研究生，已经多年没有上过编程课了。谢谢。

Answer 1

你可能并不是真的想要相反的结果。如果你试图解决一个方程系统，你最好使用LU分解。

你没有说你想用什么语言来写这个应用程序。Java有Apache Commons Math；Python有NumPy；FORTRAN有LinPack。选择一种语言并使用库；不要自己编写它。

Answer 2

如果你几年没上过编程课，我建议你用Scilab或MATLAB来学习矩阵运算。

如果是一次性计算，请查看WolframAlpha。

如果使用库是C++的一种选择，请查看Armadillo。

Answer 3

  SUBROUTINE MATINV(A,N)
  DIMENSION A(N,N)
  DO 1 I=1,N
     Z=A(I,I)
     A(I,I)=1.0
     DO 2 J=1,N
2       A(I,J)=A(1,J)/Z
     DO 1 K=1,N
        IF (K-I) 3,1,3
3          Z=A(K,I)
        A(K,I)=0.0
        DO 4 J=1,N
4          A(K,J)=A(K,J)-Z*A(I,J)
1 CONTINUE
  RETURN
  END