我想将下列产生随机二项正态分布的协方差矩阵的最小特征值的程序重复1000遍得到1000个最小特征值。应该如何实现？
proc iml;
mu={0,0};
sig={1 0.9, 0.9 1};
a=half(sig)`;
seed=2343543;
p=2;
n=30;
x=rannor(j(p,n,seed));
y=a*x+mu*j(1,n,1);
ybar=y[,:];
yc=y-ybar*j(1,n,1);
s=yc*yc`/(n-1);
val=eigval(s);
minval=min(val);
print minval;
run;