clc;
open('Test_2.slx')


opts = cmaes;
opts.TolFun = 1e-3;
% opts.TolFun = 1e-9;

opts.UBounds = [
    500 ... % k2
    100 ... % d2
    ]';

opts.LBounds = [
    10 ... % k2
    1 ...   %d2
    ]'; 

%% Optimierungsparameter

x0 = ones(size(opts.UBounds));
mysigma = 0.33 * (opts.UBounds - opts.LBounds);

%% Optimierungsaufruf

warning off
tic
[optim_params, fmin, counteval, stopflag, out, bestever] = cmaes('f_model', x0, mysigma, opts);
toc
warning on

save('my_optim_params', 'optim_params')%, 'myseed');
fprintf('k= %f; d= %f \n',optim_params(1),optim_params(2))
disp('Fertig!')