  We demonstrate the applicability of four-body exponentially correlated
functions for the accurate calculations of relativistic effects in lithium-like
atoms and present results for matrix elements of various operators which
involve negative powers of interparticle distances.
