powerrusse(n, p);

             function y = powerrusse(n, p);
             % calcul de n^p.
             
             if p == 0
               y = 1;
             elseif rem(p,2) == 1                 % p est-il impair ?
               y  = n*powerrusse(n*n, floor(p/2));% floor(a) renvoie le plus grand   
             else                                 % entier infrieur ou gal  a
               y  = powerrusse(n*n, floor(p/2));
             end;