sdpvar x1 x2;
p = 2*x1^4+3*x1^2*x2^2+2*x1*x2^3+4*x2^4;
Fs = sos(p);
[sol,monos,Q] = solvesos(Fs)