We explore the pattern of frequency-dependent linear and non-linear optical (NLO) response of electron impurity doped quantum dots harmonically confined in two dimensions. For some fixed values of transverse magnetic field strength (ωc), and harmonic confinement potential (ω0), the influence of impurity location (x0,y0) on the diagonal components of frequency dependent linear (αxx and αyy), and the first (βxxx and βyyy) NLO responses of the dot is computed through linear variational route. The non-linear polarizabilities undergo maximization at some typical oscillation frequency of the external field depending upon the impurity location.