t = sdpvar;
x = t^2+ t + 1;
y = 2*t^2 + t + 2;
% Minimum distance determination version
optimize(1<=t<=2,x^2+y^2,sdpsettings('solver','baron'))

% Feasibility version. [] is used instead of 0 for the objective, which accomplishes the same thing
optimize([1<=t<=2,x^2+y^2-33<=0],[],sdpsettings('solver','baron'))