anchor.diff:  
       end if
   
  -    f = PEN * (plim - f)
  +    f = PEN * max(0.0D0, plim - f) ** 2.0D0
   
       ! Quadratic term
   
  
positive.code: 
    open(99, FILE='ccp.out')
  
    write(99, FMT=020) prob, n, np, mi, f, feas, f - PEN * (plim - p), &
         p, fcnt, flag
  
    close(99)
  
negative.code: 
    ! np = INT(n * npfrac)
  
    me = 0
  
    ! mi = 10
  
  