anchor.diff:            beta(j)=beta(j)+ym*wt
           end do
         end do
  -      call choldc2(covar,npc,ma,a)
  +      call choldc2(covar,npc,ma,a,ierr,errstrng)
         call cholsl(covar,npc,ma,a,beta,beta)
         do j=1,ma
           a(j)=beta(j)
  
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: c
  c...  subroutine arguments
  c
        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
  