Abstract We study numerically the two-point correlation functions of height in six-vertex model with domain wall boundary conditions. The and are computed by Markov chain Monte-Carlo algorithm. Particular attention is paid to free fermionic point (Δ = 0), for which obtained analytically thermodynamic limit. A good agreement exact numerical results allows us extend calculations disordered (|Δ| &...