<abstract><p>Problems of non-linear equations to model real-life phenomena have a long history in science and engineering. One of the popular of such non-linear equations is the Duffing equation. An adapted block hybrid numerical integrator that is dependent on a fixed frequency and fixed step length is proposed for the integration of Duffing equations. The stability and convergence of the method are demonstrated; its accuracy and efficiency are also established.</p></abstract>