import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.axisartist.axislines import SubplotZero

delta = 0.01
x = np.arange(-3, 3, delta)
y = np.arange(-4, 4, delta)
X, Y = np.meshgrid(x, y)
Z = (10 * (X**3)) + (Y**3) - (15 * X * Y) + 1

fig = plt.figure(figsize=(8,8))
ax = SubplotZero(fig, 1, 1, 1)
fig.add_subplot(ax)

ax.contour(X, Y, Z, [0], colors='red')
ax.grid(True)

for border in ['left','right','top','bottom']:
    ax.axis[border].set_visible(False)

ax.axis['xzero'].set_visible(True)
ax.axis['yzero'].set_visible(True)