import numpy as np
from mpl_toolkits.mplot3d import Axes3D
x = np.arange(5, 7, 0.1)
d = (600 + 1000 * x) / 3600
d
array([ 1.55555556, 1.58333333, 1.61111111, 1.63888889, 1.66666667, 1.69444444, 1.72222222, 1.75 , 1.77777778, 1.80555556, 1.83333333, 1.86111111, 1.88888889, 1.91666667, 1.94444444, 1.97222222, 2. , 2.02777778, 2.05555556, 2.08333333])
%pylab --no-import-all
Using matplotlib backend: Qt4Agg Populating the interactive namespace from numpy and matplotlib
x = np.arange(5, 20, 0.1)
Ns = 1000. / np.arange(1, 11)
X, Y = np.meshgrid(x, Ns)
Z = (5 + X) * Y / 3600
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot_surface(X, Y, Z)
<mpl_toolkits.mplot3d.art3d.Poly3DCollection at 0x45a3550>