 and G'
  approximately equal norms.  (See METHOD below for details.) For
  norm scaling, tau = MAX(1, ||A||, ||G||, ||Q||) where ||.||
  denotes the 1-norm (column sum norm).

</PRE>
<A name="Specification"><B><FONT SIZE="+1">Specification</FONT></B></A>
<PRE>
      SUBROUTINE MB04DY( JOBSCL, N, A, LDA, QG, LDQG, D, DWORK, INFO )
C     .. Scalar Arguments ..
      INTEGER           INFO, LDA, LDQG, N
      CHARACTER         JOBSCL
C     .. Array Arguments ..
      DOUBLE PRECISION  A(LDA,*), D(*), DWORK(*)