Abstract
The calculation of orientation distribution functions (ODF) from incomplete pole figures of any crystal symmetry can be carried out by an iterative procedure taking into account the positivity condition for all pole figures. During the iterative procedure the pole figures are normalized and the pole figure inversion leads to stable results, i.e. the coefficients for higher degrees of series expansion keep small.