# # Written by: # -- # John L. Weatherwax 2009-04-21 # # email: wax@alum.mit.edu # # Please send comments and especially bug reports to the # above email address. # #----- # Exercise 7: # a = 0 m = matrix( data=c( 1.0, 0.9, a, 0.9, 1.0, 0.9, a, 0.9, 1.0 ), nrow=3, ncol=3, byrow=TRUE ) E = eigen(m) w = E$vectors[,3] # use these weights print( w ) print( t(w) %*% m %*% w )