[xintaolv]

a=linspace(-3,12.1,2000);
b=linspace(4.1,5.8,2000);
[x,y]=meshgrid(a,b);%将x,y转换为矩阵形式，网状结构
z=+21.5+x.*sin(4*pi*x)+y.*sin(20*pi*y);
surf(x,y,z);%绘制3D彩色表面
close
t=max(max(z));%求最大值
[xx,yy]=max(max(z));
[xxx,yyy]=max(z(:,yy));
t
x(yy,yyy)
y(yy,yyy) %求对应的x,y值