High order accurate, one-sided finite-difference approximations to concentration gradients at the boundaries, for the simulation of electrochemical reaction-diffusion problems in one-dimensional space geometry.