anchor.diff:          if(mhist(i).ne.izero) then
             if(mhist(i).gt.nhist.or.mhist(i).lt.0) then
               ierr=100
  +            errstrng="mathist"
               return
             end if
             ptmp(i)=histry(mhist(i),nstep+1)*prop(i)
  
positive.code:       integer nprop,imat,nstep,nhist,lastep,ierr
        integer mhist(nprop)
        logical matchg
        double precision ptmp(nprop),prop(nprop),histry(nhist,lastep+1)
  c
  c...  defined constants
  
negative.code:       end
  c
  c version
  c $Id: mathist.f,v 1.4 2004/07/05 19:40:17 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  