Is there a function which returns one of the properties accessible through cubeprop() (e.g. electron density), but at any arbitrary 3D coordinate, given a wavefunction? I assume the information itself could be easily computed, given that cubeprop() can return this property if you choose the grid correctly, but this of course comes along with the entire cubefile. I’m thinking of something like this:

```E, wfn = psi4.energy(level_of_theory, return_wfn=True)
result = pointprop(wfn, [x, y, z], tasks="density")```

oeprop() seems to provide some of this functionality, but not for instance electron density.

While oeprop can calculate the ESP at user defined grid points (http://psicode.org/psi4manual/master/oeprop.html#sec-oeprop-grid) there is no function to just give the density.
Though I imagine that it would be small task to implement given the ESP routine uses the density.