anchor.diff:  c
   c...  subroutine arguments
   c
  -      integer nprop,ierr
  +      integer ipauto,nprop,ierr
         integer iddmat(nstr,nstr)
         character errstrng*(*)
  -      double precision autofac,dmat(nddmat),prop(nprop)
  +      double precision tpois,tyoungs,dmat(nddmat),prop(nprop)
   c
   c...  local variables
   c
         double precision ptmp(10)
   c
         call dcopy(nprop,prop,ione,ptmp,ione)
  -      ptmp(2)=autofac*ptmp(2)
  +      if(ipauto.eq.ione) then
  +        ptmp(2)=tyoungs
  +        ptmp(3)=tpois
  +      end if
         call elas_mat_5(dmat,ptmp,iddmat,nprop,ierr,errstrng)
         return
         end
   c       
   
   c version
  -c $Id: mat_5.f,v 1.6 2005/01/18 19:40:36 willic3 Exp $
  +c $Id: mat_5.f,v 1.7 2005/02/24 00:14:45 willic3 Exp $
   
   c Generated automatically by Fortran77Mill on Tue May 18 14:18:50 2004
   
  
positive.code:       pr=prop(3)
        vis=prop(4)
        f1=third*e/(one-two*pr)
        f2=one/((one+pr)/e+half*alfap*deltp/vis)
        dmat(iddmat(1,1))=f1+two*third*f2
        dmat(iddmat(2,2))=dmat(iddmat(1,1))
        dmat(iddmat(3,3))=dmat(iddmat(1,1))
        dmat(iddmat(1,2))=f1-third*f2
        dmat(iddmat(1,3))=dmat(iddmat(1,2))
        dmat(iddmat(2,3))=dmat(iddmat(1,2))
        dmat(iddmat(4,4))=half*f2
        dmat(iddmat(5,5))=dmat(iddmat(4,4))
        dmat(iddmat(6,6))=dmat(iddmat(4,4))
  cdebug      write(6,*) "dmat:",(dmat(idb),idb=1,nddmat)
  
negative.code: c
  c...  subroutine arguments
  c
        integer nstate,nstate0,nprop,ierr
        integer iddmat(nstr,nstr)
        character errstrng*(*)
        logical matchg
        double precision state(nstr,nstate),dstate(nstr,nstate)
        double precision state0(nstate0),ee(nstr)
        double precision dmat(nddmat),prop(nprop),tmax
  c
  c...  included dimension and type statements
  