output of two-particle (coulomb) -integrals

[Michael Vogel]

Hello!
Does anybody know how to calculate and print some of the general 
two-particle-integrals:
                             < a b | c d >
consisting of four orbitals (e.g. natural orbitals) in Molpro?
The documentation does not yield information about this or I did
not find the answer. The Matrop- and Properties-packets seem not to contain
this feature.
If anybody knows this, please help!

Many thanks,
M. Vogel