Apprendre Maple

[gasy72]

voici un petit procedure pour calculer le pgcd de deux polynomes:

ordonner les polynômes::

ordonner_poly:=proc(poly1,poly2)
if degree(poly2) > degree(poly1) then poly1, poly2
else poly2,poly1
fi
end;

algorithme d'euclide::

algo_euclide_poly:=proc(petitdegre,granddegre)
local
reste,liste_etapes,modulo,divise;
reste:=petitdegre;
modulo:=granddegre;
liste_etapes:=[[reste,modulo]]:
while reste<>0 do
divise:=modulo:
modulo:=reste:
reste:=prem(divise,modulo,x):
liste_etapes :=[op(liste_etapes), [reste,modulo]]
od;
array(liste_etapes),modulo:
end:
procedure de pgcd::

pgcd_poly:=proc(poly1,poly2)
local poly_non_unit;
poly_non_unit := algo_euclide_poly(ordonner_poly(poly1, poly2))[2];
sort(poly_non_unit/lcoeff(poly_non_unit))
end;

[gasy72]

>x:='x';

p1:=(x-4)*(x+2)*(x-1);

p2:=(x-4)*(x-3)*(x-2)*(x-1)^2;


algo_euclide_poly(p1,p2);

donne:



matrix([[(x-4)*(x+2)*(x-1), (x-4)*(x-3)*(x-2)*(x-1)^2], [-240-60*x^2+300*x, (x-4)*(x+2)*(x-1)], [0, -240-60*x^2+300*x]]), -240-60*x^2+300*x

> pgcd_poly(x^4+2*x^3+3*x^2+2*x+1,4*x^3+6*x^2+6*x+2);
2
x + x + 1

[Henriette Dumans]

Merci, gasy72 c justement ce que je cherchais
A+

[prof_simplet]

Bonjour, vous savez aussi qu'il existe la commande gcd de Maple qui produit la mem chose.
Bon week-end.