P = defipoly({'min z', ...
      '71*x1^3 + 23*x2^2 - 25*x3^3 + 67*x1^2*x2 + 5*x1^2*x3 - 13*x2^2*x1 + 75*x2^2*x1 + 14*x2^2*x3 - 19*x3^2*x1 + 98*x3^2*x2 + 64*x1^2 - 87*x2^2 + 45*x3^2 + 80*x1*x2 - 92*x1*x3 + 46*x2*x3 + 10*x1 + 5*x2 + 8*x3 + 15 <= z', ...
      '22*x1^3 - 9*x2^2 + 17*x3^3 + 29*x1^2*x2 + 1*x1^2*x3 + 7*x2^2*x1 - 9*x2^2*x1 - 22*x2^2*x3 + 81 *x3^2*x1 - 27*x3^2*x2 + 99*x1^2 - 20*x2^2 + 69*x3^2 + 23*x1*x2 + 71*x1*x3 + 72*x2*x3 + 76*x1 + 98*x2 - 12*x3 + 23 <= z', ...
      '10*x1^3 - 18*x2^2 - 27*x3^3 + 69*x1^2*x2 + 16*x1^2*x3 + 78*x2^2*x1 - 43*x2^2*x1 + 41*x2^2*x3 - 39*x3^2*x1 - 87*x3^2*x2 - 34*x1^2 + 15*x2^2 + 67*x3^2 + 29*x1*x2 + 84*x1*x3 - 45*x2*x3 + 67*x1 +34*x2 - 68*x3 - 58 <= z', ...
     'x1^2 + x2^2 + x3^2 + z^2 <= 100'}, 'x1,x2,x3,z');
output=gloptipoly(P,3);
output
output.sol{:}