Try VerifyRealRoots Online
Edit your code:
% Here is an illustrative example of kissing number problem: Kissing[2,2]. % Define the variables: syms x0 x1 x2 x3; % The list of the equations: h1=0, h2=0 h1 = x0^2 + x1^2 - 1; h2 = x2^2 + x3^2 - 1; E = [h1, h2]; % The list of the inequality: f1>0 f1 = x0^2 + x1^2 + x2^2 + x3^2 -2*x0*x2-2*x1*x3-1; I = [f1]; % Construct the polynomial system: K = {E,I}; % Select the solver and initialize the option: option.solver = 'MMCRSolver'; option.order = 3; % Call the main function verifyrealroots(..) [out, vars] = verifyrealroots(K,option);
Result: