This presentation suggests and solves a 3D model of laser hardening of steel bodies combined with their induction pre-heating. The model consists of two partial differential equations describing the distribution of a nonlinear and nonstationary magnetic field in the system and a nonstationary temperature field in the heated object. Both equations are solved numerically in the hard-coupled formulation, respecting the temperature dependences of the material parameters. The methodology is illustrated with a typical example.