[x y] = meshgrid(0:0.01:0.5, 0:0.01:0.4); %网格生成用meshgrid命令
z=-0.013*sin(pi.*x./0.5).*sin(pi.*y./0.4);
surf(x,y,z)