function shili21
h0=figure('toolbar','none',...
    'position',[200 100 450 450],...
    'name','實例21');
[x,y,z]=peaks(30);
subplot(2,1,1)
x=x(1,:);
y=y(:,1);
i=find(y>0.8&y<1.2);
j=find(x>-0.6&x<0.5);
z(i,j)=nan*z(i,j);
surfc(x,y,z)
xlabel('X軸');
ylabel('Y軸');
zlabel('Z軸');
title('Figure1:surfc函數形成的曲面')


subplot(2,1,2)
x=x(1,:);
y=y(:,1);
i=find(y>0.8&y<1.2);
j=find(x>-0.6&x<0.5);
z(i,j)=nan*z(i,j);
surfl(x,y,z)
xlabel('X軸');
ylabel('Y軸');
zlabel('Z軸');
title('Figure2:surfl函數形成的曲面')