A segédprogramot Steib Roland írta.

A matematikai háttér:

Ez a program egy négyzetes mátrix determinánsát számolja ki.

A program használata és a szubrutinok meghívása

A program egy fõprogramból és két subroutinból áll.

Bemeneti adatok:

• A mátrix mérete (<=50)
  
• A mátrix elemei (soronként)
  

Kimeneti adatok:

• A mátrix determinánsa

A subroutine DTRM számolja ki a mátrix determinánsát a subroutine ELGS (Gauss elimináció) segítségével. A program az eredményt a képernyõre írja ki. Minden numerikus módszer hibával számol, ezért nem biztos, hogy a pontos megoldást kapjuk.

A két szubrutin:

- A mátrix determinansanak a kiszámitására.

Call DTRM(A,N,D)

Bemenõ paraméterek: N <=50 egész, A(50,50) valós

Kimenõ paraméterek: D valós

A DTRM pedig ezt a alszubrutint hivja meg az alábbi módon, ami a Gauss eliminációt számolja:

Call ELGS(A,N,IR,INDX)

Bemenõ paraméterek: N<=50 egész, A(50,50) valós, IR logikai

Kimenõ paraméterek: INDX(50) egész, A(50,50) valós

Példa a program használatára:

program deter1dr
REAL A(50,50),D
integer n
write (*,*)"Kerem a matrix sorainak a szamat:"
read(*,*)n

...

call DTRM (A,N,D)
write(*,*)Dstop
end

Próbafuttatás:

N=3

A(i,j)=1; 2; 3; 4; 5; 6; 7; 8; 10 (a mátrix elemei)

A kapott determináns:

D=-2,999997

Ha jól írtad be, ezt kell kapnod!

A segédprogramot itt leled: deter.for

[Vissza a fõoldalra]