Duffing equation can describe many important nonlinear physical systems. In this paper a coupled approach based on quasilinearization and Bessel polynomial collocation method has been suggested to solve nonlinear duffing oscillator equation. The nonlinearity in duffing oscillator can be of various variety. This approach is very efficient, stable and reliable to deal with any kind of nonlinearity. Numerical examples demonstrate the validity and applicability of the approach on various types of nonlinear duffing oscillator equation.