matlab绘制球面及柱面（急!）

x=-10:0.1:10;y=x;z=x;
[x,y,z]=meshgrid(x,y,z);
f1=x.^2+y.^2+z.^2-4;
f2=x.^2+y.^2-1;
f3=x.^2+z.^2-1;
f4=z.^2+y.^2-1;
p1=patch(isosurface(x,y,z,f1,0));
set(p1,'FaceColor','b','EdgeColor','none');
p2=patch(isosurface(x,y,z,f2,0));
set(p2,'FaceColor','r','EdgeColor','none');
p3=patch(isosurface(x,y,z,f3,0));
set(p3,'FaceColor','y','EdgeColor','none');
p4=patch(isosurface(x,y,z,f4,0));
set(p4,'FaceColor','h','EdgeColor','none');
daspect([1 1 1])
view(3); axis tight
camlight;
lighting phong
图像出来之后,注意旋转一下图像3D,否则只能看到2维平面图.