test de primalité de fermat et nombres de carmicael sous SCILAB

Citation :
// TEST DE PRIMALITE DE FERMAT
// pourquoi ce test renvoie false su les nombres de Carmichael?

function [test] = primfermat(n)
// en argument un entier supérieur à 7
// renvoie true si composé, false si probable premier

  if ~((n>7)&(ceil(n)==n)); disp("argument mal défini"); end
  m = 2:7  //les témoins sont 2, 3, ..., 7
  p=puissmatmodn(m,n-1,n);  //on calcule les puissances n-1 de chaque élément de m modulo n
  if p == ones(1,6)  //si tout est égal à 1 dans le calcul précédent, n est probable premier
    test = %T
  else
    test = %F
  end
endfunction


// FONCTION EVITANT LES GRANDS NOMBRES DANS LE CALCUL DES PUISSANCES D'UNE MATRICE
function [y] = puissmatmodn(mat,p,n)
// en argument un vecteur, la puissance, et le modulo
  y=mat
  for i=1: (p-1) // c'est ":" puis "(" pas ^^'''
    y=modulo(y.*mat,n)
  end
endfunction