TER        ( MBMAX = MMAX + PMAX )
      PARAMETER        ( LDWORK = 2*NMAX*NMAX + MBMAX*(NMAX+PMAX) +
     $                      MAX( NMAX*(MAX( NMAX, MMAX, PMAX) + 5),
     $                     2*NMAX*PMAX + MAX( PMAX*(MBMAX+2),
     $                                        10*NMAX*(NMAX+1) ) ) )
*     .. Local Scalars ..
      DOUBLE PRECISION ALPHA, BETA, TOL1, TOL2
      INTEGER          I, INFO, IWARN, J, M, N, NR, NS, P
      CHARACTER*1      DICO, EQUIL, JOB, ORDSEL
*     .. Local Arrays ..
      