anchor.diff:  c
    200                           return
         end
  -      subroutine reduc(nm,n,a,b,dl,ierr)
  +      subroutine util_reduc(nm,n,a,b,dl,ierr)
   c
         integer i,j,k,n,i1,j1,nm,nn,ierr
         double precision a(nm,n),b(nm,n),dl(n)
  
positive.code:       util_epslon = eps*dabs(x)
        return
        end
        subroutine tqlrat(n,d,e2,ierr)
  c
        integer i,j,l,m,n,ii,l1,mml,ierr
        double precision d(n),e2(n)
  
negative.code:  1000  ierr = l
   1001   return
        end
        subroutine util_tred1(nm,n,a,d,e,e2)
  c
        integer i,j,k,l,n,ii,nm,jp1
        double precision a(nm,n),d(n),e(n),e2(n)
  