“…Forty years ago, Wu, Barouch, McCoy and Tracy [4] showed that the full susceptibility of the square-lattice Ising model can be decomposed as the infinite sum of holonomic n-fold integrals [5,6,7,8,9], denoted χ (n) . In the last decade the linear differential operators corresponding to the first χ (n) 's, up to n = 6, were obtained, underlying the role of the elliptic curve parametrization [10], but showing also the emergence of (at least) one Calabi-Yau ODE, and beyond, of linear differential operators with selected differential Galois groups [11,12,13]. A complete description of the singular points of the linear differential operators corresponding to the first few χ (n) 's has also been obtained [6,14,15,16].…”