Voici la fonction d’exponentiation modulaire rapide que j'ai faite en scilab.
Elle est en mode itératif, je n'arrive pas à faire la même fonction en mode récursif, d'autant plus qu'il n'y a pas d'opérateur & ou >> comme en langage C pour manipuler un nombre en binaire (bit-à-bit)
d'où l'utilisation de la variable tableau n!!!!

e = 13;
d = 77;
u = [1, 0, 1, 1]; // Tableau avec les données 13 en binaire
m = 9;

function [mes_signe] = ExpModRapide (a, n, m)
N= length(n)
mes_signe = 1
for i = 1 : N
if (n(i) == 1) then
mes_signe = modulo (mes_signe * a, m)
end
a = modulo((a*a),m)
end
endfunction

[mes_signe] = ExpModRapide (m, u,d)