A mathematical model of Wu et al. [J. Membr. Sci 254 (2005) 119-127] of a cationic glucose-sensitive membrane is discussed. The model involves the system of non-linear steady-state reaction-diffusion equations. Analytical expres-sions pertaining to concentration of oxygen, glucose, and gluconic acid for all values of parameters are presented. We have employed Homotopy analysis method to evaluate the approximate analytical solutions of the non-linear boundary value problem. A comparison of the analytical approximation and numerical simulation is also presented. A good agreement between theoretical predictions and numerical results is observed.