options nodate nonumber ps=55 ls=80; title 'RPD -- Example 2.10'; proc iml; A = {5 2 -1,2 2 2,-1 2 5}; A = A/6; AA = A*A; rankA = trace(A); print A; print AA; print rankA; run;