-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
"John Randall" <[EMAIL PROTECTED]> writes:
> All the eigenvectors are for the matrix operating on the left. We have
> found v such that Av=lv, where A is the matrix and l is the eigenvalue.
> You need to premultiply v by A and then divide.
Showing my J skill and using your power script, something like ...
cureigenvec =: (%+/) @: power
neweigenvec =: mp cureigenvec
ev=: [: {. neweigenvec % cureigenvec
NB. Calculate the consistency ratio
n =: #
cr =: (ev - n) % <: @: n
NB. Results for all 4 of Coyle's 3x3 matrices
cr exp
0.014531883 NB. Coyle: 0.072
cr und
0.014531883 NB. Coyle: 0.026
cr rod
_4.4408921e_16 NB. Coyle: 0
cr pod
0.00061677663 NB. Coyle: 0
Bill
- --
Bill Harris http://facilitatedsystems.com/weblog/
Facilitated Systems Everett, WA 98208 USA
http://facilitatedsystems.com/ phone: +1 425 337-5541
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.1 (MingW32)
Comment: For more information, see http://www.gnupg.org
iD8DBQFEB1453J3HaQTDvd8RAsHBAJ4wrBd3XzfVzN7lUWwSGDNnpBAhmgCfcpQ1
tjvN8RPufSAmO1wGb1JMEI4=
=HUYF
-----END PGP SIGNATURE-----
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
