#you probably need to type "dev x11" to get a plotting window
#use "dev postencap 'filename'" to make .ps files instead 
#then you need to type "data cmbfast_output.cl"
#then type "input plot_cmbfast.sm"
xcolumn 1 
ycolumn 2
limits 0 1200. 0 120.
box
set y=sqrt(y)*2.73E6  
#turns unitless fluctuations into deltaT in uK.
connect x y 
xlabel  \it{l}
ylabel ({\it{l}}({\it{l}}+1)C_{\it{l}}/2Pi)^(1/2)  {\mu}K
#that's it, enjoy the plot!  ;)
#you can plot multiple theory curves by repeating these steps 
#type "erase" to get a fresh window