meshgrid(0:10,0:20,0:1);你的变量这么少能画出图来嘛?给你个提示
[x,y,z]=meshgrid(0:0.1:10,0:0.1:20,0:0.01:1);
v=(x+(1-0.2).*y.*x).*20.*z.*(1-z).^(20-1)+(1-z).^20-1;
p = patch(isosurface(x, y, z, v, 'b'));
>> isonormals(x,y,z,v, p)
set(p, 'FaceColor', 'red', 'EdgeColor', 'none');
daspect([1 1 1])
view(3)
camlight; lighting phong