SUBROUTINE PMATRIX ( RHO, RSTAR1, RSTAR2, P, N, NORDER )
Argument Definitions (+ indicates altered content)
DOUBLE PRECISION  +P(*),       RHO(*),     RSTAR1(*),  RSTAR2(*)
INTEGER            N,          NORDER
Description
Compute an (N x N) matrix P = I - RHO(1)*RSTAR1, if norder = 1, 
        or   P = I - RHO(1)*RSTAR1 - RHO(2)*RSTAR2, if norder = 2, 
where I is an identity matrix of order N, 
      RHO is the vector containing RHO(1) and RHO(2), 
      RSTAR1 and RSTAR2 are (N x N) general matrices.
Source file:pmatrix.f
Local Variables (+ indicates altered content)
INTEGER           +I,         +ID,        +IJ,        +J