mopa03r1.m

function A = mopa03r1
%mopa03r1   A 4-by-4 matrix from \cite{mopa03}.

Thalf = [3.8235*60*24 3.10 26.8 19.9]/60;  % Half lives in seconds/
a = log(2)./Thalf;  % decay constant
A = diag(-a) + diag(a(1:end-1),-1);