f( x , y, z) = (const.)
No Mathematica,
ContourPlot3D[
Abs[Sin[\[Pi] x] Sin[\[Pi] y] Sin[\[Pi] z]]^2 == 1/2,
{x, -1, 1}, {y, -1, 1}, {z, -1, 1}]

Mostre as superfícies de probabilidade constante 0.2, 0.5 e 0.8:
ContourPlot3D[
Abs[Sin[\[Pi] x] Sin[\[Pi] y] Sin[\[Pi] z]]^2,
{x, -1, 1}, {y, -1, 1}, {z, -1, 1}, Contours -> {0.2, 0.5, 0.8},
ContourStyle -> (Directive[#, Opacity[0.25]] & /@ {Yellow, Orange, Red}),
Lighting -> "Neutral", Mesh -> None]

Você pode fazer algum tipo de visualização de volume , possivelmente com recortes e fatias. Você poderá atribuir uma cor e uma opacidade a cada ponto em 3D. Ferramentas mais avançadas também permitem escolher uma função de transferência.
imgdata =
Table[Abs[Sin[\[Pi] x] Sin[\[Pi] y] Sin[\[Pi] z]]^2,
{x, -1., 1, .01}, {y, -1., 1, .01}, {z, -1., 1, .01}];
img = Image3D[imgdata, ClipRange -> {{150, 200}, {0, 100}, {0, 200}}]

O fatiamento geralmente ajuda, especialmente se você pode controlar interativamente qual fatia exibir.
Image3DSlices[img, Range[1, 200, 10]]
