title: Monopole, Dipole and Quadrupole maps with Healpy
body: <p>I am trying to create a map of the monopole (l=0, m=0), dipoles (l=1, m=-1,0,1) and quadrupoles (l=2, m=-2,-1,0,1,2) of the spherical harmonics using Healpy. So far I constructed this extremely simple code with the help of the healpy tutorial:</p>
  
  <pre><code>import numpy as np
  import healpy as hp
  
  NSIDE=128
  m=np.arange(hp.nside2npix(NSIDE))
  x=hp.pixelfunc.fit_monopole(m, nest=False)
  y=hp.pixelfunc.fit_dipole(m, nest=False)
  </code></pre>
  
  <p>but this just returns the values of the monopole (x) and the value of the monopole and dipole (y):</p>
  
  <pre><code>In [22]: print x
  98303.5
  
  In [23]: print y
  (98303.5, array([ -4.59988602e-01,  -1.81289977e+02,  -9.83040000e+04]))
  </code></pre>
  
  <p>When I try hp.mollview(x) I get the following error:</p>
  
  <pre><code>TypeError: 'float' object is not iterable
  </code></pre>
  
  <p>How can I create maps of the monopole and dipoles. Also how do I create maps of the quadrupoles?</p>