Three dimensional calculation of control dynamics for finding the optimized laser filed is formulated using an iterative method and time-dependent density functional approach. An appropriate laser pulse is designed to control the desired products in the dissociation of methane molecule. The tailored laser pulse profile, eigenstate distributions and evolution of the efficient occupation numbers are predicted and exact energy levels of this five-atomic molecule is obtained. Dissociation rates of up to 78%, 80%, 90%, and 82% for CH, CH, C and C are achieved. Based on the present approach one can reduce the controlling costs.