This paper examines the dynamics of drug concentration in the body accounting for random factors like patient and environmental variability. We develop an explicit solution for drug concentration using a Stochastic Differential Equation (SDE) model. We calculate formulas for the expected value and variance, enabling statistical evaluation and prediction of the drug’s concentration trajectory and its uncertainty. The unknown parameters in the model are estimated using the method of moments. We apply our proposed methods to a real-world dataset, providing useful insights analysis of drug concentration and the determination of its therapeutic range.