这个问题可以这样解决，增加一个meshgrid（）函数，即
x=linspace(0,3,6)；
y=linspace(0,3,6)；
[x,y]=meshgrid（x,y）
z=4*(x-1).^2-(y+1).^2
运行结果
