The paper deals with numerical and experimental solution of residual stresses, which occurs due to a contact loading. The first solved case corresponds to frictionless contact in an indent test accepted from literature. In the finite element analysis, three material options with isotropic, kinematic and combined hardening were tested. The model of Calloch and Marquis, which was implemented into the Ansys FE program by user subroutine, shows the best correlation with experiments. The experimental results were obtained by two experimental methods, namely by the Neutron diffraction method and by the Contour method, which is not well known in the Czech Republic. Hence, this paper includes also a brief description of the new destructive method. Both methods can lead to prediction of full field residual stresses distribution. Some difficulties from application of the contour method to own experiment in rolling contact fatigue domain are also discussed.