I'm trying to solve equation in the attached pdf, which describes anistropic diffusion in 3D with an additional term to account for hydrogen bonding and unbonding of the diffusing substance to the medium. I've considered Laplace transforms, then solving in the Laplace domain, then inverting numerically. Ideally I would get a fully analytical solution, but I'm not sure how to approach it. Solving it numerically in three dimensions (using Crank-Nicolson or something similar) would undoubtedly be very computationally expensive, correct? If anyone has any suggestions on how to arrive at an analytical solution to this, or an efficient way to implement it numerically, I would greatly appreciate your input. I’ve been struggling with this for some time now, and just haven’t gotten very far. Thanks in advance.