anchor.diff:            gauss(1,l)=g1
             gauss(2,l)=g1
             gauss(3,l)=g1
  -          if(l.ne.ngauss) gauss(l,l)=g2
  +          if(l.ne.ione) gauss(l-ione,l)=g2
             gauss(4,l)=fourth*tetvol
           end do
         end if
  
positive.code:       end
  c
  c version
  c $Id: pquadtet.f,v 1.3 2004/08/02 21:21:54 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code:         rr=gauss(1,l)
          ss=gauss(2,l)
          tt=gauss(3,l)
  c*        uu=g1
  c*        if(l.eq.ione) uu=g2
  c*        if(intord.eq.itwo) uu=fourth
          uu=one-rr-ss-tt
  cdebug        write(6,*) "l:",l
  cdebug        write(6,*) "rr,ss,tt,uu:",rr,ss,tt,uu
          sh(4,1,l)=uu*(two*uu-one)
          sh(1,1,l)=one-four*uu
          sh(2,1,l)=one-four*uu
          sh(3,1,l)=one-four*uu
          sh(4,2,l)=rr*(two*rr-one)
          sh(1,2,l)=four*rr-one
          sh(2,2,l)=zero
          sh(3,2,l)=zero
          sh(4,3,l)=ss*(two*ss-one)
          sh(1,3,l)=zero
          sh(2,3,l)=four*ss-one
          sh(3,3,l)=zero
          sh(4,4,l)=tt*(two*tt-one)
          sh(1,4,l)=zero
          sh(2,4,l)=zero
          sh(3,4,l)=four*tt-one
          sh(4,5,l)=four*rr*uu
          sh(1,5,l)=four*(uu-rr)
          sh(2,5,l)=-(four*rr)
          sh(3,5,l)=-(four*rr)
          sh(4,6,l)=four*rr*ss
          sh(1,6,l)=four*ss
          sh(2,6,l)=four*rr
          sh(3,6,l)=zero
          sh(4,7,l)=four*ss*uu
          sh(1,7,l)=-(four*ss)
          sh(2,7,l)=four*(uu-ss)
          sh(3,7,l)=-(four*ss)
          sh(4,8,l)=four*tt*uu
          sh(1,8,l)=-(four*tt)
          sh(2,8,l)=-(four*tt)
          sh(3,8,l)=four*(uu-tt)
          sh(4,9,l)=four*rr*tt
          sh(1,9,l)=four*tt
          sh(2,9,l)=zero
          sh(3,9,l)=four*rr
          sh(4,10,l)=four*ss*tt
          sh(1,10,l)=zero
          sh(2,10,l)=four*tt
          sh(3,10,l)=four*ss
  cdebug        dbsum=zero
  cdebug        do idb=1,nen
  cdebug          write(6,*) "sh:",(sh(jdb,idb,l),jdb=1,4)
  