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

syms a b c d
f1=a^4-b*c*d;
f2=b^4-a*c*d;
f3=c^4-a*b*d;
f4=d^4-a*b*c;
vars=[a,b,c,d];
zeros=rand(1,4)/1000;
tol= 0.005;
F=[f1,f2,f3,f4]
appsol=zeros.'

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


                             

