Apprendre Maple

[mbm]

je remarque que mon procedure Gauss retourne des resultats errones et je essier des les resoudre vous pouver m'aidez s'il vous plais?
voiala le procedure:
*************************************************************

> Euclide:=proc(a,b)
local r , u,v, q, i, u0, u1, v0, v1, r0, r1;
r0:=a; r1:=b; u0:=1; u1:=0; v0:=0; v1:=1;
while r1<>0 do q:=iquo(r0,r1);
r:=r0-r1*q; u:=u0-u1*q; v:=v0-v1*q;
r0:=r1; r1:=r; u0:=u1; u1:=u; v0:=v1; v1:=v; end do; [r0,u0,v0]
end proc:
> gauss:=proc(A,p)
local n,i,j,piv,s,d,l,B,x,ipiv,k, eu,c,sgn;
n:=linalg[coldim](A); d:=1; B:=map(x->x mod p, evalm(A));sgn:=1;
for i to n do j:=i; piv:=B[j,i] mod p;
while piv=0 and j<n do j:=j+1; piv:=B[j,i] mod p end do ;
if piv=0 then return(0) else d:=d*piv;
for k from i to n do c:=B[i,k]; B[i,k]:=B[j,k]; B[j,k]:=c;
sgn:=-1*sgn; end do;
eu:=Euclide(piv,p);
if piv < 0 then do
ipiv:=-eu[2] mod p; end do;
else
ipiv:=eu[2] mod p; end if;
for k from j+1 to n do l:=B[k,i]*ipiv mod p;
for s from i to n do B[k,s]:=B[k,s]-l*B[i,s] end do ; end do;
end if ; end do; sgn*d mod p ; end proc:
>

[prof_simplet]

Bonjour,
Déjà vous devriez expliquer ce que vous cherchiez à faire par ces procédures.
Si j'ai bien compris, la première , Euclide, retourne [r0,u0,v0] où r0=pgcd(a,b) et (u0,v0) est une solution de l'équation ax+by=pgcd(a,b).
Si a et b ne sont pas premiers enre eux, on obtient u0=0 v0=1 ou u0=1 v0=0. Donc tout semble correct dans Euclide.
Pour la deuxième, pourriez-vous expliquer votre démarche?
La fonction préexistante de Maple
Gausselim(A) mod p
doit faire le même travail, non?

[mbm]

je veut que gauss me calcule le determinant d'une matice modulo des entiers premieres.
j'espere que vous pouvez m'aidez .

[ALS]

Bonjour, j'ai modifié quelque peu la seconde procédure en normalisant la ligne concernée lorsque piv<>0: on multiplie toute la ligne par 1/piv mod p de façon à faire apparaitre un 1 en diagonale.
Ainsi la matrice B obtenue est triangulaire supérieure, son déterminant modulo p est le produit des termes diagonaux.
La procédure donne pour résultat la liste [B,det(B]] c'est à dire [evalm(B),product(B[h,h],h=1..n)]
On peut vérifier la matrice B obtenue avec le résultat donné par la procédure préexistante Gausselim(A) mod p; de Maple.
Ca semble maintenant correctement fonctionner: le vérifier plus en détail sur des matrices générées aléatoirement:
> A:=linalg[randmatrix](4,4); u:=rand(1..20): p:=nextprime(u()); gauss(A,p); Gausselim(A) mod p;

[mbm]

merci de la reponse se bien marche