%
%  Precip data: Dry/Wet
%

clear;
rand('seed',0);

P = zeros(2,2); 
S = 2;
Pww = 643/771;
Pdd = 186/309;
P(1,1) = Pdd;
P(1,2) = 1 - Pdd;
P(2,2) = Pww;
P(2,1) = 1- Pww; 

%% Finding stationary distribution 

A = P' - eye(S);
A(S,:) = ones(1,S);
b = zeros(S,1);
b(S) = 1;
pi = (A\b)';