iFit/sw_orbital

 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.