In view of the great importance of determination of hydrogen peroxide in many chemical, biological, clinical, environmental processes, and the attractive operational simplicity of potentiometric approach for the enzymatic assay, after preliminary optimization, the efficient response surface methodology (RSM) with a central composite rotatable design (CCRD), was used for modeling and optimization of the initial-rate potentiometric method of HRP enzyme catalyzed the assay of H 2 O 2 . The combined “OVAT” (one-variable-a-time), full factorial and RSM analysis was able to explain the importance of the factors, their interactions, along with their optimum values (i.e. acetate buffer concentration 300mM, 4-fluorophenol concentration=3mM, pH=6.15, enzyme activity=0.95U/ml and T=41°C). The obtained results showed that the determined second-order polynomial equation explains adequately the non-linear nature of the modeled response, as confirmed by a reasonable coefficient of determination (R 2 =78.3%). Accordingly, the ANOVA (analysis of variance) indicated, particularly, that the terms x 1 =pH and x 2 2 =pH*pH are highly significant (with p value <0.05) in this model. The performance of the optimized method including its linearity range (1–200μM), within-day reproducibility of H 2 O 2 detection in aqueous samples at 2 levels (RSD%=1.06 at 17μmoll −1 , and RSD%=1.81 at 34μmoll −1 , N=10 replica), along with recovery (or matrix effect) of the method (97.71–03.84%) in hair bleaching real samples, were all determined. The optimized method was also successfully validated by the assay of H 2 O 2 in commercial hair bleaching samples and compared to the corresponding values obtained by the reference spectrophotometric method.