options ls=80 pagesize = 60 nodate nonumber;
 
proc iml;
 
%include disc;
%include pvar;
 
s = {1.00  .376  .174  .048  .268  .352  .216  .179 -.036  .073 -.183,
     .376  1.00  .217 -.010  .119  .138  .097  .243  .070  .279  .026,
     .174  .217  1.00  .015  .078  .582  .032  .599 -.039  .682  .287,
     .048 -.010  .015  1.00 -.018 -.033 -.069 -.066  .063 -.114  .089,
     .268  .119  .078 -.018  1.00  .187  .095  .110 -.009  .132 -.023,
     .352  .138  .582 -.033  .187  1.00  .094  .531  .047  .452  .194,
     .216  .097  .032 -.069  .095  .094  1.00  .150  .035 -.006 -.266,
     .179  .243  .599 -.066  .110  .531  .150  1.00  .063  .564  .122,
    -.036  .070 -.039  .063 -.009  .047  .035  .063  1.00  .028  .084,
     .073  .279  .682 -.114  .132  .452 -.006  .564  .028  1.00  .354,
    -.183  .026  .287  .089 -.023  .194 -.266  .122  .084  .354  1.00};
 
ip = 1;   *ip = 1: input is correlation matrix;
	  *ip = 2: input is covariance matrix;
ix = 0;   *ip = 0:  Correlation matrix analysis requested;
          *ip = 1:  Covariance matrix analysis requested;
nbest = 5;
call pvar(s, ip,ix, nbest);