固有値問題(2)


つまらないことでひっかかっていましたが、ようやく動くようになりました。
すごい!
あとで実例はお見せしますが、下記のプログラムで多次元の固有値問題を解いてしまいます。
ひっかかっていたのはプログラム中にタブがあるとエラーになってしまうことに気づかず入力していました。
とりあえず解決。

Naigen =: 3 : 0
:
Y  =. y
ny =. # Y
X  =. x
if. X=2 do. Naigen2 y return. end.
dety =. det Y
my   =. minorir Y
Nmy  =: (X - 1) Naigen each my
Nmyp =. +/ > Nmy
e    =. Nmyp % ( >: i.ny )
x: Nans =.(((_1)^ny)*dety), e NB.extended precision
)

Naigen2=: 3 : 0
dety =.det y
my   =. minorir y
mdet =. detc my
e2m  =. (-mdet), 1
x:  Nans2=. (dety), e2m
)

each =: &.>
det  =: -/ . *
detc =: 3 : 0
 +/> (det each y)
)

omitr =: 3 : 0
:
(<< 0 do.
ci =.<(i,i) minor y
c  =. c, ci
i  =. i - 1
end.
c
)