%General Case -- Caprasse
clc;
clear;
format long;
format infsup;

syms x1 x2 x3 x4
f1=-x1^3*x3+4*x1*x2^2*x3+4*x1^2*x2*x4+2*x2^3*x4+4*x1^2-10*x2^2+4*x1*x3-10*x2*x4+2;
f2=-x1*x3^3+4*x2*x3^2*x4+4*x1*x3*x4^2+2*x2*x4^3+4*x1*x3+4*x3^2-10*x2*x4-10*x4^2+2;
f3=x2^2*x3+2*x1*x2*x4-2*x1-x3;
f4=2*x2*x3*x4+x1*x4^2-x1-2*x3;
F=[f1,f2,f3,f4]
vars=[x1,x2,x3,x4];
zeros=rand(1,4)/1000+rand(1,4)/1000*sqrt(-1)+[2,-sqrt(-3),2,sqrt(-3)];
appsol=zeros.'
zeros=zeros+rand(1,4)*sqrt(-1)/1000;
tol= 0.02;

[X,xs] = viss(F,vars,zeros,tol)

