import numpy as np
from mayavi import mlab

x, y, z = np.ogrid[-10:10:20j, -10:10:20j, -10:10:20j]

s = np.sin(x*y*z)/(x*y*z)

mlab.figure(bgcolor=(0.5,0.5,0.5), size=(800, 600))

mlab.contour3d(s)

#mlab.pipeline.volume(mlab.pipeline.scalar_field(s))

#mlab.pipeline.volume(mlab.pipeline.scalar_field(s), vmin=0, vmax=0.8)

#mlab.pipeline.image_plane_widget(mlab.pipeline.scalar_field(s),plane_orientation='x_axes',slice_index=10)

mlab.outline()

mlab.show()