function q=jeu(p)

// On reserve de l'espace memoire
//stacksize(1000000)

// Affection de 1-p a d pour alleger l'ecriture
d = 1-p

// Matrice identité
I = eye(18,18)

// Vecteur pour recuperer l'état le joueur 1 gagne en partant du premier etat
E = zeros(1,18)
E(1,1) = 1

// Creation de la matrice M0
M0 = zeros(18,18)

// Creation de la matrice M2
M1 = zeros(18,1)

// Creation de la matrice M3
M2 = zeros(18,1)

// On rempli a la main la matrice M0
M0(1,2) = p
M0(1,3) = d

M0(2,4) = p
M0(2,5) = d

M0(3,5) = p
M0(3,6) = d

M0(4,7) = p
M0(4,8) = d

M0(5,8) = p
M0(5,9) = d

M0(6,9) = p
M0(6,10) = d

M0(7,11) = d

M0(8,11) = p
M0(8,12) = d

M0(9,12) = p
M0(9,13) = d

M0(10,13) = p

M0(11,14) = d

M0(12,14) = p
M0(12,15) = d

M0(13,15) = p

M0(14,16) = d

M0(15,16) = p

M0(16,17) = p
M0(16,18) = d

M0(17,16) = d

M0(18,16) = p


// On rempli la matrice M1
M1(4,1) = p
M1(11,1) = p
M1(14,1) = p
M1(17,1) = p

// On rempli la matrice M2
M2(10,1) = d
M2(13,1) = d
M2(15,1) = d
M2(18,1) = d

// On calcul maintenant X1, probabilite que le joueur A gagne le jeu
X1 = inv(I-M0)*M1
q = E*X1


endfunction