lit
    JINRLIB

KIM - .

F240

: .., .. :

popaa KIM c coc a coc op co cpo ap oo o c oa ooa opa acaoo a (coc op oppoa a ).

:

: - SUBROUTINE
: - KIM

:

CALL KIM(H,N,IEGEN,U,NR,NDIM,X,NO,EPS), :

H - (REAL*8) acc, oopo acaa ca cpa apa, H(NDIM,NDIM);
N - (INTEGER) opo ap;
IEGEN - (INTEGER) c IEGEN=0, cc coc a coc op.
c IEGEN pao 0, cc oo coc a;
U - (REAL*8) apa, pacoo coc op, U(NDIM,NDIM);
NR - (INTEGER) co pa, po opopao oc aao ooc;
NDIM - (INTEGER) papoc acca H, NDIM ≥ N;
X - (REAL*8) pao acc, X(NDIM);
NO - (INTEGER) c NO=0 IEGEN=0, o p aao pao popa a co U aocc a apa.
c NO pao 0, IEGEN=0, o apa coc opo oac a ap, oopa a pacooa a c U p opa opopa;
EPS - (REAL*8) aaa ooc c coc opo.
ooc c coc a EPS2.

:

  1. opopaa po ap H aoao .
    p o aoa ap c
      c o c   a .
  2. I- coc op pacoaac ap c opao: U(K,I), K=1,N.

:

  1. B .
    . M, M, 1962, .269-278.
  2. .K., B.H. B .
    , M, 1963.
:
       IMPLICIT REAL*8 (A-H,O-Z)
       DIMENSION H ( 4,4 ), U ( 4,4 ), X ( 4 )
       EPS = 0.0000001D0
       H ( 1,1 ) = 1.00D0
       H ( 1,2 ) = 0.42D0
       H ( 1,3 ) = 0.54D0
       H ( 1,4 ) = 0.66D0
       H ( 2,1 ) = 0.42D0
       H ( 2,2 ) = 1.00D0
       H ( 2,3 ) = 0.32D0
       H ( 2,4 ) = 0.44D0
       H ( 3,1 ) = 0.54D0
       H ( 3,2 ) = 0.32D0
       H ( 3,3 ) = 1.00D0
       H ( 3,4 ) = 0.22D0
       H ( 4,1 ) = 0.66D0
       H ( 4,2 ) = 0.44D0
       H ( 4,3 ) = 0.22D0
       H ( 4,4 ) = 1.00D0
       CALL KIM ( H, 4, 0, U, NR, 4, X, 0, EPS )
       . . .
:
       EIGENVALUES AND EIGENVECTORS OF A REAL SYMMETRIC MATRIX
 
    H= 
       2.322748800        0.000000001        0.000000000        0.000000000
       0.000000001        0.638283803        0.000000000        0.000000000
       0.000000000        0.000000000        0.796706689        0.000000000
       0.000000000        0.000000000        0.000000000        0.242260708
 
    X= 
       0.579642502        0.459996664        0.433459111        0.514325614
      -0.380449881        0.850275474        0.035889606       -0.361941214
       0.050328449       -0.237226458        0.812846171       -0.529595843
      -0.718845953       -0.095698981        0.387435463        0.569206433
 
       2.322748800        0.638283803        0.796706689        0.242260708
       NUMBER ITERATIONS=   14 M=    5
      
       2.322748800        0.796706690        0.638283800        0.242260710
       EXACT EVALUE


home up e-mail