We give a description of the boundary singularity of the Szegö kernel of M-harmonic functions, i.e. functions annihilated by the invariant Laplacian, on the unit ball of the complex n-space, in terms of the Gauss hypergeometric functions.