Friends:
Does somebody know why for some matriz fail the verb
below?
Does somebody know another verb in J in order to
calculate eigenvalues?
The verb in question is:
NB. ¡Oh No!, ¡Not Eigenvalues Again!
NB. script of Eugene McDonnell
NB. Nota:
NB. Se tarda mucho con matrices mayores de 9
renglones
NB. y muchos de los resultados no concuerdan con
los
NB. del programa de las utiler¡as de APL.
NB. En este caso el resultado si coincide con el de
APL
A=: 33 16 72 , _24 _10 _57 ,: _8 _4 _17
NB. Sintaxis:
NB. eigen A
NB. 3 2 1
eigen=: 3 : 0
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=:>@{:@[EMAIL PROTECTED]@db
(<'eigenvalores con el m ximo'),:<(<<>./r),<"0 r=.cm
y.
)
¡S A L U D!
Alfonso Salazar
____________________________________________________________________________________
¡Sé un mejor ambientalista!
Encuentra consejos para cuidar el lugar donde vivimos.
http://mx.yahoo.com/promos/mejorambientalista.html
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
