供参考：
proc iml;
start kt(t);
kt=k(t) ; /*function k expression*/
return kt;
finish;

start main;
t=[] ; /**注意矩阵t的维数，以及由函数k(t)返回到kt的维数，比如10*10 **/
H=kt(t);
print H;
finish;
run;
quit;

多谢了