with(linalg):PT := matrix(5,5,[0,1/4,0,0,0,1,0,1/2,0,0,0,3/4,0,3/4,0,0,0,1/2,0,1,0,0,0,1/4,0]);Id := Matrix(5,5,shape=identity);det(PT-lambda*Id);