clc;
clear;
%定义球半径
R_t=10^4;
n=100;
fai=-pi/2:(pi/2)/n:pi/2;
theta=0:2*pi/n:2*pi;
[theta,fai]=meshgrid(theta,fai);
R_t=R_t.*cos(fai);
[x,y,z]=sph2cart(theta,fai,R_t);
figure(2)
surf(x,y,z)
clear;
%定义球半径
R_t=10^4;
n=100;
fai=-pi/2:(pi/2)/n:pi/2;
theta=0:2*pi/n:2*pi;
[theta,fai]=meshgrid(theta,fai);
R_t=R_t.*cos(fai);
[x,y,z]=sph2cart(theta,fai,R_t);
figure(2)
surf(x,y,z)