linspace(a,b,c)均匀生成介于a到b的c个值，c默认为100
如linspace(0,100,5)即[0 25 50 75 100]
linspace(-1.3,1.3)就是生成介于-1.3到1.3的100个值
meshgrid是生成网格的函数，一般是生成二维网格矩阵数据，但这里生成的是三维网格矩阵数据。
先以二维数据为例，
比如要绘制函数z=f(x,y)的图像
x方向在[1,2,3,4,5]五个点上采样，即x=1:5
y方向在[10 20 30 40]四个点上采样，即y=10:10:40
那么就一共有5*4=20个采样点。
分别为
(1,10) (2,10) (3,10) (4,10) (5,10)
(1,20) (2,20) (3,20) (4,20) (5,20)
(1,30) (2,30) (3,30) (4,30) (5,30)
(1,40) (2,40) (3,40) (4,40) (5,40)
命令
[X,Y]=meshgrid(x,y)
生成的X就是这20个采样点的横坐标，Y既是其纵坐标，即
X=
1 2 3 4 5
1 2 3 4 5
1 2 3 4 5
1 2 3 4 5
Y=
10 10 10 10 10
20 20 20 20 20
30 30 30 30 30
40 40 40 40 40
然后再对这20个采样点计算z的值，即Z=f(X,Y)，最后用mesh(X,Y,Z)即可绘制图像
三维网格数据与此类似，现在xyz三个方向都是取-1.3到1.3的100个数据，总共就是100^3即一百万个采样点，[x,y,z]=meshgrid(linspace(-1.3,1.3));中的x，y和z就分别是这一百万个采样点的三个维度上的坐标。
而
s=(x.^2 + (9/4)*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (1/9)*y.^2.*z.^3;
就得到了函数s=f(x,y,z)=(x.^2 + (9/4)*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (1/9)*y.^2.*z.^3在这一百万个采样点中的值
