Sorry to be replying so late, but I have been away.
Alfonso Salazar <[EMAIL PROTECTED]> wrote:
> consider in Eugene McDonnell's algorithm:
> IR=:@([EMAIL PROTECTED]@])
> Does somebody know why for some matriz fail the verb
below?
I think it works, unless there is some incompatibility between j504 (which
Alfonso appears to be using) and j601.
Here is the example, modified for j601:
9!:7 '+++++++++|-'
A =: 33 16 72 , _24 _10 _57 ,: _8 _4 _17
eigen=: 3 : 0
IR=:@([EMAIL PROTECTED]@]) NB. Here is yet
B1=:<"1
B0=:<"0
d=:(<0 1)&|:
dif=:-/@,:
difb=:dif&.>
ppr=:+//.@(*/)
pprb=:ppr&.>
detp=:>@(difb/ .pprb)
db=:B1@(d ,. _1:) d@([EMAIL PROTECTED]@])} B0
cm=:>@{:@p. @ detp @ db
(<'eigenvalores con el m ximo'),:<(<<>./r),<"0 r=.cm y
)
eigen A
+--------------------------+
|eigenvalores con el m ximo|
+--------------------------+
|+---+-+-+-+ |
||+-+|3|2|1| |
|||3|| | | | |
||+-+| | | | |
|+---+-+-+-+ |
+--------------------------+
John
Delete & Prev | Delete & Next
Move to:
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
