Abstract
Approximate
p
-point Leibniz derivation formulas as well as interpolatory Simpson quadrature sums adapted to oscillatory functions are discussed. Both theoretical considerations and numerical evidence concerning the dependence of the discretization errors on the frequency parameter of the oscillatory functions show that the accuracy gain of the present formulas over those based on the exponential fitting approach [L. Ixaru, "Computer Physics Communications", 105 (1997) 1--19] is overwhelming.