proc iml;
start main;


*P={19881 112.8 253.8,
   112.8  64   108,
   253.8 108   48.6};

P={19881 112.8 ,
   112.8  64   };

*G={6958 -23.3 -58.4,
   -23.3 9.6  14.1,
   -58.4 14.1 48.6};

G={6958 -23.3 -58.4,
   -23.3 9.6  14.1};

a={.022,
   .5,
   .5};

*w=inv(P)*G*a; 
w={1,
   1};      
cov=w`*G*a;
VZ=w`*p*w;
SZ=sqrt(vz);
DH=cov*inv(SZ);
DG=w`*G*inv(SZ);
print w DH DG;


finish main;
run;
quit;