[x1,x2]=meshgrid(-3:.1:3);  % 生成網格型矩陣
z=-x1.^2-x2;                % 計算出矩陣上各點的高度

i=find(x1.^2+x2.^2>9); z(i)=NaN;  % 找出 $x_1^2+x_2^2>9$ 的坐標，并置函數值為 NaN
i=find(x1+x2>1); z(i)=NaN;  % 找出 $x_1+x_2>1$ 的坐標，并置函數值為 NaN
surf(x1,x2,z); shading interp;