X=-10:10;
Y=X;
[x y]=meshgrid(X,Y);
z=pr07_fce(x,y);
surf(x,y,z);
max1=fminsearch('pr08_fce',[-10 -10])