行列の逆符号（c、C++、fortranまたはtksolver）

Question

行操作を使ってnxn行列Aの逆行列を求めるコードが必要です。 私はこのコードを書いてばかげている。何かが助けてくれる、私は大学院生で、何年もプログラミングクラスを取らなかった。ありがとう。

Answer 1

何年もプログラミングクラスを習得していない場合は、ScilabまたはMATLABを使用することをお勧めします。

1回限りの計算であれば、WolframAlphaをご覧ください。

ライブラリを使用することがC++のオプションである場合は、Armadilloを参照してください。

Answer 2

あなたは本当に逆を望んでいないかもしれません。方程式の系を解くことを試みるなら、あなたはLU分解を使う方が良いでしょう。

あなたはこのアプリケーションを書いてみたい言語を言っていません。JavaにはApache Commons Mathがあります。 PythonにはNumPyがあります。 FORTRANにはLinPackがあります。言語を選択してライブラリを使用する。自分で書いてはいけません。

Answer 3

Numerical Recipies Onlineをご覧ください。 Fortran 77、Ansi C、Fortran 90のコードは無料で、最善の方法についての手がかりを得ることができます。

Answer 4

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