# La liste des premiers nombres de Fermat premiers
fermat:=proc()
local l,n,k;

 l := [NULL];
 k := 1;
 n:=2^(2^k)+1;

 while (isprime(n)) do
   l := [op(l), n];
   k := k+1;
   n := 2^(2^k)+1;

 od;

 l;
end;

# L'algorithme de Syracusse
syracusse:=proc(nombre)
local k, l, n;

  n := nombre;
  l := [NULL];
  for k while(n<>1) do
    if type(n,odd) then
       n := 3*n+1;
    else
       n := n/2;
    fi;
    l := [op(l), n];
  od;

 l;
end;

# Approximation d'une racine de polynome par m&eacute;thode de Dichotomie
dichotomie := proc(f,a,b,p)
local a1,b1;

  a1 := evalf(a);
  b1 := evalf(b);
  while abs(b1-a1) > p do
    if f(a1)*f(b1)<0 then
      b1 := 0.5*(a1+b1);
    else
      a1 := 0.5*(a1+b1);
    fi;
  od;

 0.5*(a1+b1);
end;