%filename simulation = @run_sim; %starting point x0 = [0.1; 0.1; 50e-3; 0.2e-3]; %constraints b1 = [0.05; 0.05; 25e-3; 0.1e-3]; % Lower bound (b1 < x) b2 = [0.2; 0.2; 75e-3; 1e-3 ]; % Upper bound (x < b2) A = [-eye(4); eye(4)]; b = [-b1;b2]; %optimization x_opt = fmincon(simulation, x0, A, b); %show results [cost, out] = run_sim(x_opt,false);