A new approach to calculation of optical quantities of inhomogeneous layers is presented. In this approach, the Richardson extrapolation is used to improve the accuracy of the method, in which the inhomogeneous layer is approximated by a stack of thin homogeneous layers. The results presented in this paper are based on the assumption that the media are isotropic and the inhomogeneity is along the axis normal to the boundary. The results obtained by the new method are compared with those obtained without the Richardson extrapolation. The Richardson extrapolation brings significant improvement in accuracy, especially if the number of approximating layers is large. Moreover, the method using the Richardson extrapolation proceeds in steps with an error estimate available in each step; thus, the calculation can be stopped when the desired accuracy is reached. The use of the method is illustrated by means of the optical characterization of strongly inhomogeneous film of non-stoichiometric silicon nitride. KEYWORDS ellipsometry, inhomogeneous layers, optical characterization, silicon nitride 1 INTRODUCTION Within the optical characterization and optical synthesis of thin films, inhomogeneous layers play an important role. While there are many works devoted to methods usable in optical characterization and optical synthesis of homogeneous thin films (see, eg, previous studies 1-18 ), less attention has been devoted to the methods enabling us to characterize or perform optical synthesis of inhomogeneous thin films exhibiting refractive index profiles. Therefore, it is necessary to develop new procedures usable in the optics of inhomogeneous layers.In the literature, many procedures for calculating of the optical quantities of inhomogeneous layers have been published, eg, previous studies. 11,[19][20][21][22][23][24][25][26][27][28][29][30] The theoretical procedures can be divided into those providing the exact solution and those providing approximate solution.The exact solution is known only for few special refractive index profiles (eg, linear profile and exponential profile 24 ); thus, the approximate methods are usually used in practice. As more and more sophisticated methods of the optical characterization and optical synthesis are developed, the approximate methods with increased precision and reduced calculation time are required.The general approximate method enabling us to calculate the optical quantities of inhomogeneous layers with arbitrary refractive index profile is based on dividing these layers into a stack of thin homogeneous layers. 19,20,22,24 We will refer to this method as to homogeneous layers approximation (HLA). Unfortunately, to calculate the optical quantities with high accuracy, a very large number of approximating layers is required in this method. In this paper, we will show that the Richardson extrapolation (RE) can be used to increase the accuracy of this method. The RE 31,32 is a mathematical method that can be used to increase the rate of convergence of a series, provided certain criteria ...