Python实现将n个点均匀地分布在球面上的方法

#!/usr/bin/python
# -*- coding: utf-8 -*-
import math
class Spherical(object):
 '''球坐标系'''
 def __init__(self, radial = 1.0, polar = 0.0, azimuthal = 0.0):
 self.radial = radial
 self.polar = polar
 self.azimuthal = azimuthal
 def toCartesian(self):
 '''转直角坐标系'''
 r = math.sin(self.azimuthal) * self.radial
 x = math.cos(self.polar) * r
 y = math.sin(self.polar) * r
 z = math.cos(self.azimuthal) * self.radial
 return x, y, z
def splot(limit):
 s = Spherical()
 n = int(math.ceil(math.sqrt((limit - 2) / 4)))
 azimuthal = 0.5 * math.pi / n
 for a in range(-n, n + 1):
 s.polar = 0
 size = (n - abs(a)) * 4 or 1
 polar = 2 * math.pi / size
 for i in range(size):
 yield s.toCartesian()
 s.polar += polar
 s.azimuthal += azimuthal
for point in splot(input('')):
 print("%f %f %f" % point)