Determination of the frequency-amplitude response curves of undamped forced Duffing’s oscillators using an ancient Chinese algorithm