%这一行没有必要
rE=4226;
e=870;
alphb=linspace(-pi/6,pi/6,11);
betab=linspace(-pi/8,pi/4,10);
[Ab,Bb]=meshgrid(alphb,betab);
x=-(rE*cos(Bb)+e).*sin(Ab);
y=(rE*cos(Bb)+e).*cos(Ab);
z=rE*sin(Bb);
surf(x,y,z) %区分大小写的
%这一行没有必要
rE=4226;
e=870;
alphb=linspace(-pi/6,pi/6,11);
betab=linspace(-pi/8,pi/4,10);
[Ab,Bb]=meshgrid(alphb,betab);
x=-(rE*cos(Bb)+e).*sin(Ab);
y=(rE*cos(Bb)+e).*cos(Ab);
z=rE*sin(Bb);
surf(x,y,z) %区分大小写的