anchor.diff:          do l=1,ngauss
             gauss(1,l)=r(l)*g1
             gauss(2,l)=s(l)*g1
  -          gauss(3,l)=-two*third
  +          gauss(3,l)=-(two*third)
             gauss(4,l)=w1
           end do
           gauss(3,5)=two/five
           gauss(4,5)=c6/c2
         end if
  +c
  +      infetype(1)=ngauss
  +      infetype(2)=nen
  +      infetype(3)=nec
  +      infetype(4)=nee
  +c
         do l=1,ngauss
           do i=1,nen
             rr=half*(one+r(i)*gauss(1,l))
  
positive.code:       end
  c
  c version
  c $Id: plinpyr.f,v 1.1 2004/04/14 21:18:30 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code: c
  c...  local variables
  c
        integer nen,ngauss,nec,nee,i,l,nshsize,ngssize
        double precision rr,ss,tt,g1,w1
  c
  c...  definitions
  