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:         beta(j)=zero
        end do
        do i=1,ndat
          call funcs(x(i),xn,afunc,ma,ndat)
          ym=y(i)
          sig2i=one/(sig(i)*sig(i))
          do j=1,ma
  