NB. written by NAKANO Yoshihiro NB. JAPLA symposium 04/dec/2010 NB. Eigen Value and Eigen Vector LF0=: >@{:@p.@charn0 charn0 =: 3 : 0 In=. =@i.n=. # y X=.In i=.0 p=.1 for_k. >: i.n do. X =. y + / . * X trX =. +/(<0 1)|:X pk=.-k%~trX p=.p,pk X=.X+ pk * In i=.i + 1 end. |.p ) un=: 3 : '=@i.y' NB. unit matrix pfmt3=: 3 : '7j3 ": L:0 /: ~ |: y' sample=: 3 : 0 : ~./:"1~(i.x){"1 plist y ) plist=:i.@!A.i. nchar_lf =: 3 : 0 NB. by Nakano 2010 Sep.19 1; (LF0 y);charn0 y ) nchar_evec=: 3 : 0 NB. by Nakano 2010 Sep.19 NB. char_evec=: 3 : 0 NB. EIGEN=: {@> ; 1{ char_lf y EIGEN=: {@> ; 1{ nchar_lf y EIGEN2=: {@> L:0 EIGEN ^/L:0 |.i.# EIGEN ADJMAT=: char_lf_evec_sub y ANS=.<'' for_LF. i. # y do. NB. for_LF. 0 1 2 4 3 5 do. wr ' LF(Eig2) = ', ": LF{EIGEN2 wr ' ADJMART = ' wr ADJMAT wr 'lebel of ADJMAT = ' wr L. ADJMAT TMP=. +/> (> LF{ EIGEN2) * L:0 ADJMAT wr ' TMP = ' wr TMP ANS=. ANS,