Matlab simulation of Duffing non-linear oscillator with harmonic driving force