The creation of nano-zero-valent iron (nZVI) by implementing plant extracts is an environmentally friendly process. High antioxidant capacity and phenol content indicated the possibility of oak-nZVI synthesis using oak leaf extract as a stable material with minimal agglomeration. Simultaneously removal of Cd and phosphates and Ni and phosphates was optimized by a statistically designed experiment with a definitive screening design. In terms of significance, 4 input parameters on process productivity were monitored: initial metal concentration (1 - 9 mgL-1), initial ion concentration (1- 9 mgL-1), pH value (2 - 10) and oak-nZVI dosage (2 - 16 ml). Phosphate removal efficiency, in the presence of cadmium, was the most influenced by oak-nZVI dose and cadmium concentration, while pH gave three statistically significant interactions with oak-nZVI dosage, cadmium concentration, and phosphate concentration. Phosphate removal efficiency, in the presence of Ni, is the most influenced by nickel concentration, phosphate concentration, pH, two-factor interaction between nZVI dose and Ni concentration, and quadratic phosphate interaction. The process optimization yielded the highest simultaneous removal efficiency of 98.99% and 87.30% for cadmium and phosphate ions, respectively. Also, the highest simultaneous removal efficiency of nickel and phosphates ion was 93.44% for 96.75%, respectively. The optimization process fits into the confidence intervals, confirming the assumption that the adopted regression model describes the process well. This work demonstrated an enormous potential and prosperous application for the Cd(II), Ni(II), and phosphate removal from water matrices.