%《MATLAB及其在理工課程中的應(yīng)用指南》(十一五規(guī)劃版）第五章例5－4－5a
% 畫積分區(qū)域的程序
% 西安電子科技大學(xué)出版  陳懷琛著  2007年7月
%
%本程序給出由x=1,y=x,z=xy三個曲面以及z=0坐標(biāo)平面圍成的積分區(qū)域.
clear, close all
[x,y]=meshgrid(0:.05:1);    % 確定計(jì)算和繪圖的定義域網(wǎng)格
z1=x.*y.*(y-x<0);           % 求z1=xy并去掉定義域外網(wǎng)格區(qū)的z1
mesh(x,y,z1);hold on;       	% 畫出積分區(qū)頂部
% 畫出積分區(qū)域的幾個側(cè)面
x1=[0:0.02:1];y1=x1;sx1=length(x1);	% y=x平面,與z平行
zd=[zeros(1,sx1);x1.*y1];			% x=y平面上,與z平行的線族端點(diǎn)zd=xy
plot3([x1;x1],[y1;y1],zd,'*')		% 畫出此平行線族,端點(diǎn)畫*
plot3(ones(2,sx1),[y1;y1],[zeros(1,sx1);y1],'o') %畫出x=1平行線族端點(diǎn)o
xlabel('x'),ylabel('y'),zlabel('z'),hold off
set(gcf,'color','w')