This paper is devoted to construction of a mathematical model combining simplicity for practical usage and high accuracy. It is based on the solution of an axisymmetric contact problem on penetration of a rigid indenter into an elastic half-space with a functionally graded or homogeneous coating. The problem is reduced to solution of a dual integral equation. Asymptotically exact expressions for indentation force, depth, contact stiffness and distribution of contact pressures are obtained in simplified analytical form using one-parameter approximation of the integral equation kernel transform. Numerical calculations are provided for a number of homogeneous and functionally graded coatings. Accuracy of the solution is analyzed against ratio of Young's moduli of coating and substrate and the value of relative coating thickness.