anchor.diff:  c     solution.
   c     Note that only the upper triangle is used (or available), as
   c     dmat is assumed to always be symmetric.
  +c     Note also that
   c
         include "implicit.inc"
   c
  
positive.code: c       
  
  c version
  c $Id: mat_1.f,v 1.3 2004/06/25 21:40:09 willic3 Exp $
  
  c Generated automatically by Fortran77Mill on Tue May 18 14:18:50 2004
  
  
negative.code:       double precision state(nstr,nstate),dstate(nstr,nstate),ee(nstr)
        double precision dmat(nddmat),prop(nprop),tmax
  c
        call dcopy(nstr,ee,ione,dstate(1,2),ione)
        call dcopy(nstr,state(1,1),ione,dstate(1,1),ione)
        call dspmv("u",nstr,one,dmat,dstate(1,2),ione,one,dstate(1,1),
       & ione)
  c
        return
        end
  