分别用surf和mesh函数绘制
%先计算极限
syms x1 y1
f=sin( sqrt(x1.^2+y1.^2) ) ./sqrt(x1.^2+y1.^2);
L=limit(limit(f,x1,0),y1,0);
[x,y]=meshgrid(-8:0.1:8,-8:0.1:8);
z=sin( sqrt(x.^2+y.^2) ) ./sqrt(x.^2+y.^2);
z(isnan(z))=double(L);
figure(1);
surf(x,y,z);
figure(2);
mesh(x,y,z);