C1=circonferenza con diametro uguale 
all'asse maggiore dell'ellisse

C1

C2=circonferenza con diametro uguale
all'asse minore dell'ellisse

C2

A = (5.06, 0.22)
c = Circle(A, 7)
d = Circle(A, 4)
B = (5.06, 7.398790936919672)
a = Segment(A, B)
C = (5.06, -7.056040940758926)
b = Segment(A, C)
D = (13.263138919048, 0.22)
e = Segment(A, D)
E = (-2.6646762016281347, 0.22)
f = Segment(A, E)
F = Intersect(c, a, 1)
G_2 = Intersect(d, a, 1)
H_2 = Intersect(d, e)
I = Intersect(c, e, 1)
J = Intersect(d, b, 1)
K = Intersect(d, f, 1)
L_2 = Intersect(c, f, 1)
M_2 = (10.977596277205802, 6.137596277205802)
g = Segment(A, M_2)
N_2 = (8.9236, 3.00632)
h = Circle(G_2, N_2)
k = Circle(H_2, Radius[h])
O = Point(k)
P_2 = Point(k)
p = CircleArc(H_2, O, P_2)
Q = (8.57754, 5.66832)
q = CircleArc(G_2, N_2, Q)
R = (-0.5751388832042146, -5.415138883204215)
i = Segment(A, R)
T = (-0.8714570818644587, 6.151457081864459)
j = Segment(A, T)
U = (10.948991234498486, -5.668991234498485)
l = Segment(A, U)
S = Intersect(c, b, 1)
N = Intersect(d, g, 1)
W = (7.334385999999997, 5.314274000000002)
r = Circle(G_2, W)
s = Circle(N, Radius[r])
Z = Intersect(r, s, 1)
A_1 = Point(r)
B_1 = (6.192387999999997, 6.251298000000002)
t = CircleArc(G_2, A_1, B_1)
C_1 = Point(s)
D_1 = Point(s)
c_1 = CircleArc(N, C_1, D_1)
m, n = AngularBisector(a, j)
a_1, b_1 = AngularBisector(j, f)
M = Intersect(d, b_1, 1)
G = Intersect(c, b_1, 1)
P = Intersect(d, n, 1)
L = Intersect(c, n, 1)
I_1 = (8.272217071341974, 3.915518130045147)
d_1 = Segment(M, I_1)
J_1 = (7.738784026555628, 3.6802183262354555)
e_1 = Segment(G, J_1)
H = Intersect(c, g, 1)
L_1 = (10.48836004835104, 3.0484271247461914)
f_1 = Segment(N, L_1)
M_1 = (10.009747468305832, 2.754760134873163)
g_1 = Segment(H, M_1)
N_1 = (12.18292470609304, 1.7507337294603589)
h_1 = Segment(P, N_1)
O_1 = (11.527156727579007, 1.1571469106800585)
i_1 = Segment(L, O_1)
P_1 = Intersect(d_1, e_1)
Q_1 = Intersect(f_1, g_1)
R_1 = Intersect(h_1, i_1)
k_1 = Circle(G_2, Segment[L_2, A])
F1 = Intersect(k_1, f, 1)
F2 = Intersect(k_1, e, 1)
U_1 = Intersect(d, a_1, 1)
V_1 = Intersect(c, a_1, 1)
W_1 = (12.854168426382618, -1.3107337294603587)
j_1 = Segment(U_1, W_1)
Z_1 = (11.527156727579007, -0.49193152514941896)
l_1 = Segment(V_1, Z_1)
A_2 = Intersect(j_1, l_1)
V = Intersect(d, l, 1)
E_1 = Intersect(c, l, 1)
F_1 = (10.698640852673268, -2.589921231707026)
m_1 = Segment(V, F_1)
G_1 = (9.998543322027457, -2.228580570728543)
n_1 = Segment(E_1, G_1)
H_1 = Intersect(m_1, n_1)
K_1 = Intersect(d, m, 1)
S_1 = Intersect(c, m, 1)
T_1 = (8.34408251780129, -3.475518130045147)
p_1 = Segment(K_1, T_1)
B_2 = (7.738784026555628, -3.1485327496558426)
q_1 = Segment(S_1, B_2)
C_2 = Intersect(p_1, q_1)
r_1 = Ellipse(F1, F2, G_2)