i=0;
a=-0.5:0.02:0.5;
for x=a
  i=i+1; j=0;
  for y=a
    j=j+1;
    z(i,j)=corana([x y]);
  end
end
clf
%x平面內的圖形
plot(z(:,1)) 
xlabel('x')  
%x的范圍為[250.0-250.25]
clf
%y平面內的圖形
plot(z(1,:))
xlabel('y') 
%y的范圍為[0-250]
%三維圖形
mesh(a,a,z);
view(30,60);
grid;