The aim of this study is the use of kn type factorial programs for the modeling and optimization of the milk coagulation process in the presence of hydrotalcite type clays. By using the experiment programming method, the mathematical model that corresponds to the answer function (the coagulum yield) is obtained, allowing simulating individual or interacting effects between parameters [4]. According to the used experimental program, and the variation of parameters interval, the real optimal values of the parameters can be obtained, in the initial domain or not [7].