anchor.diff:    ! output solver statistics
     print *, '  '
     print *, 'Final Solver Statistics:'
  -  print '(3(A,i7),A)', '   Total internal solver steps =',iout(3),&
  -       ' (acc =',iout(5),',  conv =',iout(6),')'
  -  print '(2(A,i7))', '   Total RHS evals:  Fe =',iout(7),',  Fi =',iout(8)
  +  print '(3(A,i7),A)', '   Total internal solver steps =', iout(3), &
  +       ' (acc =', iout(5), ',  conv =', iout(6), ')'
  +  print '(2(A,i7))', '   Total RHS evals:  Fe =', iout(7), &
  +       ',  Fi =', iout(8)
     print '(A,i7)', '   Total linear solver setups =', iout(9)
  -  print '(A,i7)', '   Total number of Jacobian evaluations =',iout(18)
  -  print '(A,i7)', '   Total number of Newton iterations =',iout(11)
  -  print '(A,i7)', '   Total number of linear solver convergence failures =',iout(12)
  -  print '(A,i7)', '   Total number of error test failures =',iout(10)
  -
  -  ! check final solution
  +  print '(A,i7)', '   Total Jacobian evaluations =', iout(18)
  +  print '(A,i7)', '   Total Newton iterations =', iout(11)
  +  print '(A,i7)', '   Total linear solver convergence failures =', &
  +       iout(12)
  +  print '(A,i7)', '   Total error test failures =', iout(10)
  +  print *, '  '
  +
  +  ! check final solution against reference values for problem
     call solut(neqn, Tf, ytrue)
     print *, '     y(Tf) =', y
     print *, '  yref(Tf) =', ytrue
  
positive.code:   ipar(2) = imex
    if (dense_order == -1)  denseout = .false.
  
    ! call problem setup routine
    call prob(fullnm, problm, type, neqn, ndisc, t, &
              numjac, nummas, mlmas, mumas, ind)
    print *,'full problem name: ',fullnm
  
negative.code: !-----------------------------------------------------------------
  
  
  subroutine farkdjac(neq,t,y,fy,DJac,h,ipar,rpar,wk1,wk2,wk3,ier)
  !-----------------------------------------------------------------
  ! Jacobian computation routine
  !-----------------------------------------------------------------
  