proc iml;
/*
id = i(5);
print id;

jd = j(5);
print jd;

a ={2 0 0, 0 4 0, 0 0 8};
b = inv(a);
print b;

a1 ={1 2 3, 0 4 0, 2 4 6};
b1 = inv(a1);
print b1;



A =  {1 -5, -5 1};
call eigen(lambda, e, A);
print lambda;
print e;

e1 = eigvec(A);
print e1;

l1= eigval(A);
print l1;

*/

x = {1 ,2 ,3};
A = diag(x);
print A;