O reator nuclear natural de Oklo foi alvo de pesquisas por vários cientistas ao longo de décadas, devido à sua condição de criticalidade. Admitindo-se que esse reator possa ser tratado com dimensão infinita a dois grupos de energia na condição crítica, a equação da difusão de nêutrons multigrupo é dada por
!$ \Bigl ( \textstyle \sum_{a1} + \textstyle \sum_{s}^{1 \rightarrow 2} - \nu_1 \textstyle \sum_{f1} \Bigr ) \phi_1 - \nu_2 \textstyle \sum_{f2} \phi_2 = 0 !$
!$ \textstyle \sum_{s}^{1 \rightarrow 2} \phi_1 - \textstyle \sum_{a2} \phi_2 = 0 !$
A equação acima admite a seguinte representação matricial: !$ M \varPhi = 0 !$, onde !$ \varPhi = \begin{pmatrix} \phi_1 \\ \phi_2 \end{pmatrix} !$
Nessa condição, a matriz !$ M !$ é representada por
