anchor.diff:        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
  +c #     subroutine lowesr(yy,iv,liv,lv,wv)
  +c      integer liv,lv
  +c      integer iv(*)
  +c      DOUBLE PRECISION yy(*),wv(*)
  +c
  +c      external ehg182,ehg192
  +c      if(.not.(iv(28).ne.172))then
  +c         call ehg182(172)
  +c      end if
  +c      if(.not.(iv(28).eq.173))then
  +c         call ehg182(173)
  +c      end if
  +c      call ehg192(yy,iv(2),iv(3),iv(19),iv(6),iv(14),wv(iv(13)),
  +c     +     wv(iv(34)),iv(iv(25)))
  +c      return
  +c      end
   
         subroutine lowesw(res,n,rw,pi)
   c Tranliterated from Devlin's ratfor
  
positive.code:       return
        end
  
        subroutine ehg140(iw,i,j)
        integer i,j
        integer iw(i)
        iw(i)=j
        return
        end
  
        subroutine ehg141(trl,n,deg,k,d,nsing,dk,delta1,delta2)
        double precision trl,delta1,delta2
  
negative.code:       return
        end
  
  c #     subroutine ehg192(y,d,n,nf,nv,nvmax,vval,lf,lq)
  c      integer d,i,i1,i2,j,n,nf,nv,nvmax
  c      integer lq(nvmax,nf)
  c      DOUBLE PRECISION i3
  c      DOUBLE PRECISION lf(0:d,nvmax,nf),vval(0:d,nvmax),y(n)
  
  c      do 3 i2=1,nv
  c         do 4 i1=0,d
  c            vval(i1,i2)=0
  c    4    continue
  c    3 continue
  c      do 5 i=1,nv
  c         do 6 j=1,nf
  c            i3=y(lq(i,j))
  c            do 7 i1=0,d
  c               vval(i1,i)=vval(i1,i)+i3*lf(i1,i,j)
  c    7       continue
  c    6    continue
  c    5 continue
  c      return
  c      end
  
        DOUBLE PRECISION function ehg128(z,d,ncmax,vc,a,xi,lo,hi,c,v,
       +     nvmax,vval)
  