NB.---Leverrier-Faddeev Algorithm--- NB. input: nxn array A NB. output: char. polynomial of A char=: 3 : 0 [EMAIL PROTECTED] for_k. >:i.n do. X=.y +/ . * X p=.p,pk=.-k%~+/(<0 1)|:X X=.X+pk*I end. |.p )
eval2=:[EMAIL PROTECTED]
;{:eval2 A
3 2 1
----
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
