Вот именно, я не понимаю, что за числа я нашел.
И что с ними делать, тоже не понимаю.
Разве это не Вы писали:
вероятности нахождения системы в различных состояниях
...
после первого выстрела соответствующие вероятности будут равны первой строке матрицы ...
Найдите определение, что у Вас обозначено для цепи Маркова через

, которое Вы вычисляете по выписанной выше формуле.
Как мне представляется, вероятности

- относятся к первому выстрелу,

- ко второму, а

к третьему. Но тогда не понимаю, зачем еще две матрицы, вот в чем беда.
Разве при вычислении

Вы не использовали вторую матрицу? А при вычислении

- третью?
Давайте так. Вот был вектор начальных вероятностей системе находиться в её состояниях. Умножаете его на матрицу перехода из начального состояния в первое. Что за вектор, из вероятностей каких событий получаете? Дальше, умножаете его на матрицу вероятностей перехода из первого состояния во второе - что за вектор получаете, из вероятностей системе
что делать? Дальше умножаете его на матрицу вероятностей перехода из второго состояния в третье - что за новый вектор получаете?
Надо бы с базовыми понятиями-то разобраться.