optionization tool

function [c, ceq] = constrain(x) c(1) = (x(1) - 5) .^ 2 + (x(2) - 5) .^ 2 - 25; c(2) = -(x(1) - 8) .^ 2 - (x(2) + 3) .^ 2 + 7.7; c(3) = -x(1) + 2; ceq = []; end