To prevent the numerical instabilities in topology optimization, continuous approximation of material distribution (CAMD) is proposed to the homogenization design method (HDM) and the simple isotropic material with penalization (SIMP) method. The continuous FE approximation of design variables including high order elements is applied to the formulation of SIMP method. Numerical examples are presented to compare the efficiency of CAMD both in HDM and SIMP.