# Cercle tangent à deux paraboles
R.<a,r,x> = PolynomialRing(QQ,"a,r,x", order="lex")
CerPar1 = (x-a)^2+(x^2)^2-r^2
CerPar2 = (x-a)^2+(x^2-1/2)^2-r^2
Tan1= CerPar1.polynomial(x).discriminant()
Tan2= CerPar2.polynomial(x).discriminant()
S=R.remove_var(x)
I=S.ideal([S(Tan1),S(Tan2)])
P=I.elimination_ideal(S(a)).gens()[0]
Q=factor(P)[1][0]
print("Solutions réelles du système :\n",I.variety(AA))
print("Polynôme Q minimal de r :",Q)
rsols=solve(SR(Q),SR(r))
print("Racines de Q :\n",rsols)
rsolsnum=[rs.rhs().n() for rs in rsols]
print("Valeurs numériques approchées des racines :\n",rsolsnum)
Solutions réelles du système :
[{r: -0.1781313487737803?, a: -0.4843785328691515?}, {r: -0.1781313487737803?, a: 0.4843785328691515?}, {r: 0.1781313487737803?, a: -0.4843785328691515?}, {r: 0.1781313487737803?, a: 0.4843785328691515?}]
Polynôme Q minimal de r : 5038848*r^8 - 1306368*r^6 + 204704*r^4 - 6192*r^2 + 27
Racines de Q :
[
r == -1/54*sqrt(3)*sqrt(26*sqrt(3) + 486*sqrt(-152/6561*sqrt(3) - 776/19683) + 63),
r == 1/54*sqrt(3)*sqrt(26*sqrt(3) + 486*sqrt(-152/6561*sqrt(3) - 776/19683) + 63),
r == -1/54*sqrt(3)*sqrt(26*sqrt(3) - 486*sqrt(-152/6561*sqrt(3) - 776/19683) + 63),
r == 1/54*sqrt(3)*sqrt(26*sqrt(3) - 486*sqrt(-152/6561*sqrt(3) - 776/19683) + 63),
r == -1/54*sqrt(3)*sqrt(12*sqrt(2/3)*sqrt(57*sqrt(3) - 97) - 26*sqrt(3) + 63),
r == 1/54*sqrt(3)*sqrt(12*sqrt(2/3)*sqrt(57*sqrt(3) - 97) - 26*sqrt(3) + 63),
r == -1/162*sqrt(3)*sqrt(-108*sqrt(2/3)*sqrt(57*sqrt(3) - 97) - 234*sqrt(3) + 567),
r == 1/162*sqrt(3)*sqrt(-108*sqrt(2/3)*sqrt(57*sqrt(3) - 97) - 234*sqrt(3) + 567)
]
Valeurs numériques approchées des racines :
[-0.381250423062243 - 0.184949974623226*I, 0.381250423062243 + 0.184949974623226*I, -0.381250423062242 + 0.184949974623226*I, 0.381250423062242 - 0.184949974623226*I, -0.178131348773780, 0.178131348773780, -0.0723719413811676, 0.0723719413811676]
Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 3 invités
Tu pars déja ?
Identification
Pas encore inscrit ?
Ou identifiez-vous :