clear
clc
function [f]=fst(x)
//f = -sin(0.2*x(1)*x(2))-sin(x(2));
//f = -sin(0.2*x(1)*x(2))-4*sin(x(2)).*cos(x(2))+0.1*x(1)*x(2);
f = -sin(0.2*x(1)*x(2))-sin(x(2))-cos(x(2));
endfunction

function f=f1(x,y)
z =[0;0];
z(1) = x;
z(2) = y;
f = fst(z);
endfunction

x=0:0.1:10;
y=0:0.1:10;
xset('window',1);
xbasc();
subplot(2,1,1);
fplot3d(x,y,f1);
subplot(2,1,2);
xset("fpf"," ");
fcontour2d(x,y,f1,20);

[min_wartosc,wspolrzedne]=optim(list(NDcost,fst),[8;8])
suma_wspolrzednych=sum(wspolrzedne)

subplot(2,1,2);
fcontour2d(x,y,f1,[-1.99, -1.99], style=[-1,-1]);