options nodate nonumber ps=55 ls=80; title 'RPD -- Example 2.11'; proc iml; A = {10 3, 3 8}; lambda = eigval(A); Z = eigvec(A); lambda1i = lambda(|1,|)*i(2); lambda2i = lambda(|2,|)*i(2); detA_lambda1 = det(A-lambda1i); detA_lambda2 = det(A-lambda2i); print A; print lambda; print Z; print detA_lambda1 detA_lambda2; run;