sdpvar x y[p,c] = polynomial([x y],4);solvesos([sos(p), replace(p,[x y],[pi 2]) == 3],[],[],c);f = sdisplay(replace(p,c,value(c)));f = f{1};f = replace(f,'*','.*');f = replace(f,'^','.^');f = eval(['@(x,y)' f]);
[X,Y] = meshgrid(-1:0.01:1,-1:0.01:1);F = f(X,Y);
mesh(X,Y,F);