denom_f := 21*X[3]^2+12*X[2]^2+4*X[1]^2;
motzkin_diff := 9*X[2]^2*X[1]^4+81*X[1]^2*X[2]^4-243*X[1]^2*X[2]^2*X[3]^2+729*
X[3]^6;
iters_num := 19;
t_newton := .304;
vec_d := Vector(9, [1/15309,1/8748,1/2916,1/6075,0,100/26973,0,0,0]);
q_list := Vector(9, [15309*X[3]^4-1701*X[1]^2*X[2]^2,8748*X[2]*X[3]^3-972*X[1]
^2*X[2]*X[3],2916*X[1]*X[3]^3-2916*X[1]*X[2]^2*X[3],6075*X[1]*X[2]*X[3]^2-4131
/2*X[1]*X[2]^3-891/2*X[1]^3*X[2],0,26973/100*X[1]*X[2]^3-2997/100*X[1]^3*X[2],
0,0,0]);