[sol,m,Q,residuals,everything] = solvesos(F,obj,options,params,candidateMonomials)
sdpvar x t a
p = x^2-a*x; % Check lower bound of this w.r.t a
[F,obj,monoms,etc] = compilesos(sos(p-t),-t,sdpsettings('sos.model',2,'sos.scale',0),[a;t]);
sossolver = optimizer(F,obj,sdpsettings('solver','mosek'),a,{t,etc.BlockedQ{1}{1}});
sol = sossolver{0};
coefficients(replace(p,a,0)-sol{1})-monoms{1}'*sol{2}'*monoms{1},x)