Calcul numérique de la fonction harmonique sphérique généralisée dotée des propriétés de symétrie du cristal cubique ou hexagonal

The three-dimensional quantitative analysis of the distribution function f(g) of crystalline orientations requires knowledge of the numerical values of the generalized spherical harmonic function endowed with crystal and sample symmetries. A Fortran program has been written in order to compute this function for cubic or hexagonal crystal symmetry, without a limiting assumption about sample symmetry - except Friedel's law which imposes only a centre of symmetry on the sample. The entire utilization of the relations arising from symmetries allows one to limit the amount of calculation and to obtain the results quickly and easily. This program can be used on a 32-bit-word computer of which the compiler does not work with complex numbers.