程序代码：

clear;
clc;
yalmip;
Cplex;
De = sdpvar(1,1,'full');
x = sdpvar(3,1,'full');
Ld = sdpvar(2,1,'full');
st = [];
st = st + [3*x(1,1)^2 - 2*x(2,1)^3 +3*x(3,1)^2 + 3 - Ld(1,1) * (x(1,1) - x(2,1)^2 +x(3,1)^2 - 20) - Ld(2,1) * (x(1,1)^2 - x(2,1) +x(3,1)^2 ) < De]; 
ops = sdpsettings('solver','cplex');
result = solvesdp(st,De);