%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第五章例5-2-3程序
% 空間二次曲面的繪制
% 西安電子科技大學出版社出版  陳懷琛編著 2007年5月
%
clear
a=input('a=');b=input('b=');	%輸入參數，N為網格線的數目
c=input('c=');N=input('N=');
xgrid= linspace(-abs(a),abs(a),N);		% 建立x網格坐標
ygrid= linspace(-abs(b),abs(b),N);		% 建立y網格坐標
[x,y]=meshgrid(xgrid,ygrid);
z=c*sqrt(1-y.*y/b/b-x.*x/a/a);u=1;		% u=1,表示z要取正負值
z1=real(z);										% 取z的實部z1
for k=2:N-1							% 以下七行程序為取消z中含虛數的點
   for j=2:N-1
      if imag(z(k,j))~=0 z1(k,j)=0; end
      if all(imag(z([k-1:k+1],[j-1:j+1])))~=0 z1(k,j)=NaN; 			
      end
   end
end
surf(x,y,z1),hold on				% 畫空間曲面
if u==1 z2=-z1;surf(x,y,z2);	% u=1時加畫負半面，并加負坐標軸
	axis([-abs(a),abs(a),-abs(b),abs(b),-abs(c),abs(c)]);
end
xlabel('x'),ylabel('y'),zlabel('z')
axis equal,hold off,
set(gcf,'color','w')