Abstract
Exponentially fitted versions of the extended one-step methods for first-order ODEs are constructed and their properties are examined. It is shown how formulae for optimal frequencies can be constructed. Some numerical experiments illustrate the use of the developed algorithms.