谁来解释下这个matlab程序中的算法

谁来解释下这个matlab程序中的算法
clf
R = 100;   % 球的半径
n = 21;    % 圆的数量
N = 30;    % 每个平面圆的点数
t = linspace(0,2*pi,N);        %这个lispace什么意思?
for el = linspace(-pi/2, pi/2, n)
    r = R * cos(el); %半径r怎么是按Rcos（el）变化的?  
    x = r * cos(t);   %下面的xyz坐标如何多加几个,比如0,25,50,75,100
    y = r * sin(t);
    z = R * sin(el) * ones( size(x) );
    plot3(x, y, z, 'k');
    hold on
end
axis equal
结果如下：