数学建模与数学实验（第3版）课件+matlab

				clear;
				x=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5];
				y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5];
				z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9];
				
				
				[fnodes,a,rnw,b,c]=e01sef(x,y,z);
				
				nx=100;
				cx=linspace(75,200,nx);
				ny=200;
				cy=linspace(-50,150,ny);
				for i=1:ny
				   for j=1:nx
				      [cz(i,j),ifail]=e01sff(x,y,z,rnw,fnodes,cx(j),cy(i));
				   end
				end
				
				meshz(cx,cy,cz),rotate3d
				xlabel('X'),ylabel('Y'),zlabel('Z')
				pause
				
				figure(2),meshz(cx,cy,cz),rotate3d
				xlabel('X'),ylabel('Y'),zlabel('Z')
				hold on
				dx=30:220;
				dy=-90:170;
				dz=ones(length(dx),length(dy))*(-5);
				meshz(dx,dy,dz)
				pause
				
				figure(3),contour(cx,cy,cz,[-5 -5]);grid
				hold on 
				plot(x,y,'+');
				xlabel('X'),ylabel('Y')