docs: def ellipse(n=1000, adaptive=False):
      """
      Get a parameterized set of vectors defining
      ellipse for a major and minor axis length.
      Resulting vector bundle has major axes
      along axes given.
      """
      u = N.linspace(0,2*N.pi,n)
      # Get a bundle of vectors defining
      # a full rotation around the unit circle
      return N.array([N.cos(u),N.sin(u)]).T