The analysis of x-ray line profile broadening can be considered as solving a series of inverse problems. There are usually two steps-removing the instrumental contribution (deconvolution), and determining the broadening contribution in terms of crystallite size and microstrain. Here we are concerned with quantifying only the size broadening, in terms of the shape and size distributions of the crystallites. We present a method that removes the instrumental broadening and determines the particle size distribution in a single step. The general theoretical framework developed makes it possible to determine the crystallite shape and average dimensions, and to fully quantify these results by also assigning uncertainties to them.In general, there are two approaches that can be adopted. The first assumes functional forms for the size distribution and shape of the crystallites, and applies least squares fitting to determine the parameters defining the size distribution [1,2]. For pragmatic reasons, this approach is often used to ensure numerical stability; however, it is based on an explicit assumption for the crystallite size distribution and does not take into account the non-uniqueness of the solution.The second approach takes into account the nonuniqueness of the problem of determining the size distribution P(D) from the experimental data, by assigning a probability to the solutions and enabling an average solution to be determined from the set of solutions; moreover, it also allows any a priori information and assumptions to be included and tested. This approach is