lmax=6;
lmin=2;
for l=lmin:lmax-1
    for m=0:l
        subplot(lmax-lmin,lmax,lmax*(l-lmin)+m+1)
        sh(l, m, 1);
        title(' ')
    end
end
set(gcf,'Color',[1,1,1])