%Special Case -- Ojika4
clc;
clear;
format long;
format infsup;

syms x y z
f1=x+z*x^3+x*z*y^2-x*z;
f2=10*y-2*y*z*x^2-z*y^3-y*z;
f3=-6*z^2*x^4-3*x^2*y^2*z^2-z^2*x^2+28*z*x^2-3*z^2*y^4+2*z^2*y^2+7*z*y^2+z^2-11*z+10;

vars = [x,y,z];
zeros = rand(1,3)/1000 + [0,0,1];
tol = 0.005;
F =[f1,f2,f3]
appsol=zeros.'

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




