anchor.diff:  c
   c...  subroutine arguments
   c
  -      integer ndat,ma,npc
  +      integer ndat,ma,npc,ierr
         double precision x(ndat),y(ndat),sig(ndat),a(ma),covar(npc,npc)
         double precision xn(3,5)
  +      character errstrng*(*)
   c
   c...  local variables
   c
  
positive.code:       end
  c
  c version
  c $Id: lfit.f,v 1.2 2004/07/07 15:55:35 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code:           beta(j)=beta(j)+ym*wt
          end do
        end do
        call choldc2(covar,npc,ma,a,ierr,errstrng)
        call cholsl(covar,npc,ma,a,beta,beta)
        do j=1,ma
          a(j)=beta(j)
  