In []:

import numpy as np

import mayavi.mlab as mlab

Consider the three equations:

\(y=x^2+\delta\)
\(z=x^2-\delta\)
\(y=z^2+\delta\)

In []:

delta = 0.5

In []:

res = 200j

x, z = np.mgrid[-3:3:res,-3:3:res]
y = x**2 + delta

mlab.mesh(x, y, z)

if 1:
    y, x = np.mgrid[-3:3:res,-3:3:res]
    z = x**2 - delta
    
    mlab.mesh(x, y, z)

if 1:
    x, z = np.mgrid[-3:3:res,-3:3:res]
    y = z**2 + delta
    
    mlab.mesh(x, y, z)

In []:

mlab.orientation_axes()
mlab.show()