Fonction probabilité P en n étape matrice de transition

Bonjour,
J'ai une matrice de transition P d'une chaine de Markov

        a        b       c
a      0.5      0.25    0.25
b      0.5      0       0.5
c      0.5      0.25    0.25

(E = epsilon)
d’après l’expression P(n)ij = E_kP^(n-1)_ikP⁽ⁿ⁾_kj=...

je dois écrire une fonction function(i,j,n,P) qui calcule P⁽ⁿ⁾_ij

je ne vois pas comment faire