proc iml;

use lastAccumImpulse;

read all where(ric='B' )into x;

x1=x[1,1];

x2=x[1,2];

x3=x[1,3];

x4=x[1,4];

use CovInnovation;

read all where(ric='B' )into y ;

y1=y[1:4,1];

y2=y[1:4,2];

y3=y[1:4,3];

proc
iml;

use lastAccumImpulse;

read all into x;

x1=x[1,1];

x2=x[1,2];

x3=x[1,3];

x4=x[1,4];

use CovInnovation;

read all into y ;

y1=y[1:4,1];

y2=y[1:4,2];

y3=y[1:4,3];

y4=y[1:4,4];

z=T(x);

PV=x*y*z;

print x y z PV;