Simulation of the vibratory motion of the Duffing Oscillator using several Fourth Order Runge Kutta methods