anchor.diff:  340     dd = rho * rho
           sd = 0.0_wp
           ss = 0.0_wp
  -        do 350 i = 1, n
  +        do i = 1, n
               if (abs(dx(i)) >= 1.0e-6_wp*rho) dd = dd - dx (i) ** 2
               sd = sd + dx (i) * sdirn (i)
  -350     ss = ss + sdirn (i) ** 2
  +            ss = ss + sdirn (i) ** 2
  +        end do
           if (dd <= 0.0_wp) go to 490
           temp = sqrt (ss*dd)
           if (abs(sd) >= 1.0e-6_wp*temp) temp = sqrt (ss*dd+sd*sd)
  
positive.code:                             denom = min (denom, temp)
                          end if
                      end if
  570             continue
                  if (denom == 0.0_wp) then
                      parmu = 0.0_wp
                  else if (cmax-cmin < parmu*denom) then
  
negative.code:                 xopt (2) = 1.0_wp
                  xopt (3) = 2.0_wp
                  xopt (4) = - 1.0_wp
                  
              else if (nprob == 9) then
      !
      !     This problem is taken from page 111 of Hock and Schittkowski's
  