The conformal covariance of correlation functions is checked in the second-order transition induced by random bonds in the two-dimensional eight-state Potts model. The decay of correlations is obtained via transfer matrix calculations in a cylinder geometry, and large-scale Monte Carlo simulations provide access to the correlations and the profiles inside a square with free or fixed boundary conditions. In both geometries, conformal transformations constrain the form of the spatial dependence, leading to accurate determinations of the order parameter scaling index, in good agreement with previous independent determinations obtained through standard techniques. The energy density exponent is also computed.PACS numbers: 64.60. Cn, 05.50.+q,05.70.Jk, 64.60.Fr It is well known that quenched randomness can deeply affect the critical properties at second-order phase transitions [1], and is liable to smooth first-order transitions, eventually leading to continuous transitions [2,3].In random systems, owing to strong inhomogeneities inherent in disorder, the usual symmetry properties required by conformal invariance (CI) [4] do not hold. However, by averaging over disorder realizations (denoted by [. . .] av ), translation and rotation invariance are restored, and it is generally believed that conformal invariance techniques should apply in principle. Recent results based on this assumption have recently been obtained at randomness-induced second-order transitions [5,6], but clear evidence of the validity of conformal invariance is still missing. The question is of both fundamental and practical interest. From the fundamental point of view, situations are known where a diverging correlation length does not guarantee the validity of CI. A few years ago, lattice animals were indeed found to be not conformally invariant although they display isotropic critical behavior with correlation lengths satisfying the usual scaling ξ ∼ L with the system size [7]. If conformal invariance works, on the other hand, its powerful techniques might be applied with no restriction to investigate the critical behavior of 2D random systems [8].In the 2D random-bond Ising Model, both analytic [9] and numerical results [10] are available. Transfer matrix (TM) calculations were also used to study the correlation function decay along strips and to compute the conformal anomaly (defined as a universal amplitude in finite-size corrections to the free energy) [11] and, since disorder is marginally irrelevant in the 2D Ising model, conformal invariance techniques were indeed efficient. At randomness-induced second-order phase transitions, a direct comparison between the results deduced from conformal invariance, and standard techniques, such as finitesize scaling (FSS) Monte Carlo (MC) simulations, have nevertheless not yet been made. After the pioneering work of Imry and Wortis, the first large-scale MC simulations devoted to the influence of quenched randomness, in a system whose pure version undergoes a strong firstorder phase transition, were a...