;
0022 npix = 12*nside*nside;
0023 Spin = 2;
0024 
0025 
0026 TT_l = rand(L,1);
0027 EE_l = rand(L,1);
0028 TE_l = rand(L,1);
0029 BB_l = rand(L,1);
0030 <span class="keyword">for</span> l = 0:Spin-1
0031     TT_l(l+1) = 0;
0032     EE_l(l+1) = 0;
0033     TE_l(l+1) = 0;
0034     BB_l(l+1) = 0;
0035 <span class="keyword">end</span>
0036 [T_lm, E_lm, B_lm] = ebsep_sim_maps(TT_l, EE_l, TE_l, BB_l, L);
0037 [QpU_lm] = ebsep_compute_QpU_fullsky(E_lm, B_lm);
0038 [QmU_lm] = ebsep_compute_QmU_fullsky(E_lm, B_lm);