On this page you can calculate the pair correlation function g(r*) of the Lennard Jones fluid by numerically solving the Ornstein-Zernike equation with the Percus-Yevick approximation. The PY-OZ equation is solved with the method of Gillan. The source code of the program can be found in Lee's book.

You can choose two variables: the reduced temperature T* and the reduced density RHO*. The critical values for the Lennard Jones-fluid are: T*=1.35, RHO*=0.35.

The lower the temperature or the higher the density the longer takes the calculation. For T*<1.0 RHO* should be <0.01. For RHO*> 0.1 T* should be > 1.4. If the message "Document contains no data" appears the calculation took too long e.g. it may not converge. You may begin with T*=1.4, RHO*=0.1.