>> clear
>> [x,y]=meshgrid(0.1:0.1:5);
>> r=x+y+x.*y;
>> z=((r-(r.^2-4*x.*y)^0.5)./y/2)^0.5;
>> surf(x,y,abs(z))
>> shading interp
>> xlabel('x');
>> ylabel('y');
>> zlabel('z');
注意上面那個z,里面有復數(shù),所以mesh(x,y,z)或者surf(x,y,z)都會有點小問題.你最好再檢查一下,你的z的表達式.
這里做出它的模,就是長度.
哈哈,上面我粗心了,上面是錯的,你用下面的就OK了.
>> clear
>> [x,y]=meshgrid(0.1:0.1:5);
>> r=x+y+x.*y;
>> z=((r-(r.^2-4*x.*y).^0.5)./y/2).^0.5;
>> surf(x,y,z)
>> shading interp
>> xlabel('x');
>> ylabel('y');
>> zlabel('z');