In this study, a traditional optimization technique Fuzzy Logic and a deterministic optimization technique filled function method are employed in the modeling of surface hardening problems in physics.