先用meshgird 将x和y扩成89*89的矩阵
[x,y]=meshgrid(1/180*pi:1/180*pi:89/180*pi);
然后z=f(x,y);就可以了
注意f(x,y)你面的向量乘法,除发,乘方要用 .* ./ .^2