anchor.diff:        return
         end
   
  -      subroutine ehg140(iw,i,j)
  -      integer i,j
  -      integer iw(i)
  -      iw(i)=j
  -      return
  -      end
  +c#      subroutine ehg140(iw,i,j)
  +c      integer i,j
  +c      integer iw(i)
  +c      iw(i)=j
  +c      return
  +c      end
   
         subroutine ehg141(trl,n,deg,k,d,nsing,dk,delta1,delta2)
         double precision trl,delta1,delta2
  
positive.code:       return
        end
  
        subroutine lowesr(yy,iv,liv,lv,wv)
        integer liv,lv
        integer iv(*)
        DOUBLE PRECISION yy(*),wv(*)
  
        external ehg182,ehg192
        if(.not.(iv(28).ne.172))then
           call ehg182(172)
        end if
        if(.not.(iv(28).eq.173))then
           call ehg182(173)
        end if
        call ehg192(yy,iv(2),iv(3),iv(19),iv(6),iv(14),wv(iv(13)),
       +     wv(iv(34)),iv(iv(25)))
        return
        end
  
        subroutine lowesw(res,n,rw,pi)
  c Tranliterated from Devlin's ratfor
  
negative.code:       return
        end
  
  c      integer function ehg138(i,z,a,xi,lo,hi,ncmax)
  c      logical i1
  c      integer i,j,ncmax
  c      integer a(ncmax),hi(ncmax),lo(ncmax)
  c      DOUBLE PRECISION xi(ncmax),z(8)
  cc     descend tree until leaf or ambiguous
  c      j=i
  cc     top of while loop
  c    3 if(a(j).ne.0)then
  c         i1=(z(a(j)).ne.xi(j))
  c      else
  c         i1=.false.
  c      end if
  c      if(.not.(i1))goto 4
  c         if(z(a(j)).le.xi(j))then
  c            j=lo(j)
  c         else
  c            j=hi(j)
  c         end if
  c         goto 3
  cc     bottom of while loop
  c    4 ehg138=j
  c      return
  c      end
  
        subroutine ehg106(il,ir,k,nk,p,pi,n)
  
  