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

#!/usr/bin/python# -*- coding: utf-8 -*-import mathclass 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, zdef 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 += azimuthalfor point in splot(input('')):  print("%f %f %f" % point)