“…To approximate µ i,j,k (i = 1, 2), we use the methods described in [40] to express µ 1,j,k as the following ordinary differential equation that can be solved numerically using the Matlab function ode113: A similar method is used to approximate µ 2,j,k .…”