% Matrices for 18 killer whale pods in British Columbia;
% From Appendix of Brault and Caswell (1993);
% Format of matrices for each pod = [0 F2 F3 0; G1 P2 0 0; 0 G2 P3 0; 0 0 G3 P4];

J01=[0 0.0067 0.1632 0; 0.9535 0.8827 0 0; 0 0.0802 0.9586 0; 0 0 0.0414 0.9752];
K01=[0 0.0062 0.1737 0; 1.0000 0.9020 0 0; 0 0.0694 0.9582 0; 0 0 0.0418 0.9855];
L01=[0 0.0037 0.0988 0; 0.9562 0.9030 0 0; 0 0.0722 0.9530 0; 0 0 0.0406 0.9798];
A01=[0 0.0043 0.1148 0; 1.0000 0.9015 0 0; 0 0.0727 0.9515 0; 0 0 0.0485 0.9667];
A04=[0 0.0042 0.1054 0; 0.8165 0.8903 0 0; 0 0.0774 0.9515 0; 0 0 0.0485 0.9810];
A05=[0 0.0027 0.0732 0; 1.0000 0.9123 0 0; 0 0.0730 0.9515 0; 0 0 0.0485 0.9545];
B01=[0 0.0025 0.0651 0; 1.0000 0.9254 0 0; 0 0.0746 0.9515 0; 0 0 0.0485 0.9810];
C01=[0 0.0047 0.1159 0; 1.0000 0.9200 0 0; 0 0.0800 0.9706 0; 0 0 0.0294 0.9608];
D01=[0 0.0068 0.1761 0; 1.0000 0.9241 0 0; 0 0.0759 0.9562 0; 0 0 0.0438 1.0000];
G01=[0 0.0061 0.1418 0; 1.0000 0.9167 0 0; 0 0.0833 0.9286 0; 0 0 0.0714 1.0000];
G12=[0 0.0050 0.1251 0; 1.0000 0.9216 0 0; 0 0.0784 0.9515 0; 0 0 0.0485 0.9810];
H01=[0 0.0021 0.0542 0; 1.0000 0.9254 0 0; 0 0.0746 0.9515 0; 0 0 0.0485 0.9810];
I01=[0 0.0027 0.0732 0; 1.0000 0.9286 0 0; 0 0.0714 0.9515 0; 0 0 0.0485 0.9810];
I02=[0 0.0045 0.1220 0; 1.0000 0.9286 0 0; 0 0.0714 0.9515 0; 0 0 0.0485 1.0000];
I11=[0 0.0052 0.1428 0; 1.0000 0.9286 0 0; 0 0.0714 0.9515 0; 0 0 0.0485 0.9810];
I18=[0 0.0037 0.0998 0; 1.0000 0.9286 0 0; 0 0.0714 0.9515 0; 0 0 0.0485 0.9810];
I31=[0 0.0047 0.1273 0; 1.0000 0.9286 0 0; 0 0.0714 0.9515 0; 0 0 0.0485 0.9810];
R01=[0 0.0024 0.0797 0; 1.0000 0.8929 0 0; 0 0.0595 0.9515 0; 0 0 0.0485 1.0000];

% Sensitivity matrix for a mean matrix of all the matrices above;
S = [0.0422 0.3608 0.3686 0.3699;
     0.0443 0.3785 0.3867 0.3881;
     0.0663 0.5670 0.5793 0.5813;
     0      0      0      0     ];