《MATLAB及其在理工课程中的应用指南》书上的程序集

				%《MATLAB及其在理工课程中的应用指南》第五章例5－2－2程序
				% 空间曲面的交线
				% 西安电子科技大学出版社出版  陈怀琛编著 1999年10月初版，2004年11月第二版
				%
				% s1=input('输入方程1 ','s');		% z1方程
				% s2=input('输入方程2 ','s');		% z2方程
				% xmax=input('xmax=ymax= ');dx=input('dx=dy= ');
				% [x,y]=meshgrid(-xmax:dx:xmax);
				% z1 = eval(s1);
				% z2 = eval(s2);
				% 若要使输入函数有较大灵活性,将后面的三行换成前面的六行
				[x,y]=meshgrid(-2:.1:2);        % 确定计算和绘图的定义域网格
				z1=x.*x-2*y.*y;	                % 第一个曲面方程
				z2=2*x-3*y;		                % 第二个曲面方程（平面）
				mesh(x,y,z1);hold;mesh(x,y,z2);	% 分别画出两个曲面（保持坐标比例不变）
				set(gcf,'color','w')            % 置图形背景色为白色
				r0=abs(z1-z2)				zz=r0.*z1;yy=r0.*y;xx=r0.*x; 	% 求这些网格上的坐标值,即交线坐标值
				plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0)','x');	% 画出这些点
				colormap(gray),hold off