clear;
tic;
liczba_krokow = 20000;
pokaz_co = 1000;
p=0.01;
T = 100;
Tmin = 0.01;
wT = 0.995;
c = 0.1;
r=0.1;
w=0.995;
[X1 X2] = meshgrid(-10:p:10,-10:p:10);
Z = [];
i=1;
j=1;
while i<=length(X1)
    j=1;
    while j<=length(X2)
        Z(i,j)= liczY(X1(i),X2(j));
        j=j+1;
    end
    i=i+1;
end

x10=round((rand()*20-10)*100)/100;
x20=round((rand()*20-10)*100)/100;
f0=liczY(x10,x20);
kro=1;
while T>Tmin
    x1n=round((x10+r*randn())*100)/100;
    x2n=round((x20+r*randn())*100)/100;
    fn=liczY(x1n,x2n);
    if fn<f0
        x10=x1n;
        x20=x2n;
    else if rand < 1/(1+exp(fn-f0/(c*T)))
            x10=x1n;
            x20=x2n;
        end
    end
    T=T*w;
end
mesh(X1,X2,Z);