%% Harj. 2, teht. 13 mplV0009.m  
f=@(x,y) 1./x + 1./y + sin((x.^2).*y.^2)
close all
%%
%N=200;
N=2*N
x=linspace(1.6,2,N);y=x;
[X,Y]=meshgrid(x,y);
Z=f(X,Y);
%surf(x,y,Z)
%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)
%%
