isosurface函数可以画隐函数的等值曲线或曲面
如画笛卡尔心形函数，其方程为（x^2 + (9/4)y^2 + z^2 - 1)^3 - x^2z^3 - (9/80)y^2z^3 = 0
由于是隐函数，不好用一般的方法画出其图形来，这就可以用isosurface来实现，代码如下
figure
[x,y,z]=meshgrid(linspace(-3,3)); %做出网格meshgrid
p=(x.^2+(9/4)*y.^2+z.^2-1).^3-x.^2.*z.^3-(9/80)*y.^2.*z.^3; %实现结果的表达
isosurface(x,y,z,p,0,x);
axis equal;
axis off;
view(0,45); %视角的控制
colormap([1 0 0]); %绘图颜色红色
brighten(0.5); %增亮
camlight right; %光源位置
lighting phong; %光照模式
代码参考http://blog.sina.com.cn/s/blog_4b013fb10100llxh.html（待看）
另外这个函数还有一个好处是，可以得到等值面的顶点和面，然后可以直接调用patch画出来
fv = isosurface(x,y,z,p,0);
patch(fv)
这一点很好，可以将得到的顶点和面保存到文件，用opengl画出来，可以更好的渲染图形，这点已经实现。
下面简单翻译一下help里的解释吧：
isosurface函数用于由体积数据中提取等值数据。
使用方法如下
fv = isosurface(X,Y,Z,V,isovalue)
fv = isosurface(V,isovalue)
fvc = isosurface(...,colors)
[f,v] = isosurface(...)
[f,v,c] = isosurface(...)
isosurface(...)
参数意义为，从由X,Y,Z构造的体积V数据中提取由isovalue指定的等值数据，返回结果fv为一个结构体，包含了等值面的顶点和面（顶点的次序），这些参数可以直接传给patch命令画出图形来。
另外，可以在最后加上colors，返回值中也会得到每个面的颜色插值，这样能够使你用不同于计算等值面的数据控制等值面的颜色映射。
如果不设置返回值，就会自动在当前坐标系里用计算得到的顶点和面创建三维patch对象。
