function [] = prob_10_7 % PROB_10_7 - % % Written by: % -- % John L. Weatherwax 2005-08-28 % % email: wax@alum.mit.edu % % Please send comments and especially bug reports to the % above email address. % %----- % Pt (a): R = sym( rosser ) e = eig( R ) [ ignore, k ] = sort(double(e)); es = e(k) % sort(eig(R)) seems to sort alphebetically. % Pt (b): p = poly(R) f = factor(p) pretty(f) % Pt (c): e = eig( sym(rosser) ) r = eig( rosser ) sort(double(e))-sort(r) double(e(k)-r) % Pt (d): R(1,1) = 612; eig( R )