怎么用matlab画以下函数的三维图（以x,y,z为坐标）

x=1:0.1:10;y=1:0.1:10;z=0:0.1:10;
[x,y,z]=meshgrid(x,y,z);
f=((x-y).^2-1).*((y-z).^2-1).*((z-x).^2-1);
p=patch(isosurface(x,y,z,f,0));
set(p, 'FaceColor'橡仔, 'red', 'EdgeColor', 'none'余团);
daspect([1 1 1])
view(3)
camlight; lighting phong

圆柱面的话，方程有问题，应该竖如橘是((x-y)^2-1)+((y-z)^2-1)+((z-x)^2-1)=0吧？

x=-3:0.1:3;y=-3:0.1:3;z=-3:0.1:3;
[x,y,z]=meshgrid(x,y,z);
f=((x-y).^2-1)+((y-z).^2-1)+((z-x).^2-1);
p=patch(isosurface(x,y,z,f,0));
set(p, 'FaceColor', 'red', 'EdgeColor', 'none');
daspect([1 1 1])
view(3)
camlight; lighting phong