Exercise 3.1.1 B=Array[Min,{2,2}]; B//MatrixForm x*IdentityMatrix[2]-B//MatrixForm p=Det[%] lambda=x/.Solve[p==0,x] NullSpace[lambda[[1]]IdentityMatrix[2]-B] NullSpace[lambda[[2]]IdentityMatrix[2]-B] Exercise 3.1.2 A=Array[1&,{9,9}]; A//MatrixForm lambda=x/.Solve[Det[x*IdentityMatrix[9]-A]==0,x] NullSpace[lambda[[1]]IdentityMatrix[9]-A] NullSpace[lambda[[2]]IdentityMatrix[9]-A] Exercise 3.1.3 Q=Table[i/j,{i,3},{j,3}]; Q//MatrixForm lambda=x/.Solve[Det[x*IdentityMatrix[3]-Q]==0,x] All eigenvalues are obviously real. NullSpace[lambda[[1]]IdentityMatrix[3]-Q] NullSpace[lambda[[2]]IdentityMatrix[3]-Q] <