%% Plot and study minimum.
%
f=@(x,y) 1./x + 1./y + sin((x.^2).*y.^2)
close all
%%
N=100;
%N=2*N
x=linspace(1.6,2,N);y=x;
[X,Y]=meshgrid(x,y);
Z=f(X,Y);
surf(x,y,Z);shg
%figure
%[c,h]=contour(x,y,Z,0.08:.01:0.12);clabel(c,h);shg
%%
format long
fmin=min(min(Z))
[i,j]=find(Z==fmin)
hold on
plot3(X(i,j),Y(i,j),fmin,'*r')
