Algorithme de Faddeev-Leverrier

Bonjour j'aimerais savoir si quelqu'un saurait coder ça en langage algorithmique ou fortran 90 svp ?

http://fr.wikipedia.org/wiki/Algorithme_de_Faddeev-Leverrier

Salut cet valgo est simple si tu as défini auparavant ,la trace,la multiplication et l'addition des matrices carrées:

ALGO (du à souriau).
entrée matrice A,sortie déterminant,adjointe(en référence).

en C++ cela donne :

double souriau(const matrice& A,matrice& B)
{
   matrice B=A;
   double tr=trace(A);
   B=B-tr*I   (I matrice unité);
  for (i=2;i  {
      B=B*A;
      tr=-trace(B)/i;
      B=B+tr;
  }
  tr=0;
  for (i=0;i     tr=tr+B[0][i]*A[i][0];
  double detA=-tr;
  return detA;
}

pour info  B est au signe près l'adjointe de A;
adj(A)=(-1)^(dimA-1)*B