from mpl_toolkits.mplot3d import Axes3D from matplotlib import pyplot as plt import numpy as np from skimage import measure # Set up mesh n = 100 x = np.linspace(-3, 3, n) y = np.linspace(-3, 3, n) z = np.linspace(-3, 3, n) X, Y, Z = np.meshgrid(x, y, z) # Create cardioid function def f_heart(x, y, z): F = (2*x**2 + 2*y**2 + z**2 - 1)**3 - 1/10*x**2*z**3 - y**2*z**3 return F # Obtain value at ev..