In this present study, Box-Behnken response design was employed to optimize and investigate the effect of process variables such as electrolyte dose, stirring speed, and treatment time on total solids (TS), and chemical oxygen demand (COD) removal in electrocoagulation process (EC) to treat grey wastewater using iron electrodes. Regression analysis showed good fit of the experimental data to the second-order polynomial model with coefficient of determination (R 2 ) value greater than 0.95 for both TS and COD removal. Threedimensional response surface plots were generated from mathematical models in order to study the interactive effect of process variables on responses. Under the optimum operating conditions such as electrolyte dose of 1.4 g/l, stirring speed of 300 rpm, and treatment time of 8 min show the predicted removal efficiencies (96.56 and 93.60% for TS and COD removal, respectively) which were close to the experimental values.