Summary: We present a new structured wavelet algorithm to solve the Ornstein-Zernike integral equation for simple liquids. This algorithm is based on the discrete wavelet transform of radial distribution functions and different low-rank approximations of the obtained convolution matrices. The fundamental properties of wavelet bases such as the interpolation properties and orthogonality are employed to improve the convergence and speed of the algorithm. In order to solve the integral equation we have applied a combined scheme in which the coarse part of the solution is calculated by the use of wavelets and Newton-Raphson algorithm, while the fine part is solved by the direct iteration. Tests have indicated that the proposed procedure is more effective than the conventional method based on hybrid algorithms.