Scaling functions in the square Ising model

We show and give the linear differential operators ${\cal L}^{scal}_q$ of order q= n^2/4+n+7/8+(-1)^n/8, for the integrals $I_n(r)$ which appear in the two-point correlation scaling function of Ising model $ F_{\pm}(r)= \lim_{scaling} {\cal M}_{\pm}^{-2} < \sigma_{0,0} \, \sigma_{M,N}> = \sum_{n} I_{n}(r)$. The integrals $ I_{n}(r)$ are given in expansion around r= 0 in the basis of the formal solutions of $\, {\cal L}^{scal}_q$ with transcendental combination coefficients. We find that the expression $ r^{1/4}\,\exp(r^2/8)$ is a solution of the Painlev\'e VI equation in the scaling limit. Combinations of the (analytic at $ r= 0$) solutions of $ {\cal L}^{scal}_q$ sum to $ \exp(r^2/8)$. We show that the expression $ r^{1/4} \exp(r^2/8)$ is the scaling limit of the correlation function $ C(N, N)$ and $ C(N, N+1)$. The differential Galois groups of the factors occurring in the operators $ {\cal L}^{scal}_q$ are given.