frommpl_toolkits.mplot3d import Axes3Dfrommatplotlib import cmfrommatplotlib.ticker import LinearLocatorfrommatplotlib.font_manager import FontProperties
font= FontProperties(fname=r"c:\windows\fonts\simsun.ttc", size=14)
import matplotlib.pyplotasplt
import numpyasnp
fig=plt.figure()
ax= fig.gca(projection='3d')
[x , t]= np.meshgrid(np.array(range(25))/24.0,np.arange(0, 575.5,0.5)/575*17*np.pi - 2*np.pi)
p= np.pi/2 * np.exp(-t/(8*np.pi))
u= 1 - (1 - np.mod(3.6*t, 2 * np.pi)/np.pi) **4/2y= 2 * (x ** 2-x)**2*np.sin(p)
r= u * (x*np.sin(p) + y *np.cos(p))
surf= ax.plot_surface(r * np.cos(t),r * np.sin(t), u *(x*np.cos(p)-y*np.sin(p)),\
rstride= 1,cstride=1,cmap = cm.gist_heat,linewidth=0,antialiased=True)
plt.title(u'祝愿单身尽早脱单',fontproperties=font)
ax.set_xticks([])
ax.set_yticks([])
ax.set_zticks([])
plt.show()