37,716
社区成员
发帖
与我相关
我的任务
分享
# -*- coding: utf-8 -*-
from numpy import *
from enthought.mayavi import mlab
#当科学以直观形式出现在我眼前的时候,我被惊呆了
# Create the data.
dphi, dtheta = pi/150, pi/150
[phi,theta] = mgrid[0:pi+dphi*1.5:dphi,0:2*pi+dtheta*1.5:dtheta]
m0 = 4; m1 = 3; m2 = 2; m3 = 3; m4 = 6; m5 = 2; m6 = 6; m7 = 4;
r = sin(m0*phi)**m1 + cos(m2*phi)**m3 + sin(m4*theta)**m5 + cos(m6*theta)**m7
x = r*sin(phi)*cos(theta)
y = r*cos(phi)
z = r*sin(phi)*sin(theta)
# View it.
s = mlab.mesh(x, y, z, representation="wireframe", line_width=1.0)
mlab.show()