proc iml;

start main;


P={19881 112.8 ,

   112.8  64   };


G={6958 -23.3 -58.4,

   -23.3 9.6  14.1};


a={.022,

   .5,

   .5};

w=inv(P)*G*a;

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;