怎么写一个宏程序

shiyiming

proc iml;
x={30    0     105,
      45   100   35,
       0    75    70};
z=diag({150,250,350});
a=x*inv(z);
zt=diag({200,300,400});
ut=diag({170,180,120});
vt=diag({90,200,180});
y={1,1,1};

x1=a*zt;

**step 1**;
u1=diag(x1*y);
r1=ut*inv(u1);
x2=r1*x;
v1=diag(t(y)*x2);
s1=vt*inv(v1);
x3=x2*s1;

**step 2**;
u2=diag(x3*y);
r2=ut*inv(u2);
x4=r2*x3;
e1=r2*y;
v2=diag(t(y)*x4);
s2=vt*inv(v2);
x5=x4*s2;

**step 3**;
u3=diag(x5*y);
r3=ut*inv(u3);
x6=r3*x5;
e2=r3*y;
v3=diag(t(y)*x6);
s3=vt*inv(v3);
x7=x6*s3;

**step 4**;
u4=diag(x7*y);
r4=ut*inv(u4);
x8=r4*x7;
e3=r4*y;
v4=diag(t(y)*x8 );
s4=vt*inv(v4);
x9=x8*s4;

print  e1 e2 e3 ;


如上程序，以后还有很多步要做
怎么写一个宏程序以提高效率呢？
诚谢！

shiyiming

%macro calc(p);
%do n=2 %to %eval(&p+1);
u&n=diag(x%eval(2*&n-1)*y);
r&n=ut*inv(u&n);
x%eval(2*&n)=r&n*x%eval(2*&n-1);
e%eval(&n-1)=r&n*y;
v&n=diag(t(y)*x%eval(2*&n));
s&n=vt*inv(v&n);
x%eval(2*&n+1)=x%eval(2*&n)*s&n;
%end;
%mend;
proc iml;
x={30 0   105,
   45 100 35,
    0 75  70};
z=diag({150,250,350});
a=x*inv(z);
zt=diag({200,300,400});
ut=diag({170,180,120});
vt=diag({90,200,180});
y={1,1,1};
x1=a*zt;

**step 1**;
u1=diag(x1*y);
r1=ut*inv(u1);
x2=r1*x;
v1=diag(t(y)*x2);
s1=vt*inv(v1);
x3=x2*s1;

%calc(3);
print e1 e2 e3 ;
quit;

shiyiming

谢谢jerry老兄了