f := -4*a^9*A*y^3+7*a^10*A^4*y^2+16*a^6*A^2*y^2+6*a^4*B^4*x^4+38*a^8*A^2*x^2+
48*a^5*B^3*x-4*a^5*B*x^5-4*a^5*B^3*x^5+a^10*A^4*B^2*x^2+3*a^2*B^4*y^2+8*a^4*A^
2*y^4+16*a^9*A*y+a^12*A^4+a^12*A^6+a^8*y^4+a^6*y^6+a^6*x^6+a^4*x^4-8*a^10*A^4-\
8*a^10*A^2+16*a^10*x^2-8*a^8*A^4+8*a^8*A^2-8*a^8*y^2+32*a^8*x^2+16*a^8*B^2-8*a
^8*x^4+32*a^6*A^2+8*a^6*y^2+8*a^6*x^2+32*a^6*B^2-8*a^6*y^4-8*a^6*x^4+16*a^4*A^
2+32*a^4*y^2-8*a^4*x^2+8*a^4*B^2-8*a^4*y^4+B^4+16*a^8+a^10*A^6*y^2+a^6*A^2*y^6
+a^6*B^2*x^6+a^2*B^6*x^2-4*a^11*A^3*y-4*a^11*A^5*y+16*a^10*B^2*x^2+3*a^10*A^4*
x^2+6*a^10*A^2*y^2+20*a^10*A^2*x^2+8*a^9*A^3*y-8*a^9*A^5*y-4*a^9*A^5*y^3-8*a^9
*A^3*y^3+32*a^9*B^3*x^3+32*a^9*B*x^3+a^6*A^2*x^4*y^2+a^6*B^2*x^2*y^4+a^2*A^2*B
^4*y^2+8*a^11*A^2*B*x+10*a^11*A^4*B*x+32*a^10*A^2*B^2*x^2+14*a^9*A^4*B*x+32*a^
6-2*a^10*A^5*B*x*y-24*a^10*A^3*B*x*y-16*a^10*A*B*x*y+8*a^9*B*x*y^2+10*a^9*A^3*
x^2*y+8*a^9*A^2*B^3*x^3+2*a^9*A^2*B*x^3+8*a^9*A*x^2*y+44*a^8*A^2*B^2*x^2+2*a^8
*A^4*x^2*y^2+32*a^8*B^2*x^2*y^2+16*a^4+B^6+16*a^9*A^4*B*x*y^2-20*a^9*A^3*B^2*x
^2*y+28*a^9*A^2*B*x*y^2+32*a^8*B^2*x^2+8*a^8*A^4*y^2+20*a^8*x^2*y^2+3*a^8*A^4*
B^2+6*a^8*A^4*y^4+8*a^8*B^2*x^4+20*a^8*A^2*B^2+7*a^8*A^2*y^4+3*a^8*A^2*x^4+16*
a^8*B^4*x^4+32*a^7*B*x+48*a^7*A^3*y-8*a^7*A^5*y^3+48*a^7*B^3*x^3+48*a^7*B*x^3+
48*a^7*A*y+8*a^7*A*y^3-4*a^7*A^3*y^5+32*a^7*B^3*x-32*a^9*A*B^2*x^2*y-28*a^8*A^
3*B*x*y-4*a^8*A^3*B*x^3*y-28*a^8*A^3*B*x*y^3+22*a^8*A^2*x^2*y^2+2*a^8*A^2*B^2*
x^4-4*a^7*A^2*B*x+2*a^7*A^3*x^2*y+12*a^7*A^2*B^3*x^3+10*a^7*A^2*B*x^3-4*a^7*A*
x^2*y+10*a^7*A^3*B^2*y+12*a^7*A^3*x^2*y^3+2*a^7*B*x^3*y^2+10*a^7*B*x*y^4+2*a^7
*A^2*B^3*x+8*a^7*A*B^2*y+10*a^7*A*x^2*y^3+14*a^7*A*x^4*y+8*a^7*B^3*x^3*y^2+22*
a^6*A^2*B^2*x^2+44*a^6*B^2*x^2*y^2+44*a^6*A^2*x^2*y^2+2*a^6*A^4*B^2*y^2+22*a^6
*A^2*B^2*y^2+2*a^6*A^2*B^4*x^2+2*a^6*A^2*x^2*y^4+2*a^6*B^2*x^4*y^2+8*a^5*A^2*B
*x-4*a^5*B*x*y^2+2*a^5*A^3*B^2*y+8*a^5*A^3*x^2*y^3+10*a^5*B*x^3*y^2+14*a^5*B*x
*y^4+10*a^5*A^2*B^3*x-4*a^5*A*B^2*y+2*a^5*A*x^2*y^3+10*a^5*A*x^4*y+12*a^5*B^3*
x^3*y^2+12*a^5*A^3*B^2*y^3+2*a^5*B^3*x*y^2+10*a^5*A*B^2*y^3+22*a^4*B^2*x^2*y^2
+32*a^4*A^2*x^2*y^2+44*a^4*A^2*B^2*y^2+2*a^4*A^2*B^2*y^4+2*a^4*B^4*x^2*y^2+8*a
^3*B*x*y^2+8*a^3*A*x^2*y+8*a^3*A^3*B^2*y^3+10*a^3*B^3*x*y^2+2*a^3*A*B^2*y^3+14
*a^3*A*B^4*y+32*a^2*A^2*B^2*y^2+8*a*A*B^2*y+10*a*A*B^4*y+54*a^8*A^2*B^2*x^2*y^
2-12*a^8*A*B*x^3*y-24*a^8*A*B*x*y^3-48*a^8*A*B^3*x^3*y+28*a^7*A^4*B*x*y^2-32*a
^7*A^3*B^2*x^2*y+28*a^7*A^2*B*x*y^2-40*a^7*A*B^2*x^2*y-4*a^7*A*y^5-8*a^7*B^3*x
^5-8*a^7*B*x^5+16*a^6*B^2*x^2+8*a^6*A^4*y^2+20*a^6*A^2*x^2+38*a^6*x^2*y^2+16*a
^6*A^4*y^4+8*a^6*B^2*x^4+38*a^6*A^2*B^2+8*a^6*A^2*y^4+16*a^6*B^4*x^4+20*a^6*B^
2*y^2+8*a^6*B^4*x^2+3*a^6*x^4*y^2+3*a^6*x^2*y^4+48*a^5*B*x+32*a^5*A^3*y+48*a^5
*A^3*y^3+8*a^5*B*x^3-32*a^7*A^2*B*x^3*y^2+16*a^7*A^2*B*x*y^4-20*a^7*A*B^2*x^2*
y^3+28*a^7*A*B^2*x^4*y-12*a^6*A^3*B*x*y+24*a^6*A*B*x*y-56*a^6*A^3*B*x*y^3+80*a
^6*A^2*B^2*x^2*y^2-28*a^6*A*B*x^3*y-28*a^6*A*B*x*y^3-56*a^6*A*B^3*x^3*y-4*a^6*
A^3*B^3*x*y-12*a^6*A*B^3*x*y-2*a^6*A*B*x^5*y-4*a^6*A*B*x^3*y^3-2*a^6*A*B*x*y^5
+32*a^5*A*y+48*a^5*A*y^3-8*a^5*A^3*y^5-8*a^5*A*y^5-8*a^5*B^5*x^3+32*a^4*A^2*y^
2+20*a^4*x^2*y^2+16*a^4*A^4*y^4+7*a^4*B^2*x^4+20*a^4*A^2*B^2+38*a^4*B^2*y^2+8*
a^4*B^4*x^2+3*a^4*B^2*y^4+3*a^4*A^2*B^4-40*a^5*A^2*B*x*y^2+28*a^5*A*B^2*x^2*y-\
20*a^5*A^2*B*x^3*y^2+28*a^5*A^2*B*x*y^4-32*a^5*A*B^2*x^2*y^3+16*a^5*A*B^2*x^4*
y-32*a^5*A^2*B^3*x*y^2+28*a^5*A*B^4*x^2*y-48*a^4*A^3*B*x*y^3+54*a^4*A^2*B^2*x^
2*y^2-24*a^4*A*B*x^3*y-12*a^4*A*B*x*y^3-28*a^4*A*B^3*x^3*y-8*a^4*B^4+16*a^2*y^
2-8*a^2*B^2-8*a^2*B^4-28*a^4*A*B^3*x*y-4*a^4*A*B^3*x*y^3+16*a^3*B*x+32*a^3*A^3
*y^3-8*a^3*B^3*x^3-4*a^3*B*x^3+32*a^3*A*y^3+8*a^3*B^3*x-4*a^3*B^5*x^3-8*a^3*B^
5*x+6*a^2*B^2*x^2+16*a^2*A^2*y^2+20*a^2*B^2*y^2+7*a^2*B^4*x^2-4*a*B^3*x-4*a*B^
5*x-32*a^3*A^2*B*x*y^2+28*a^3*A*B^2*x^2*y-20*a^3*A^2*B^3*x*y^2+16*a^3*A*B^4*x^
2*y-16*a^2*A*B*x*y-24*a^2*A*B^3*x*y-2*a^2*A*B^5*x*y;
t_newton := 4.202;
iters := 78;
r_min := 0;
vec_d := Vector(76, [1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1/16,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1/28,0,0,7/
27,0,0,0,0,0,0,0,0,0,0,0,0]);
q_list := Vector(76, [4*a^4*B^2*x^2+a^6*A^2+a^4*y^2-2*a*B*x+4*a*A*y-4*a^4-4*a^
2+a^5*A^2*B*x+a^3*B*x*y^2+a^3*A*x^2*y+a*A*B^2*y+4*a^5*B*x-a^5*A^3*y-2*a^5*A*y+
2*a^4*A^2*y^2+a^2*x^2-a^3*B*x^3-a^3*A*y^3+2*a^2*B^2*x^2+4*a^2*A^2*y^2-a*B^3*x+
B^2-6*a^2*A*B*x*y-6*a^4*A*B*x*y,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
-4*a^2*A*B^2-4*a^3*x^2*y+16*a^2*A*y^2+4*a^4*A*y^2+16*a^3*A^2*y-24*a^2*B*x*y-24
*a^3*A*B*x-16*a^4*B*x*y-16*a^5*A*B*x+16*a*y+16*a^2*A+16*a^3*y+16*a^4*A-4*a^3*y
^3-4*a^6*A^3-4*a*B^2*y-4*a^4*A*x^2+4*a^5*A^2*y,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
,0,4*a*A*B*y+4*a^2*A*x*y+6*a^3*A*B*y+6*a^4*A*x*y-4*a^2*B-4*a^3*x-4*a^4*B+a^3*x
^3-4*a^5*x-a*B^2*x+a^2*B*y^2-a^2*B*x^2+a^3*x*y^2-4*a^3*B^2*x-4*a^4*B*x^2+a^4*A
^2*B+a^5*A^2*x+B^3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2*a^2*B^2*x*y+2*a^2*A*B*
y^2-2*a^3*B*x^2*y+2*a^3*A*x*y^2-2*a^3*A*B^2*x+2*a^3*A^2*B*y-2*a^4*A*B*x^2+2*a^
4*A^2*x*y+28*a*B*y+28*a^2*x*y+28*a^2*A*B+26*a^3*B*y+28*a^3*A*x+26*a^4*x*y+26*a
^4*A*B+26*a^5*A*x,0,0,27/7*a^2*B^2*x*y-27/7*a^2*A*B*y^2+27/7*a^3*B*x^2*y-27/7*
a^3*A*x*y^2+27/7*a^3*A*B^2*x-27/7*a^3*A^2*B*y+27/7*a^4*A*B*x^2-27/7*a^4*A^2*x*
y+27/7*a^3*B*y+27/7*a^4*x*y+27/7*a^4*A*B+27/7*a^5*A*x,0,0,0,0,0,0,0,0,0,0,0,0]
);