The paper presents analysis of surface induction hardening of a cylindrical element made of steel 38Mn6. The mathematical model of the non-stationary process is elaborated. Calculations of coupled electromagnetic and temperature fields are provided by means of the Flux 2D software. Computations are compared with the measurements realized at the laboratory stand located in the Silesian University of Technology. The expected hardness distribution within the surface layer is noticed.