Ordinary differential equations arise in many physical applications. The conventional means to calculate a numerical solution of such equations is by using Runge-Kutta or linear multistep methods. This approach is satisfactory for many problems and probably the method of choice for general equations. However, for some classes of problems the classical codes are too unspecialized and methods tuned on the characteristic features of a problem receive better results. In this project we will construct some specialized methods fo rhighly-oscillatory problems (e.g. Sturm-Liouville problems).