from mpl_toolkits.mplot3d import axes3d
import matplotlib.pyplot as plt

fig = plt.figure()  # fig - figure, yep, logical
ax = fig.add_subplot(111, projection='3d')  # three digit integer, where the first digit is the number of rows,
# the second the number of columns, and the third the index of the subplot.
# Projection - custom projection (rectilinear is default)

# Grab some test data.          (gotta work there, duh)
X, Y, Z = axes3d.get_test_data(0.05)    # return a tuple X, Y, Z with a tested data

# X и Y заданы кортежом от -30 до 29,5 с шагом 0,5 (видимо ещё происходит деление на 0.05, отсюда такие малые значения)
# Z высчитывается благодаря заданной формуле
# Plot a basic wireframe.
ax.plot_wireframe(X, Y, Z, rstride=10, cstride=10)  # rstride and cstride - step size

# т.к. X, Y, Z - кортежи, то делать цикл for с перебором значений здесь не имеет смысла

plt.show()