ObjectiveFunction = @simple_objective;
X0 = [0.5 0.5];   % Starting point
[x,fval,exitFlag,output] = simulannealbnd(ObjectiveFunction,X0)

%%
X1 = x; %second iteration
options = saoptimset('PlotFcns',{@saplotbestx,@saplotbestf,@saplotx,@saplotf}, ...
'TemperatureFcn',@temperaturefast,'TolFun',1e-6,'MaxFunEvals',1e40,'MaxIter',10000, ...
'ReannealInterval',100,'InitialTemperature',5000);
[x,fval,exitFlag,output] = simulannealbnd(ObjectiveFunction,X1,[],[],options)
