anchor.diff:  c
   cdebug      write(6,*) "Hello from elas_mat_5_f!"
   c
  +      call fill(dmat,zero,nddmat)
         e=prop(2)
         pr=prop(3)
         pr1=one-pr
  
positive.code: c
  cdebug      write(6,*) "Hello from td_matinit_5_f!"
  c
        tmax=big
        e=prop(2)
        pr=prop(3)
  
negative.code:       include "rgiter_def.inc"
        include "ntimdat_def.inc"
  c
        call td_matinit_5(state,dstate,dmat,prop,rtimdat,
       & rgiter,ntimdat,iddmat,tmax,nstate,nprop,matchg,ierr,errstrng)
        call td_strs_5(state,dstate,ee,dmat,prop,rtimdat,rgiter,ntimdat,
       & iddmat,tmax,nstate,nprop,matchg,ierr,errstrng)
  