A new numerical method of frequency analysis is described, designed mainly to search for discrete frequencies in a time series. An integral transform is applied twice to the data for different reference times. A complex amplitude within a selected narrow frequency band is obtained for each transform. The frequency is then determined from the phase change of the complex amplitude over the difference of the two reference times. Very high precision is obtained, which is demonstrated in two examples.