options nodate nonumber ps=55 ls=80;

title 'RPD -- Example 2.1';

proc iml;

T={1 2, 4 5, 3 0};
W={-1, 3};

TW=T*W;

print T W TW;

run;