Self-consistent method to extract non-linearities from pulsating star light curves – I. Combination frequencies Article in Monthly Notices of the Royal Astronomical Society (August 2020)