P:=matrix(2,2, [0.5, 0.5,0.25,.75]);

Pn := proc(n::nonnegint)

if n<2 then

P;

else

evalm(Pn(1)&*Pn(n-1));

fi;

end:

seq( Pn(i), i=1..15);