The phenomena of bifurcation and chaos in forced vibration systems containing a square nonlinear term x+kx+ω02x-ax2 =μcosωt have been investigated by using a 2-D iterated map and by integrating the differential equation. The critical values of control parameter μ are inversely proportional to a.