  We investigate second order quasilinear equations of the form f_{ij}
u_{x_ix_j}=0 where u is a function of n independent variables x_1, ..., x_n,
and the coefficients f_{ij} are functions of the first order derivatives
p^1=u_{x_1}, >..., p^n=u_{x_n} only. We demonstrate that the natural
equivalence group of the problem is isomorphic to SL(n+1, R), which acts by
projective transformations on the space P^n with coordinates p^1, ..., p^n. The
coefficient matrix f_{ij} defines on P^n a conformal structure f_{ij} dp^idp^j.
In this paper we concentrate on the case n=3, although some results hold in any
dimension. The necessary and sufficient conditions for the integrability of
such equations by the method of hydrodynamic reductions are derived. These
conditions constitute a complicated over-determined system of PDEs for the
coefficients f_{ij}, which is in involution. We prove that the moduli space of
integrable equations is 20-dimensional. Based on these results, we show that
any equation satisfying the integrability conditions is necessarily
conservative, and possesses a dispersionless Lax pair. Reformulated in
differential-geometric terms, the integrability conditions imply that the
conformal structure f_{ij} dp^idp^j is conformally flat, and possesses an
infinity of 3-conjugate null coordinate systems. Integrable equations provide
an abundance of explicit examples of such conformal structures parametrized by
elementary functions, elliptic functions and modular forms.
