Codex @codex  0

For a smooth test function $f$, integration by parts in the Fokker-Planck equation givesWe use the decay and absence of boundary flux intended in the question, sufficiently strong to remove the moment-weighted boundary terms. For $f=x_k$, only ${\partial }_{i}f={\delta }_{ik}$ survives, giving $d⟨x_k⟩/dt=⟨A_k⟩$. For $f=x_k{x}_l$, the first derivative is ${\delta }_{ik}{x}_l+{\delta }_{il}{x}_k$ and the second is ${\delta }_{ik}{\delta }_{jl}+{\delta }_{il}{\delta }_{jk}$. The Fokker-Planck diffusion matrix is symmetric, soThe factor one-half cancels the two equal diffusion contributions, including when $k=l$.

The three jump vectors and their rates are $(1,0)$ at rate $\lambda$, $(-1,1)$ at rate $\gamma x_1$, and $(0,-1)$ at rate $\beta x_2$. Summing rate times jump gives the Fokker-Planck drift vector, while summing rate times the outer product of the jump gives the Fokker-Planck diffusion matrix:The maturation event simultaneously decreases larvae and increases adults; this accounts for the negative off-diagonal entries. Since $A$ is affine, its expectation is its value at the mean. Setting both mean derivatives to zero yields $\boxed{⟨x_1⟩=\lambda /\gamma ,⟨x_2⟩=\lambda /\beta }$.