The vectorized operation will be,

proc iml;

n=5;

a1=i(n);
b1=loc(a1);
b2=b1+1;
*b3=b1-1;

a1(|b1(|2:n-1|)|)=2;
a1(|b2(|1:n-1|)|)=-1;
*a1(|b3(|2:n  |)|)=-1;



print a1 ;
quit;

谢谢，满足分析的要求。

谢谢了，还有一个问题：现在有个列向量比如a=(a1,a2,a3,a4,a5,a6)，现抽取该向量中的元素形成一个协方差矩阵如下(对称矩阵)：
A={a1 a2 a3 a4 a5 a6;
     a2  a1 a2 a3 a4 a5;
     a3  a2 a1 a2 a3 a4;
     a4  a3 a2 a1 a2 a3;
     a5  a4 a3 a2 a1 a2;
     a6  a5 a4 a3 a2 a1}
   谢谢了。

There is no problem to form a symmetric matrix as below. BUT a symmetric matrix is not necessary a covariance matrix.

proc iml;
a={10 2 3 4 5 6};



n=ncol(a);

symmetric_matrix=i(n);
b=loc(symmetric_matrix);
b2=loc(symmetric_matrix);

do i=0 to n-1;

   k=b+i;
   k2=b2-i;
   symmetric_matrix(|k(|1:n-i|)|)=a(|i+1|);
   if i>=1 then symmetric_matrix(|k2(|1+i:n |)|)=a(|i+1|);

end;

print symmetric_matrix;
quit;



quit;

谢谢你的程序，你说的很对，对称矩阵未必是协方差矩阵，不过我的结果的确是协方差矩阵，叩谢老兄了。