options ls=76 ps=55;
proc iml;
c={1 9 3 4, 1 0 3 12};
d={1 0 0 4, 1 0 0 1};
a=c+d;
b=c-d;
print c;
print d;
print a;
print b;
aa=c//d;
print aa;
aat=aa`;
print aat;
aat2=aat*2;
dt=d`;
print dt;
cdt=c*dt;
print cdt;
rowcdt=nrow(cdt);
print rowcdt;
colcdt=ncol(cdt);
print colcdt;
print aat2;
mdaat2= diag(aat2);
print mdaat2;
vdaat2= vecdiag(aat2);
print vdaat2;
ident=i(3);
print ident;
unity=j(3, 3, 1);
print unity;

run;