set pointsize 2
set xlabel 'x = cos theta'
set ylabel 'Spherical Harmonic'
set xrange [-1:1]
set yrange [-0.5:1]
set grid
plot 'leg2.dat' u 1:2 title 'sphPlm(2,0)', '' u 1:3 title 'sphPlm(2,1)', '' u 1:4 title 'sphPlm(2,2)', sqrt(5.0/4/pi)*(1.5*x*x-0.5) title 'Y20', -sqrt(15.0/8/pi)*sqrt(1-x*x)*x title 'Y21', sqrt(15.0/2/pi)/4*(1-x*x) title 'Y22'
