PROC IML;                                           
                                             

V0={34         6.77504    0         -2.19881,
     6.77504  21.60000    0         -1.75256,
     0         0         13        -10.87709,
    -2.19881  -1.75256  -10.87709   56.88000};

G0={15.30000   3.85503    0          1.97825,
     3.85503  10.80000    0          1.66260,
     0         0          5.20      -3.84567,
     1.97825   1.66260   -3.84567   11.37600};
                                         

W= {-2.50,
     7.20,
     0.00,
    10.80};
                  

INV_V0= INV(V0);                         
b = INV(V0)*G0*w;   PRINT b;
                         

G0w=G0*w;                           PRINT G0w;                            
V_I = b`*V0*b;                      PRINT V_I;  
V_H = w`*G0*w;                      PRINT V_H;  
R_IH  = SQRT (V_I/V_H);             PRINT R_IH; 

DLT_A = G0*b*(1 / SQRT(V_I) );      PRINT DLT_A;                                                
                                          

DUMY=J(2,2,0);k2={0,0};  PRINT DUMY k2;   
                                          

y=G0w//k2;          PRINT y;    
G1=G0(| ,1|);       PRINT G1;   
G2=G0(| ,2|);       PRINT G2;   
G02=G1||G2;         PRINT G02;  

C1=V0||G02;                     
C2=G02`||DUMY;                  
C =C1//C2;          PRINT C;    

INV_C=INV(C);
b0=INV_C*y;         PRINT INV_C;PRINT b0;  

kc=G0*b;                           
yc=G0w//kc;            PRINT kc yc;

O22=J(4,4,0);          PRINT O22;  
CC1=V0||G0;
CC2=G0||O22;                       
CC =CC1//CC2;          PRINT CC;   
                                   

INV_CC=INV(CC);
b02=INV_CC*yc;         PRINT INV_CC;PRINT b02;                            

ratio={0.25, 0.5,  1,  1};    
kp=ratio#kc;                  
PRINT ratio kc kp;
yp=G0w//kp;            PRINT G0w; PRINT kp yp;

bd1=INV_CC*yp;                                       
bd2=INV(G0)*kp;        PRINT INV_CC;PRINT bd1 bd2;   
RUN;