denom_f := (X[1]^2+X[2]^2+X[3]^2)^2;
motzkin_diff := X[2]^2*X[1]^4-27*X[1]^2*X[2]^2*X[3]^2+9*X[1]^2*X[2]^4+81*X[3]^
6;
vec_d := Vector(34, [0,0,0,0,0,0,1/324,2182758400/284650799,831180333080/
3308566689,0,0,0,1/324,2182758400/388911398639,1135621284025880/4397949999,0,0
,0,0,0,0,0,1/81,0,0,0,0,0,0,0,0,0,0,0]);
q_list := Vector(34, [0,0,0,0,0,0,-687563/23360*X[1]^3*X[2]^2-1893231/23360*X[
1]*X[2]^2*X[3]^2+256329/11680*X[1]*X[2]^4+324*X[1]*X[3]^4,-2444707357/
19644825600*X[1]^3*X[2]^2+284650799/2182758400*X[1]*X[2]^2*X[3]^2+1080028279/
1091379200*X[1]*X[2]^4,-367618521/831180333080*X[1]^3*X[2]^2+3308566689/
831180333080*X[1]*X[2]^4,0,0,0,-314289/23360*X[1]^2*X[2]*X[3]^2-63459/11680*X[
1]^4*X[2]+615591/23360*X[1]^2*X[2]^3+324*X[2]*X[3]^4,388911398639/2182758400*X
[1]^2*X[2]*X[3]^2-13888467971/1091379200*X[1]^4*X[2]-138918975161/2182758400*X
[1]^2*X[2]^3,-488661111/1135621284025880*X[1]^4*X[2]+4397949999/
1135621284025880*X[1]^2*X[2]^3,0,0,0,0,0,0,0,81*X[1]^2*X[2]^2*X[3]+81*X[3]^5-\
81*X[1]^2*X[3]^3-81*X[2]^2*X[3]^3,0,0,0,0,0,0,0,0,0,0,0]);
t_newton := .47e-1;