function plotsphfun
[lats,lons]=meshgrid((-90:3:90)*pi/180,(-180:5:180)*pi/180);
x=cos(lats).*cos(lons); y=cos(lats).*sin(lons); z=sin(lats);
label={};
for n=0:10
    for m=0:n
        s=sprintf('Y%d%d',n,m);
        if m==0, label={label{:},s}; else label={label{:},[s,'c'],[s,'s']}; end
    end
end
figure, k=1;
for i=1:9
    for j=1:9
        axes('position',[(j-1)/9,(9-i)/9,1/9,1/9])
        hold on, axis off, axis equal, view(0,30); lightangle(90,40);
        enm=zeros(1,121); enm(k)=1;
        for n=1:size(lats,1)
            for m=1:size(lats,2)
                v(n,m)=tecmodel(0,[lats(n,m),lons(n,m)],enm,1);
            end
        end
        surf(x,y,z,v,'edgecolor','none','facecolor','interp','facelighting','phong',...
             'diffusestrength',0.3,'ambientstrength',0.6,'specularstrength',0.5)
        text(-1.2,1,0.7,label(k)); k=k+1;
    end
end

