File Applications/SpinW/m_files/sw_orbital.m

  returns a polygon of selected hydrogen orbitals
 
  hSurf = SW_ORBITAL(qNum, 'Option1', Value1, ...)
 
  Input:
 
  qNum      Quantum numbers in vector (n, l, m, real). If the length of the
            vector is 4, real orbitals will be produced, if real==1, the
            sum of the +m and -m orbitals is calculated, however if
            real==2, the difference is calculated.
 
  Options:
 
  surfRes   Number of points in the surface mesh along every axis,
            default is 30.
  rLim      Limits of the axes, default is 32.
  P         Constant probability surface, default is 1E-5.
  rBohr     Bohr radius, default is 1.
  norm      Whether to normalise the axes, default is true.
 
  Radius is normalised to 1.